<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=gb2312">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi Matt,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I have checked the diagnosed wind stress, the wind stress pattern and value are OK.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Strangely, I can successfully run the 4-degree configuration of global_with_exf, and the currents seem OK, but when I increased the resolution to 1grid*1grid, there are little currents outside the equator (even at the time of 100 years).</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Here is my data.exf and <span style="font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255); display: inline !important">EXF_OPTIONS.h.</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255); display: inline !important"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255); display: inline !important"><span># *********************<br>
</span>
<div># External Forcing Data<br>
</div>
<div># *********************<br>
</div>
<div> &EXF_NML_01<br>
</div>
<div> exf_iprec         = 32,<br>
</div>
<div> exf_debugLev      = 3,<br>
</div>
<div> &<br>
</div>
<div><br>
</div>
<div># *********************<br>
</div>
<div> &EXF_NML_02<br>
</div>
<div>#hfluxstartdate1   = 19920101,<br>
</div>
<div>#hfluxstartdate2   = 000000,<br>
</div>
<div> hfluxperiod       = -12.,<br>
</div>
<div>#<br>
</div>
<div>#sfluxstartdate1   = 19920101,<br>
</div>
<div>#sfluxstartdate2   = 000000,<br>
</div>
<div> sfluxperiod       = -12.,<br>
</div>
<div>#<br>
</div>
<div>#ustressstartdate1 = 19920101,<br>
</div>
<div>#ustressstartdate2 = 000000,<br>
</div>
<div> ustressperiod     = -12.,<br>
</div>
<div>#<br>
</div>
<div>#vstressstartdate1 = 19920101,<br>
</div>
<div>#vstressstartdate2 = 000000,<br>
</div>
<div> vstressperiod     = -12.,<br>
</div>
<div>#<br>
</div>
<div>#climsststartdate1 = 19920101,<br>
</div>
<div>#climsststartdate2 = 000000,<br>
</div>
<div> climsstperiod     = -12.,<br>
</div>
<div># 2 months restoring timescale for temperature<br>
</div>
<div> climsstTauRelax   = 5184000.,<br>
</div>
<div>#<br>
</div>
<div>#climsssstartdate1 = 19920101,<br>
</div>
<div>#climsssstartdate2 = 000000,<br>
</div>
<div> climsssperiod     = -12.,<br>
</div>
<div># 6 months restoring timescale for salinity<br>
</div>
<div> climsssTauRelax   = 15552000.,<br>
</div>
<div>#<br>
</div>
<div> hfluxfile    = 'qnet.bin',<br>
</div>
<div> sfluxfile    = 'emp.bin',<br>
</div>
<div># here we deviate from the general pattern and prescribe a forcing fiels<br>
</div>
<div># with (i,j)=(1,1) at the NW (.flipped) instead of the SW (.bin) corner<br>
</div>
<div># to test the correpsonding code, see also below<br>
</div>
<div> ustressfile  = 'taux.bin',<br>
</div>
<div> vstressfile  = 'tauy.bin',<br>
</div>
<div> climsstfile  = 'SST.bin',<br>
</div>
<div> climsssfile  = 'SSS.bin',<br>
</div>
<div>#<br>
</div>
<div> &<br>
</div>
<div><br>
</div>
<div># *********************<br>
</div>
<div> &EXF_NML_03<br>
</div>
<div> &<br>
</div>
<div><br>
</div>
<div># *********************<br>
</div>
<div> &EXF_NML_04<br>
</div>
<div> hflux_lon0    = 0.,<br>
</div>
<div> hflux_lon_inc = 1.,<br>
</div>
<div> hflux_lat0    = -80.,<br>
</div>
<div> hflux_lat_inc = 159*1.,<br>
</div>
<div> hflux_nlon    = 360,<br>
</div>
<div> hflux_nlat    = 160,<br>
</div>
<div>#<br>
</div>
<div> sflux_lon0    = 0.,<br>
</div>
<div> sflux_lon_inc = 1.,<br>
</div>
<div> sflux_lat0    = -80.,<br>
</div>
<div> sflux_lat_inc = 159*1.,<br>
</div>
<div> sflux_nlon    = 360,<br>
</div>
<div> sflux_nlat    = 160,<br>
</div>
<div>#<br>
</div>
<div> ustress_lon0    = 0.,<br>
</div>
<div> ustress_lon_inc = 1.,<br>
</div>
<div> ustress_lat0    = -80.,<br>
</div>
<div> ustress_lat_inc = 159*1.,<br>
</div>
<div> ustress_nlon    = 360,<br>
</div>
<div> ustress_nlat    = 160,<br>
</div>
<div>#<br>
</div>
<div> vstress_lon0    = 0.,<br>
</div>
<div> vstress_lon_inc = 1.,<br>
</div>
<div> vstress_lat0    = -80.,<br>
</div>
<div> vstress_lat_inc = 159*1.,<br>
</div>
<div> vstress_nlon    = 360,<br>
</div>
<div> vstress_nlat    = 160,<br>
</div>
<div>#<br>
</div>
<div> climsst_lon0    = 0.,<br>
</div>
<div> climsst_lon_inc = 1.,<br>
</div>
<div> climsst_lat0    = -80.,<br>
</div>
<div> climsst_lat_inc = 159*1.,<br>
</div>
<div> climsst_nlon    = 360,<br>
</div>
<div> climsst_nlat    = 160,<br>
</div>
<div>#<br>
</div>
<div> climsss_lon0    = 0.,<br>
</div>
<div> climsss_lon_inc = 1.,<br>
</div>
<div> climsss_lat0    = -80.,<br>
</div>
<div> climsss_lat_inc = 159*1.,<br>
</div>
<div> climsss_nlon    = 360,<br>
</div>
<div> climsss_nlat    = 160,<br>
</div>
<div>#<br>
</div>
<div> &<br>
</div>
<span></span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255); display: inline !important"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span style="font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255); display: inline !important"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>C $Header: /u/gcmpack/MITgcm/verification/global_with_exf/code/EXF_OPTIONS.h,v 1.8 2013/10/05 19:53:17 jmc Exp $<br>
</span>
<div>C $Name: checkpoint65o $<br>
</div>
<div><br>
</div>
<div>CBOP<br>
</div>
<div>C !ROUTINE: EXF_OPTIONS.h<br>
</div>
<div>C !INTERFACE:<br>
</div>
<div>C #include "EXF_OPTIONS.h"<br>
</div>
<div><br>
</div>
<div>C !DESCRIPTION:<br>
</div>
<div>C *==================================================================*<br>
</div>
<div>C | CPP options file for EXternal Forcing (EXF) package:<br>
</div>
<div>C | Control which optional features to compile in this package code.<br>
</div>
<div>C *==================================================================*<br>
</div>
<div>CEOP<br>
</div>
<div><br>
</div>
<div>#ifndef EXF_OPTIONS_H<br>
</div>
<div>#define EXF_OPTIONS_H<br>
</div>
<div>#include "PACKAGES_CONFIG.h"<br>
</div>
<div>#include "CPP_OPTIONS.h"<br>
</div>
<div><br>
</div>
<div>#ifdef ALLOW_EXF<br>
</div>
<div>#ifdef ECCO_CPPOPTIONS_H<br>
</div>
<div><br>
</div>
<div>C-- When multi-package option-file ECCO_CPPOPTIONS.h is used (directly included<br>
</div>
<div>C    in CPP_OPTIONS.h), this option file is left empty since all options that<br>
</div>
<div>C   are specific to this package are assumed to be set in ECCO_CPPOPTIONS.h<br>
</div>
<div><br>
</div>
<div>#else /* ndef ECCO_CPPOPTIONS_H */<br>
</div>
<div><br>
</div>
<div>C-- Package-specific Options & Macros go here<br>
</div>
<div><br>
</div>
<div>C   pkg/exf CPP options:<br>
</div>
<div>C   --------------------<br>
</div>
<div>C<br>
</div>
<div>C   > ( EXF_VERBOSE ) < replaced with run-time, logical parameter "exf_verbose".<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_ATM_WIND <<<<br>
</div>
<div>C       If defined, 10-m wind fields can be read-in from files.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_ATM_TEMP <<<<br>
</div>
<div>C       If defined, atmospheric temperature and specific<br>
</div>
<div>C       humidity fields can be read-in from files.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_DOWNWARD_RADIATION <<<<br>
</div>
<div>C       If defined, downward long-wave and short-wave radiation<br>
</div>
<div>C       can be read-in form files or computed from lwflux and swflux.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_ZENITHANGLE <<<<br>
</div>
<div>C       If defined, ocean albedo varies with the zenith angle, and<br>
</div>
<div>C       incoming fluxes at the top of the atmosphere are computed<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_BULKFORMULAE <<<<br>
</div>
<div>C       Allows the use of bulk formulae in order to estimate<br>
</div>
<div>C       turbulent and radiative fluxes at the ocean surface.<br>
</div>
<div>C<br>
</div>
<div>C   >>> EXF_READ_EVAP <<<<br>
</div>
<div>C       If defined, evaporation fields are read-in, rather than<br>
</div>
<div>C       computed from atmospheric state.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_RUNOFF <<<<br>
</div>
<div>C       If defined, river and glacier runoff can be read-in from files.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_RUNOFTEMP <<<<br>
</div>
<div>C       If defined, river and glacier runoff temperature<br>
</div>
<div>C       can be read-in from files.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ATMOSPHERIC_LOADING <<<<br>
</div>
<div>C       If defined, atmospheric pressure can be read-in from files.<br>
</div>
<div>C   WARNING: this flag is set (define/undef) in CPP_OPTIONS.h<br>
</div>
<div>C            and cannot be changed here (in EXF_OPTIONS)<br>
</div>
<div>C<br>
</div>
<div>C   >>> EXF_SEAICE_FRACTION <<<<br>
</div>
<div>C       If defined, seaice fraction can be read-in from files (areaMaskFile)<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_CLIMSST_RELAXATION <<<<br>
</div>
<div>C       Allow the relaxation to a monthly climatology of sea surface<br>
</div>
<div>C       temperature, e.g. the Reynolds climatology.<br>
</div>
<div>C<br>
</div>
<div>C   >>> ALLOW_CLIMSSS_RELAXATION <<<<br>
</div>
<div>C       Allow the relaxation to a monthly climatology of sea surface<br>
</div>
<div>C       salinity, e.g. the Levitus climatology.<br>
</div>
<div>C<br>
</div>
<div>C   >>> USE_EXF_INTERPOLATION <<<<br>
</div>
<div>C       Allows specification of arbitrary Cartesian input grids.<br>
</div>
<div>C<br>
</div>
<div>C   ====================================================================<br>
</div>
<div>C<br>
</div>
<div>C       The following CPP options:<br>
</div>
<div>C<br>
</div>
<div>C          ALLOW_ATM_WIND              (WIND)<br>
</div>
<div>C          ALLOW_ATM_TEMP              (TEMP)<br>
</div>
<div>C          ALLOW_DOWNWARD_RADIATION    (DOWN)<br>
</div>
<div>C          ALLOW_BULKFORMULAE          (BULK)<br>
</div>
<div>C          EXF_READ_EVAP               (EVAP)<br>
</div>
<div>C<br>
</div>
<div>C       permit the ocean-model forcing configurations listed in the<br>
</div>
<div>C       table below.  The first configuration is the default,<br>
</div>
<div>C       flux-forced, ocean model.  The next four are stand-alone<br>
</div>
<div>C       configurations that use pkg/exf, open-water bulk formulae to<br>
</div>
<div>C       compute the missing surface fluxes from atmospheric variables.<br>
</div>
<div>C       The last four configurations can be used in conjunction with<br>
</div>
<div>C       pkg/seaice to model ice-covered regions.  The forcing fields<br>
</div>
<div>C       in the rightmost column are defined in exf_fields.<br>
</div>
<div>C<br>
</div>
<div>C<br>
</div>
<div>C    WIND |TEMP |DOWN |BULK |EVAP |            actions<br>
</div>
<div>C    -----|-----|-----|-----|-----|-------------------------------------<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C      -  |  -  |  -  |  -  |  -  | Read-in ustress, vstress, hflux,<br>
</div>
<div>C         |     |     |     |     | swflux, and sflux.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def | def | def | def |  -  | Read-in uwind, vwind, atemp, aqh,<br>
</div>
<div>C         |     |     |     |     | swdown, lwdown, precip, and runoff.<br>
</div>
<div>C         |     |     |     |     | Compute ustress, vstress, hflux,<br>
</div>
<div>C         |     |     |     |     | swflux, and sflux.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def | def |  -  | def |  -  | Read-in uwind, vwind, atemp, aqh,<br>
</div>
<div>C         |     |     |     |     | swflux, lwflux, precip, and runoff.<br>
</div>
<div>C         |     |     |     |     | Compute ustress, vstress, hflux,<br>
</div>
<div>C         |     |     |     |     | and sflux.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def |  -  |  -  | def |  -  | Read-in uwind, vwind, hflux,<br>
</div>
<div>C         |     |     |     |     | swflux, and sflux.<br>
</div>
<div>C         |     |     |     |     | Compute ustress and vstress.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C      -  | def |  -  | def |  -  | Read-in ustress, vstress, atemp,<br>
</div>
<div>C         |     |     |     |     | aqh, swflux, lwflux, precip, and<br>
</div>
<div>C         |     |     |     |     | runoff.  Compute hflux and sflux.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def | def |  -  |  -  | def | Read-in uwind, vwind, atemp, aqh,<br>
</div>
<div>C         |     |     |     |     | swflux, lwflux, precip, runoff,<br>
</div>
<div>C         |     |     |     |     | and evap.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def | def |  -  | def |  -  | Read-in uwind, vwind, atemp, aqh,<br>
</div>
<div>C         |     |     |     |     | swflux, lwflux, precip, and runoff.<br>
</div>
<div>C         |     |     |     |     | Compute open-water ustress, vstress,<br>
</div>
<div>C         |     |     |     |     | hflux, swflux, and evap.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def | def | def |  -  | def | Read-in uwind, vwind, atemp, aqh,<br>
</div>
<div>C         |     |     |     |     | swdown, lwdown, precip, runoff,<br>
</div>
<div>C         |     |     |     |     | and evap.<br>
</div>
<div>C         |     |     |     |     |<br>
</div>
<div>C     def | def | def | def |  -  | Read-in uwind, vwind, atemp, aqh,<br>
</div>
<div>C         |     |     |     |     | swdown, lwdown, precip, and runoff.<br>
</div>
<div>C         |     |     |     |     | Compute open-water ustress, vstress,<br>
</div>
<div>C         |     |     |     |     | hflux, swflux, and evap.<br>
</div>
<div>C<br>
</div>
<div>C   ====================================================================<br>
</div>
<div><br>
</div>
<div>C   Bulk formulae related flags.<br>
</div>
<div>#undef  ALLOW_ATM_TEMP<br>
</div>
<div>#undef  ALLOW_ATM_WIND<br>
</div>
<div>#undef  ALLOW_DOWNWARD_RADIATION<br>
</div>
<div>#undef  ALLOW_RUNOFF<br>
</div>
<div>#undef  ALLOW_RUNOFTEMP<br>
</div>
<div>#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))<br>
</div>
<div># define ALLOW_BULKFORMULAE<br>
</div>
<div># undef  ALLOW_BULK_LARGEYEAGER04<br>
</div>
<div>#endif<br>
</div>
<div><br>
</div>
<div>C   Zenith Angle/Albedo related flags.<br>
</div>
<div>#ifdef ALLOW_DOWNWARD_RADIATION<br>
</div>
<div># undef ALLOW_ZENITHANGLE<br>
</div>
<div>#endif<br>
</div>
<div><br>
</div>
<div>C   Use ocean_emissivity*lwdown in lwFlux. This flag should be defined<br>
</div>
<div>C   unless to reproduce old results (obtained with inconsistent old code)<br>
</div>
<div>#ifdef ALLOW_DOWNWARD_RADIATION<br>
</div>
<div># define EXF_LWDOWN_WITH_EMISSIVITY<br>
</div>
<div>#endif<br>
</div>
<div><br>
</div>
<div>C   Relaxation to monthly climatologies.<br>
</div>
<div>#define ALLOW_CLIMSST_RELAXATION<br>
</div>
<div>#define ALLOW_CLIMSSS_RELAXATION<br>
</div>
<div><br>
</div>
<div>C   Allows to read-in seaice fraction from files (areaMaskFile)<br>
</div>
<div>#undef EXF_SEAICE_FRACTION<br>
</div>
<div><br>
</div>
<div>C   Use spatial interpolation to interpolate<br>
</div>
<div>C   forcing files from input grid to model grid.<br>
</div>
<div>#define USE_EXF_INTERPOLATION<br>
</div>
<div>C   for interpolated vector fields, rotate towards model-grid axis<br>
</div>
<div>C   using old rotation formulae (instead of grid-angles)<br>
</div>
<div>#undef EXF_USE_OLD_VEC_ROTATION<br>
</div>
<div>C   for interpolation around N & S pole, use the old formulation<br>
</div>
<div>C   (no pole symmetry, single vector-comp interp, reset to 0 zonal-comp @ N.pole)<br>
</div>
<div>#undef EXF_USE_OLD_INTERP_POLE<br>
</div>
<div><br>
</div>
<div>#undef EXF_INTERP_USE_DYNALLOC<br>
</div>
<div>#if ( defined (EXF_INTERP_USE_DYNALLOC) && defined (USING_THREADS) )<br>
</div>
<div># define EXF_IREAD_USE_GLOBAL_POINTER<br>
</div>
<div>#endif<br>
</div>
<div><br>
</div>
<div>#endif /* ndef ECCO_CPPOPTIONS_H */<br>
</div>
<div>#endif /* ALLOW_EXF */<br>
</div>
<div>#endif /* EXF_OPTIONS_H */<br>
</div>
<span></span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thank you so much.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Fei</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>发件人:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> 代表 mitgcm-support-request@mitgcm.org <mitgcm-support-request@mitgcm.org><br>
<b>发送时间:</b> 2020年1月8日 1:00<br>
<b>收件人:</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>主题:</b> MITgcm-support Digest, Vol 199, Issue 5</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Send MITgcm-support mailing list submissions to<br>
        mitgcm-support@mitgcm.org<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
