FREAD1 error while generating lwdflux_096a.a

[Jeongin Kim]

Hi Alan, 

I was able to successfully run the model for 2 years.
I'm now running a simulation starting from January 1st, 1996, but an error appeared while generating lwdflux_096.a.When I looked through the log file of the other forcing files, it looks like the FDATE is different.
Other forcing files generate files at FDATE=1.208/1996, while lwdflux_096.a generates a line at FDATE = 1.207/1996.

 

 KPAR: MIN= 126.06987762 MAX= 462.87863159 AVE= 337.10989380 RMS= 343.81692505
          WRITING KPAR RECORD    1    FDAY = 34699.000  FDATE =   1.000/1996

 KPAR: MIN= 125.06996155 MAX= 442.03073120 AVE= 320.16439819 RMS= 324.78378296
          WRITING KPAR RECORD    2    FDAY = 34699.042  FDATE =   1.042/1996

 KPAR: MIN= 124.42317200 MAX= 466.64166260 AVE= 337.39642334 RMS= 344.18167114
          WRITING KPAR RECORD    3    FDAY = 34699.083  FDATE =   1.083/1996

 KPAR: MIN= 309.44931030 MAX= 497.95101929 AVE= 352.56890869 RMS= 353.56066895
          WRITING KPAR RECORD    4    FDAY = 34699.125  FDATE =   1.125/1996

 KPAR: MIN= 121.61235809 MAX= 501.66824341 AVE= 338.14874268 RMS= 344.93386841
          WRITING KPAR RECORD    5    FDAY = 34699.167  FDATE =   1.167/1996

 KPAR: MIN= 120.99027252 MAX= 503.64739990 AVE= 337.64880371 RMS= 344.47372437
          WRITING KPAR RECORD    6    FDAY = 34699.207  FDATE =   1.207/1996



                    *****  ERROR IN FREAD1  -  INPUT FLUX DAY IS 34699.25, BUT SHOULD BE     0.00  *****

And it generates a weird dlwflx line in the log file. 

 12491641 (1996/001 00) region-wide mean Density Dev:        2.6311533016
  airprs: day,span,range = 34699.08333  0.041667  -3.8452048E+03   3.1462170E+03
  dlwflx: day,span,range = 34699.08333  0.041667   1.2442317E+02   4.6664166E+02
  dswflx: day,span,range = 34699.08333  0.041667   0.0000000E+00   4.1666672E-02
 tau_ewd: day,span,range = 34699.12500  0.041667  -1.4571665E+01   2.4243092E+01
 tau_nwd: day,span,range = 34699.12500  0.041667  -2.0041187E+01   1.6874462E+01
  airtmp: day,span,range = 34699.16667  0.041667  -3.1170046E+01   2.8475142E+01
  vapmix: day,span,range = 34699.16667  0.041667   1.9856653E-04   1.8959370E-02
  precip: day,span,range = 34699.16667  0.041667   0.0000000E+00   1.6513003E-06
  airprs: day,span,range = 34699.20833  0.041667  -3.8564065E+03   3.2238669E+03
  dlwflx: day,span,range = 34699.20703**********   1.2099027E+02   5.0364740E+02
forrtl: severe (64): input conversion error, unit -5, file Internal List-Directed Read
Image              PC                Routine            Line        Source
hycom              00000000007140BB  Unknown               Unknown  Unknown
hycom              00000000007136EB  Unknown               Unknown  Unknown
hycom              00000000005E3BC9  rdpall1_                 3131  forfun.F90
hycom              00000000005E2B84  rdpall_                  2994  forfun.F90
hycom              00000000005DBCCE  forfunh_                 1194  forfun.F90
hycom              00000000005827C3  mod_hycom_mp_hyco        2382  mod_hycom.F90
hycom              0000000000414D35  MAIN__                    161  hycom.F90
hycom              0000000000414CBD  Unknown               Unknown  Unknown
libc.so.6          00001555536295D0  Unknown               Unknown  Unknown
libc.so.6          0000155553629680  __libc_start_main     Unknown  Unknown
hycom              0000000000414BD5  Unknown               Unknown  Unknown
forrtl: severe (64): input conversion error, unit -5, file Internal List-Directed Read
Image              PC                Routine            Line        Source
hycom              00000000007140BB  Unknown               Unknown  Unknown
hycom              00000000007136EB  Unknown               Unknown  Unknown
hycom              00000000005E3BC9  rdpall1_                 3131  forfun.F90
hycom              00000000005E2B84  rdpall_                  2994  forfun.F90
hycom              00000000005DBCCE  forfunh_                 1194  forfun.F90
hycom              00000000005827C3  mod_hycom_mp_hyco        2382  mod_hycom.F90
hycom              0000000000414D35  MAIN__                    161  hycom.F90
                                                                                   

Would this happen because 1996 is a leap year?
Is there be anything I'm missing?

Thank you,
Jeongin

[Alan Wallcraft]

The cfsr-sec_1996_01hr_dlwsfc.nc file has the wrong MT time values.

Compared to dswsf:

narwhal01 161> diff cfsr-sec_1996_01hr_dlwsfc.cdl cfsr-sec_1996_01hr_dswsfc.cdl
1c1
< netcdf cfsr-sec_1996_01hr_dlwsfc {
---
> netcdf cfsr-sec_1996_01hr_dswsfc {

40c40
< :title = "CFSR 0.313g 1996 01hrly sec-CERES dlwsfc" ;
---
> :title = "CFSR 0.313g 1996 01hrly sec-CERES dswsfc" ;
47,48c47,48
<     34699.2083333333, 34699.25, 0, 34699.3333333333, 34699.375,
<     34699.4166666667, 34699.4583333333, 34699.5, 34699.5416666667,
---
>     34699.2083333333, 34699.25, 34699.2916666667, 34699.3333333333,
>     34699.375, 34699.4166666667, 34699.4583333333, 34699.5, 34699.5416666667,

I will get this fixed, but in the mean time you can fix you local copy using NCO:

ncap2 -O -s 'MT(7)=34699.2916666667' cfsr-sec_1996_01hr_dlwsfc.nc cfsr-sec_1996_01hr_dlwsfc_newMT.nc

I have not tested that this file works.  Let me know either way, since I won't replace the existing file until this is confirmed to work.

Alan.

[Jeongin Kim]

Hi Alan,

Revising MT time worked!

Thank you, 

Jeongin

2026년 2월 9일 월요일 AM 10시 39분 41초 UTC-5에 Alan Wallcraft님이 작성: