[MITgcm-support] EXF

小平翼 kodaira at edu.k.u-tokyo.ac.jp
Mon Dec 11 03:49:01 EST 2017


Hello Martin,

I copied and pasted my data.exf below. The atmospheric temperature field
(atemp) of my input file for EXF package  was intentionally set to 300 K
for a test.  When I set the realistic temperature fields, the max/min
values are similar (i.e., about 1200/600). The wind velocity is also
strange, showing more than 100.

Best,

Tsubasa


# *********************

# External Forcing Data

# *********************

 &EXF_NML_01

# useExfCheckRange  = .TRUE.,

 useExfCheckRange  = .FALSE.,

# readStressOnCgrid = .TRUE.,

 repeatPeriod      = 0,

 exf_iprec         = 64,

 &


# *********************

 &EXF_NML_02

#

# hfluxstartdate1   = 20170901,

# hfluxstartdate2   = 000000,

# hfluxperiod       = 3600.0,

#

# sfluxstartdate1   = 20170901,

# sfluxstartdate2   = 000000,

# sfluxperiod       = 3600.0,

#

# ustressstartdate1 = 20170901,

# ustressstartdate2 = 000000,

# ustressperiod     = 3600.0,

#

# vstressstartdate1 = 20170901,

# vstressstartdate2 = 000000,

# vstressperiod     = 3600.0,

#

 atempstartdate1   = 20170901,

 atempstartdate2   = 000000,

 atempperiod       = 3600.0,


# apressurestartdate1 =

# apressurestartdate2 =

# apressureperiod     =


 aqhstartdate1     = 20170901,

 aqhstartdate2     = 000000,

 aqhperiod         = 3600.0,

#

#evapstartdate1    = 20170901,

#evapstartdate2    = 000000,

#evapperiod        = 3600.0,

#

 precipstartdate1  = 20170901,

 precipstartdate2  = 000000,

 precipperiod      = 3600.0,

#

 uwindstartdate1   = 20170901,

 uwindstartdate2   = 000000,

 uwindperiod       = 3600.0,

#

 vwindstartdate1   = 20170901,

 vwindstartdate2   = 000000,

 vwindperiod       = 3600.0,

#

# swfluxstartdate1  = 20170901,

# swfluxstartdate2  = 000000,

# swfluxperiod      = 3600.0,

#

# lwfluxstartdate1  = 20170901,

# lwfluxstartdate2  = 000000,

# lwfluxperiod      = 3600.0,

#

 swdownstartdate1  = 20170901,

 swdownstartdate2  = 000000,

 swdownperiod      = 86400.0,

#

 lwdownstartdate1  = 20170901,

 lwdownstartdate2  = 000000,

 lwdownperiod      = 86400.0,

#

# climsststartdate1  = 20170901,

# climsststartdate2  = 000000,

# climsstperiod      = 3600.0,

# climsstTauRelax    = 0.0,

#

# climsssstartdate1  = 20170901,

# climsssstartdate2  = 000000,

# climsssperiod      = 3600.0,

# climsssTauRelax    = 0.0,

#

# hfluxfile         = 'hsbc.bin',

 hfluxfile         = ' ',

 sfluxfile         = ' ',

 ustressfile       = ' ',

 vstressfile       = ' ',

 atempfile         = './20170901/temp.surf',

 aqhfile           = './20170901/rh.surf',

 uwindfile         = './20170901/u.surf',

 vwindfile         = './20170901/v.surf',

 evapfile          = ' ',

# apressurefile     = './20170901/psea.surf'

 precipfile        = './20170901/r1h.surf',

 lwfluxfile        = ' ',

 swfluxfile        = ' ',

 lwdownfile        = './20170901/dlwrf.surf',

 swdownfile        = './20170901/dswrf.surf',

 runoffFile        = ' '

 climsstfile       = ' ',

 climsssfile       = ' ',

#

 &


# *********************

 &EXF_NML_03

 &


# *********************

# old open64 compiler (4.2.1) cannot skip this namelist to read in the next
one;