or, via email, send a message with subject or body 'help' to<br>
        mitgcm-support-request@mitgcm.org<br>
<br>
You can reach the person managing the list at<br>
        mitgcm-support-owner@mitgcm.org<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of MITgcm-support digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Re: a question about the currents (Matthew Mazloff)<br>
   2. A question about the currents (pan fei)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 6 Jan 2020 19:26:50 -0800<br>
From: Matthew Mazloff <mmazloff@ucsd.edu><br>
To: mitgcm-support@mitgcm.org<br>
Subject: Re: [MITgcm-support] a question about the currents<br>
Message-ID: <4E25C96D-EEFE-412B-A45D-B5C4F914E625@ucsd.edu><br>
Content-Type: text/plain;       charset=us-ascii<br>
<br>
Hi Fei<br>
<br>
How does your wind stress look? That is the first thing I would check.<br>
<br>
Matt<br>
<br>
<br>
> On Jan 6, 2020, at 7:21 PM, 979583702 <979583702@qq.com> wrote:<br>
> <br>
> Hi everyone,<br>
> <br>
> Recently, I ran a global ocean model configuration which is similar with the global_with_exf but with 1 grid*1 grid resolution. Strangely, I found the currents  mainly located on the equator (e.g., strong westward currents over eastern equatorial Pacific),
 and there are no western boundary currents and  ACC. <br>
