[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