[MITgcm-support] error with exf "use u, v_wind components but not wind-stress statement executed "
李志远
oceanlizy at 163.com
Wed Jan 8 13:20:42 EST 2014
From: lizhiyuan <oceanlizy at 163.com>
To: MITgcm <mitgcm-support at mitgcm.org>
Subject:Re: [MITgcm-support] error with exf "use u, v_wind components but not wind-stress statement executed "thank you for your suggestion , Jean-Michel ! The run-time error is gone now ! but,there seems some trouble in the output STDOUT file:PID.TID 0009.0001) // =======================================================
(PID.TID 0009.0001) // End MONITOR dynamic field statistics
(PID.TID 0009.0001) // =======================================================
(PID.TID 0009.0001) MDS_READ_FIELD: opening global file: u2.bin
(PID.TID 0009.0001) MDS_READ_FIELD: opening global file: u2.bin
(PID.TID 0009.0001) // =======================================================
(PID.TID 0009.0001) // Begin MONITOR EXF statistics
(PID.TID 0009.0001) // =======================================================
(PID.TID 0009.0001) %MON exf_tsnumber = 0
(PID.TID 0009.0001) %MON exf_time_sec = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_ustress_max = 4.5739727770205E-02
(PID.TID 0009.0001) %MON exf_ustress_min = 4.5739727770205E-02
(PID.TID 0009.0001) %MON exf_ustress_mean = 4.5739727770205E-02
(PID.TID 0009.0001) %MON exf_ustress_sd = 4.1633363423443E-17
(PID.TID 0009.0001) %MON exf_ustress_del2 = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vstress_max = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vstress_min = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vstress_mean = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vstress_sd = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vstress_del2 = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_hflux_max = 8.9001119177706E+02
(PID.TID 0009.0001) %MON exf_hflux_min = 8.9001119177706E+02
(PID.TID 0009.0001) %MON exf_hflux_mean = 8.9001119177706E+02
(PID.TID 0009.0001) %MON exf_hflux_sd = 2.5011104298756E-12
(PID.TID 0009.0001) %MON exf_hflux_del2 = 2.0007645728653E+00
(PID.TID 0009.0001) %MON exf_sflux_max = 2.4150356729736E-07
(PID.TID 0009.0001) %MON exf_sflux_min = 2.4150356729736E-07
(PID.TID 0009.0001) %MON exf_sflux_mean = 2.4150356729736E-07
(PID.TID 0009.0001) %MON exf_sflux_sd = 3.7057691442376E-22
(PID.TID 0009.0001) %MON exf_sflux_del2 = 5.4290528718452E-10
(PID.TID 0009.0001) %MON exf_uwind_max = 5.0000000000000E+00
(PID.TID 0009.0001) %MON exf_uwind_min = 5.0000000000000E+00
(PID.TID 0009.0001) %MON exf_uwind_mean = 5.0000000000000E+00
(PID.TID 0009.0001) %MON exf_uwind_sd = 3.6415315207705E-14
(PID.TID 0009.0001) %MON exf_uwind_del2 = 1.1240109064643E-02
(PID.TID 0009.0001) %MON exf_vwind_max = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vwind_min = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vwind_mean = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vwind_sd = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_vwind_del2 = 0.0000000000000E+00
(PID.TID 0009.0001) %MON exf_wspeed_max = 5.0000000000000E+00
(PID.TID 0009.0001) %MON exf_wspeed_min = 5.0000000000000E+00
(PID.TID 0009.0001) %MON exf_wspeed_mean = 5.0000000000000E+00
(PID.TID 0009.0001) %MON exf_wspeed_sd = 3.6415315207705E-14
(PID.TID 0009.0001) %MON exf_wspeed_del2 = 1.1240109064643E-02
(PID.TID 0009.0001) %MON exf_evap_max = 2.4150356729736E-07
(PID.TID 0009.0001) %MON exf_evap_min = 2.4150356729736E-07
(PID.TID 0009.0001) %MON exf_evap_mean = 2.4150356729736E-07
(PID.TID 0009.0001) %MON exf_evap_sd = 3.7057691442376E-22
(PID.TID 0009.0001) %MON exf_evap_del2 = 5.4290528718452E-10
(PID.TID 0009.0001) // =======================================================
(PID.TID 0009.0001) // End MONITOR EXF statistics
(PID.TID 0009.0001) // =======================================================
my data.exf file is :
#
# *********************
# External Forcing Data
# *********************
&EXF_NML_01
#
useExfCheckRange = .TRUE.,
# repeatPeriod = 31104000.0,
exf_iprec = 64,
#
&
# *********************
&EXF_NML_02
#
uwindstartdate1 = 20061215,
uwindstartdate2 = 0000000,
uwindperiod = 2592000.0,
#
vwindstartdate1 = 20061215,
vwindstartdate2 = 00000000,
vwindperiod = 2592000.0,
#
uwindfile = 'u2.bin',
# vwindfile = 'v14.bin',
&
in my data.exf , I only input a uwind file ='u2.bin' , without 'hflux' and any other external forces .But in the STDOUT file: exf_hflux_max = 8.9001119177706E+02 , it is obviously out of the range of (-250,600) .why exf_hflux is not zero ? I think the value of it is concerned with initial sea surface temperature ,namely initial temperature in the data file :
&PARM01
tRef= 26.5, 25.5, 24.4, 23.1, 21.9,
20.5, 18.7, 16.4, 13.6, 11.3,
8.9, 6.5, 4.6, 4.6, 4.6,
4.6, 4.6, 4.6, 4.6, 4.6,
4.6, 4.6, 4.6, 4.6, 4.6,
4.6, 4.6, 4.6, 4.6, 4.6,
4.6, 4.6, 4.6, 4.6, 4.6,
4.6,
as the model running for some steps ,the exf_hflux will reduce and tend to be a reasonable value ?
thanks again !
lizhiyuan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20140109/2dddfc9d/attachment.htm>
More information about the MITgcm-support
mailing list