> <br>
> Attached are the currents at the time of 40 years.<br>
> <br>
> Here is my data file:<br>
> <br>
> # ====================<br>
> # | Model parameters |<br>
> # ====================<br>
> #<br>
> # Continuous equation parameters<br>
>  &PARM01<br>
>  tRef= 13.72 , 12.82 , 11.57 , 10.35 , 9.38 , <br>
>        8.62 , 7.99 , 7.45 , 6.91 , 6.45 , <br>
>        5.99 , 5.63 , 5.26 , 4.97 , 4.68 , <br>
>        4.44 , 4.20 , 4.06 , 3.85 , 3.68 , <br>
>        3.51 , 3.22 , 2.99 , 2.80 , 2.62 , <br>
>        2.47 , 1.93 , 1.61 , 1.27 , 1.23 ,<br>
>  sRef = 30*34.,     <br>
>  viscAr=1.E-4,<br>
> # viscAh=2.E2,<br>
>  viscAh=1.E0,<br>
>  viscAhGrid=2.E-2,<br>
>  diffKhT=0.,<br>
>  diffKrT=2.E-5,<br>
>  diffKhS=0.,<br>
>  diffKrS=2.E-5,<br>
>  rhonil=1035.,<br>
>  rhoConstFresh=1000.,<br>
>  eosType = 'JMD95Z',<br>
>  ivdc_kappa=100.,<br>
>  implicitDiffusion=.TRUE.,<br>
>  allowFreezing=.TRUE.,<br>
>  exactConserv=.TRUE.,<br>
>  useRealFreshWaterFlux=.TRUE.,<br>
>  useCDscheme=.TRUE.,<br>
> # turn on looped cells<br>
>  hFacMin=.05,<br>
>  hFacMindr=50.,<br>
>  useSingleCpuIO=.TRUE.,<br>
> # set precision of data files<br>
>  readBinaryPrec=32,<br>
>  &<br>
> <br>
> # Elliptic solver parameters<br>
>  &PARM02<br>
>  cg2dMaxIters=500,<br>
>  cg2dTargetResidual=1.E-13,<br>
>  &<br>
> <br>
> # Time stepping parameters<br>
>  &PARM03<br>
>  nIter0 =      0,<br>
>  nTimeSteps = 864000,<br>
> # 100 years of integration will yield a reasonable flow field<br>
> # startTime  =          0.,<br>
> # endTime    = 3110400000.,<br>
>  deltaTmom = 300.0,<br>
>  tauCD =     321428.,<br>
>  deltaTtracer= 1800.0,<br>
>  deltaTClock = 3600.0,<br>
> # if you are using a version later than checkpoint45d on the main branch<br>
> # you can uncomment the following line and increase the time step<br>
> # deltaTtracer and deltaTClock to 172800.0 as well to speed up the<br>
> # asynchronous time stepping<br>
> # deltaTfreesurf = 172800.0,<br>
>  abEps = 0.1,<br>
>  pChkptFreq= 31104000.,<br>
> # dumpFreq=   311040000.,<br>
> # dumpFreq=   864000.,<br>
> # taveFreq=   311040000.,<br>
>  taveFreq=   2592000.,<br>
>  monitorFreq=720.,<br>
>  &<br>
> <br>
> # Gridding parameters<br>
>  &PARM04<br>
>  usingSphericalPolarGrid=.TRUE.,<br>
>  delR= 10., 15., 20., 30., 40., <br>
>       50., 50., 60., 60., 60.,<br>
>       60., 60., 60., 60., 60.,<br>
>       60., 60., 60., 70., 70.,<br>
>       70., 100., 140., 200., 250.,<br>
>       320., 420., 550., 650., 800.,<br>
>  ygOrigin=-80.,<br>
>  xgOrigin=0.,<br>
>  dySpacing=1.,<br>
>  dxSpacing=1.,<br>
>  &<br>
> <br>
> # Input datasets<br>
>  &PARM05<br>
>  bathyFile=      'bathymetry.bin',<br>
> # hydrogThetaFile='lev_t.bin',<br>
> # hydrogSaltFile= 'lev_s.bin',<br>
>  &<br>
> Thank you so much.<br>
> Best,<br>
> Fei<br>
> <currents.png>_______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 7 Jan 2020 04:26:52 +0000<br>
From: pan fei <feiamanda@hotmail.com><br>
To: "mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
Subject: [MITgcm-support] A question about the currents<br>
Message-ID:<br>
        <PS2PR06MB34649FCB268ADE0F58A8CDBDA33F0@PS2PR06MB3464.apcprd06.prod.outlook.com><br>
        <br>