# comment out this namelist (not read).

 &EXF_NML_04

 atemp_lon0      = 124.0,

 atemp_lon_inc   = 0.0625,

 atemp_lat0      = 22.45,

 atemp_lat_inc   = 400*0.05,

 atemp_nlon      = 401,

 atemp_nlat      = 401,

 aqh_lon0      = 124.0,

 aqh_lon_inc   = 0.0625,

 aqh_lat0      = 22.45,

 aqh_lat_inc   = 400*0.05,

 aqh_nlon      = 401,

 aqh_nlat      = 401,

 precip_lon0      = 124.0,

 precip_lon_inc   = 0.0625,

 precip_lat0      = 22.45,

 precip_lat_inc   = 400*0.05,

 precip_nlon      = 401,

 precip_nlat      = 401,

 uwind_lon0      = 124.0,

 uwind_lon_inc   = 0.0625,

 uwind_lat0      = 22.45,

 uwind_lat_inc   = 400*0.05,

 uwind_nlon      = 401,

 uwind_nlat      = 401,

 vwind_lon0      = 124.0,

 vwind_lon_inc   = 0.0625,

 vwind_lat0      = 22.45,

 vwind_lat_inc   = 400*0.05,

 vwind_nlon      = 401,

 vwind_nlat      = 401,

 swdown_lon0      = 114.375,

 swdown_lon_inc   = 1.875,

 swdown_lat0      = 12.3808,

 swdown_lat_inc   = 1.9047,

 swdown_nlon      = 16,

 swdown_nlat      = 14,

 lwdown_lon0      = 114.375,

 lwdown_lon_inc   = 1.875,

 lwdown_lat0      = 12.3808,

 lwdown_lat_inc   = 1.9047,

 lwdown_nlon      = 16,

 lwdown_nlat      = 14,

 &


# *********************

 &EXF_NML_OBCS

 obcsSstartdate1   = 20170901,

 obcsSstartdate2   = 000000,

 obcsSperiod       = 86400.0,

#

 obcsNstartdate1   = 20170901,

 obcsNstartdate2   = 000000,

 obcsNperiod       = 86400.0,

#

 obcsWstartdate1   = 20170901,

 obcsWstartdate2   = 000000,

 obcsWperiod       = 86400.0,

#

 obcsEstartdate1   = 20170901,

 obcsEstartdate2   = 000000,

 obcsEperiod       = 86400.0,

#

# siobSstartdate1   = 20170901,

# siobSstartdate2   = 000000,

# siobSperiod       = 86400.0,

#

# siobNstartdate1   = 20170901,

# siobNstartdate2   = 000000,

# siobNperiod       = 86400.0,

#

# siobWstartdate1   = 20170901,

# siobWstartdate2   = 000000,

# siobWperiod       = 86400.0,

 &



2017-12-11 17:42 GMT+09:00 Martin Losch <Martin.Losch at awi.de>:

> Tsubasa,
>
> can you please provide your complete data.exf. Why would you have a
> min/max of 600/1200 K, what are the max/min of your input fields?
> this is probably not related to the interpolation but rather to either
> your input files rr the settings in data.exf.
>
> Martin
>
> > On 11. Dec 2017, at 04:32, 小平翼 <kodaira at edu.k.u-tokyo.ac.jp> wrote:
> >
> > Hello all,
> >
> > I was trying to use EXF package with the on-the-fly function.
> Unfortuntately, it didn't work as I expected. The values are set much
> larger.  As a quick check, I set the atmospheric temperature fields as
> constant, 300. The monitor function however shows that the maximum of atemp
> is 4times larger and the mean is also 2 times larger, as follows.
> >
> > (PID.TID 1000.0001) %MON exf_atemp_max                =
>  1.2000000000000E+03
> > (PID.TID 1000.0001) %MON exf_atemp_min                =
>  6.0000000000027E+02
> > (PID.TID 1000.0001) %MON exf_atemp_mean               =
>  9.1028890280222E+02
> > (PID.TID 1000.0001) %MON exf_atemp_sd                 =
>  1.7286119388930E+02
> > (PID.TID 1000.0001) %MON exf_atemp_del2               =
>  3.4626736833088E-01
> >
> > I'm wondering I may be missing something.
> > Any advice or help would be appreciated.
> > Related part of the data.exf is following
> >
> > &EXF_NML_04
> >  atemp_lon0      = 124.0,
> >  atemp_lon_inc   = 0.25,
> >  atemp_lat0      = 22.45,
> >  atemp_lat_inc   = 100*0.2,
> >  atemp_nlon      = 101,
> >  atemp_nlat      = 101,
> >
> > Thank you,
> >
> > Tsubasa
> > --
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>



-- 
///////////////////////////////////////////////
////////////////////////////////////////////////////////
小平翼  助教、海洋技術環境学専攻、東大院新領域
居室 〒277-8563 千葉県柏市柏の葉5-1-5  環境棟 214
Tsubasa Kodaira  (Ph.D), Assistant Professor
Department of Ocean Technology, Policy and Environment
Graduate School of Frontier Sciences, The University of Tokyo
TEL 04-7136-4672, E-mail kodaira at edu.k.u-tokyo.ac.jp
///////////////////////////////////////////////
////////////////////////////////////////////////////////
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20171211/25180b8b/attachment-0001.html>


More information about the MITgcm-support mailing list