about reading surface forcing file

[wrh]

I am sure that my surface forcing variable sustr have minimum value -1.731250825225399 N.m-2 and maximum value 0.753999437586370 N.m-2,but why did roms read the sustr and present different min and max values ?? svstr , surface net heat flux,surface net freswater flux also have the same problems
GET_2DFLD - surface u-momentum stress, t = 2 06:00:00
(Rec=0000006, Index=2, File: cfsr_smsu_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.7500)
(Min = -4.63317233E-05 Max = 3.56971141E-05)
GET_2DFLD - surface v-momentum stress, t = 2 06:00:00
(Rec=0000006, Index=2, File: cfsr_smsv_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.7500)
(Min = -6.59743110E-05 Max = 6.19260676E-07)
GET_2DFLD - surface air pressure, t = 2 06:00:00
(Rec=0000006, Index=2, File: cfsr_Pair_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.7500)
(Min = 1.02890894E+03 Max = 1.02998242E+03)
GET_2DFLD - solar shortwave radiation flux, t = 2 06:00:00
(Rec=0000006, Index=2, File: cfsr_swrad_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.7500)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface net heat flux, t = 2 06:00:00
(Rec=0000006, Index=2, File: cfsr_shflux_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.7500)
(Min = -1.27329486E-04 Max = -4.98705482E-06)
GET_2DFLD - sea surface temperature climatology, t = 3 00:00:00
(Rec=0000003, Index=1, File: cfsr_SST_2015_roms.nc)
(Tmin= 1.0000 Tmax= 540.0000)
(Min = -8.60936344E-01 Max = 1.68873787E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 3 00:00:00
(Rec=0000003, Index=1, File: cfsr_dqdsst_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.0000)
(Min = -6.21000215E-06 Max = -8.72371429E-07)
GET_2DFLD - surface net freswater flux, (E-P), t = 2 06:00:00
(Rec=0000006, Index=2, File: cfsr_swflux_2015_roms.nc)
(Tmin= 1.0000 Tmax= 425.7500)

[kate]

ROMS scales these things internally because that's how it uses them. There's this code in inp_par.F:

Code: Select all

        cff=1.0_r8/rho0
        Fscale(idUsms,ng)=cff*Fscale(idUsms,ng)
        Fscale(idVsms,ng)=cff*Fscale(idVsms,ng)

and

Code: Select all

        cff=1.0_r8/(rho0*Cp)
        Fscale(idTsur(itemp),ng)=cff*Fscale(idTsur(itemp),ng)
        Fscale(idTbot(itemp),ng)=cff*Fscale(idTbot(itemp),ng)
        Fscale(idSrad,ng)=cff*Fscale(idSrad,ng)
        Fscale(idLdwn,ng)=cff*Fscale(idLdwn,ng)
        Fscale(idLrad,ng)=cff*Fscale(idLrad,ng)
        Fscale(idLhea,ng)=cff*Fscale(idLhea,ng)
        Fscale(idShea,ng)=cff*Fscale(idShea,ng)
        Fscale(iddQdT,ng)=cff*Fscale(iddQdT,ng)

Should your solar radiation be zero? How does your run look?

[wrh]

Kate, thank you so much for your answers. now I got it. Actually, I should read the code carefully.