Content-Type: text/plain; charset="gb2312"<br>
<br>
I'm really sorry for the disturbance. I found the text of my previous email consisted of error codes, maybe it is related to my email. So I use another email to resend the question. Thank you so much.<br>
<br>
Hi everyone,<br>
<br>
Recently, I ran a global ocean model configuration which is similar with the global_with_exf but with 1 grid*1 grid resolution. Strangely, I found the currents were mainly located on the equator (e.g., strong westward currents over eastern equatorial Pacific),
 and there are no western boundary currents and  ACC.<br>
<br>
Attached are the currents at the time of 40 years.<br>
<br>
Here is my data file:<br>
<br>
# ====================<br>
# | Model parameters |<br>
# ====================<br>
#<br>
# Continuous equation parameters<br>
 &PARM01<br>
 tRef= 13.72 , 12.82 , 11.57 , 10.35 , 9.38 ,<br>
       8.62 , 7.99 , 7.45 , 6.91 , 6.45 ,<br>
       5.99 , 5.63 , 5.26 , 4.97 , 4.68 ,<br>
       4.44 , 4.20 , 4.06 , 3.85 , 3.68 ,<br>
       3.51 , 3.22 , 2.99 , 2.80 , 2.62 ,<br>
       2.47 , 1.93 , 1.61 , 1.27 , 1.23 ,<br>
 sRef = 30*34.,<br>
 viscAr=1.E-4,<br>
