<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>