# viscAh=2.E2,<br>
 viscAh=1.E0,<br>
 viscAhGrid=2.E-2,<br>
 diffKhT=0.,<br>
 diffKrT=2.E-5,<br>
 diffKhS=0.,<br>
 diffKrS=2.E-5,<br>
 rhonil=1035.,<br>
 rhoConstFresh=1000.,<br>
 eosType = 'JMD95Z',<br>
 ivdc_kappa=100.,<br>
 implicitDiffusion=.TRUE.,<br>
 allowFreezing=.TRUE.,<br>
 exactConserv=.TRUE.,<br>
 useRealFreshWaterFlux=.TRUE.,<br>
 useCDscheme=.TRUE.,<br>
# turn on looped cells<br>
 hFacMin=.05,<br>
 hFacMindr=50.,<br>
 useSingleCpuIO=.TRUE.,<br>
# set precision of data files<br>
 readBinaryPrec=32,<br>
 &<br>
<br>
# Elliptic solver parameters<br>
 &PARM02<br>
 cg2dMaxIters=500,<br>
 cg2dTargetResidual=1.E-13,<br>
 &<br>
<br>
# Time stepping parameters<br>
 &PARM03<br>
 nIter0 =      0,<br>
 nTimeSteps = 864000,<br>
# 100 years of integration will yield a reasonable flow field<br>
# startTime  =          0.,<br>
# endTime    = 3110400000.,<br>
 deltaTmom = 300.0,<br>
 tauCD =     321428.,<br>
 deltaTtracer= 1800.0,<br>
 deltaTClock = 3600.0,<br>
# if you are using a version later than checkpoint45d on the main branch<br>
# you can uncomment the following line and increase the time step<br>
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the<br>
# asynchronous time stepping<br>
# deltaTfreesurf = 172800.0,<br>
 abEps = 0.1,<br>
 pChkptFreq= 31104000.,<br>
# dumpFreq=   311040000.,<br>
# dumpFreq=   864000.,<br>
# taveFreq=   311040000.,<br>
 taveFreq=   2592000.,<br>
 monitorFreq=720.,<br>
 &<br>
<br>
# Gridding parameters<br>
 &PARM04<br>
 usingSphericalPolarGrid=.TRUE.,<br>
 delR= 10., 15., 20., 30., 40.,<br>
      50., 50., 60., 60., 60.,<br>
      60., 60., 60., 60., 60.,<br>
      60., 60., 60., 70., 70.,<br>
      70., 100., 140., 200., 250.,<br>
      320., 420., 550., 650., 800.,<br>
 ygOrigin=-80.,<br>
 xgOrigin=0.,<br>
 dySpacing=1.,<br>
 dxSpacing=1.,<br>
 &<br>
<br>
# Input datasets<br>
 &PARM05<br>
 bathyFile=      'bathymetry.bin',<br>
# hydrogThetaFile='lev_t.bin',<br>
# hydrogSaltFile= 'lev_s.bin',<br>
 &<br>
Thank you so much.<br>
Best,<br>
Fei<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200107/f02b4a90/attachment-0001.html">http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200107/f02b4a90/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
<br>
------------------------------<br>
<br>
End of MITgcm-support Digest, Vol 199, Issue 5<br>
**********************************************<br>
</div>
</span></font></div>
</body>
</html>