<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto">Hi Pan,<div><br></div><div>Your case reminds me of a similar experience. Have you enabled the KPP package which will transfer surface momentum downward efficiently through the mixing layer, and thus modify your surface current to some extent?</div><div><br><div dir="ltr">发自我的iPad</div><div id="AppleMailSignature" dir="ltr">Yu-Kun Qian</div><div dir="ltr"><br><blockquote type="cite">在 2020年1月9日，22:55，pan fei <feiamanda@hotmail.com> 写道：<br><br></blockquote></div><blockquote type="cite"><div dir="ltr">
<meta http-equiv="Content-Type" content="text/html; charset=gb2312">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi Matt and ED,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thanks for your reply. Today I have checked the output again. Interesingly, I found that only the currents of the first layer are strange (large value near the equator, while little value outside the equator). But the currents of the second to 30th layer (there
are 30 layers of the model) look OK (there are strong western boundary currents and ACC in these layers). I'm quite confused about why the currents of the first layer are so strange.</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);">
Fei</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日 10:30<br>
<b>收件人:</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>主题:</b> MITgcm-support Digest, Vol 199, Issue 6</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. Postdoctoral Research Associate in Physical Oceanography<br>
(Maarten Buijsman)<br>
2. ??: MITgcm-support Digest, Vol 199, Issue 5 (pan fei)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 7 Jan 2020 17:07:55 +0000<br>
From: Maarten Buijsman <Maarten.Buijsman@usm.edu><br>
To: "mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
Subject: [MITgcm-support] Postdoctoral Research Associate in Physical<br>
Oceanography<br>
Message-ID:<br>
<BY5PR10MB384445FA3FEF612591895E60803F0@BY5PR10MB3844.namprd10.prod.outlook.com><br>
<br>
Content-Type: text/plain; charset="windows-1252"<br>
<br>
A postdoctoral position in physical oceanography is available in the Division of Marine Science in the School of Ocean Science and Engineering (SOSE) at The University of Southern Mississippi. The successful candidate will work on the three-year Task Force
Ocean (TFO) project ?Modeling, characterizing, and predicting effects of internal gravity waves on acoustic propagation on basin to global scales?, which is supported by a grant from the Office of Naval Research (ONR).<br>
<br>
<br>
<br>
Applications must be submitted through the jobs.usm.edu candidate portal (<a href=""></a>https://usm.csod.com/ats/careersite/JobDetails.aspx?id=1271&site=1; job#: req1271 ). For full consideration, send by 1/31/2020: 1) letter of interest, 2) vitae, 3) names
and contact information of three references. All attachments must be uploaded at the time of submission. The initial appointment will be for one year with the possibility of an extension depending on demonstrated progress and continued funding. For questions
regarding this position, please contact Dr. Maarten Buijsman by email: maarten.buijsman@usm.edu<file:///C:/Users/w944461/Documents/work/usm_admin/postdoc/TFO/maarten.buijsman@usm.edu><br>
<br>
<br>
<br>
The main goal of this position is to advance our understanding of (tidal) internal gravity waves in global HYbrid Coordinate Ocean Model (HYCOM) simulations: i.e., their generation, dispersion, interactions with the background flow, and dissipation. Stationary
internal tides ?dissipate? due to the scattering of energy to the nonstationary tide and to higher harmonics. The scattering to higher harmonics can be due to the nonlinear steepening of energetic internal tides. One objective is to determine how well this
process is simulated in hydrostatic global ocean model simulations with various horizontal and vertical resolutions.<br>
<br>
The successful candidate is expected to diagnose internal wave energetics in HYCOM simulations, validate these simulations with existing observations, assist in performing HYCOM simulations (if needed), collaborate with project scientists and students, write
first-author manuscripts for peer-review publications (~1/year) related to the ONR-TFO project and participate as co-author on other manuscripts, disseminate scientific results at national/international conferences and ONR meetings, contribute to annual report
preparation, as required by the funding agency, and perform other duties as assigned.<br>
<br>
<br>
<br>
Applicants must have a Ph.D. in physical oceanography or similar field awarded no later than May 31, 2020. Applicants must also have a publication record that illustrates ability to conduct and communicate independent research as well as experience and proficiency
in coding, performing (nonhydrostatic) numerical model simulations, time series analysis, internal wave processes, and ocean mixing. The successful candidate will be required to pass a NASA background security check to work at Stennis, and a USM employment
background check. The preferred candidate has experience in coding in Matlab and Fortran, analyzed HYCOM simulations, diagnosed internal wave energy equations, understanding of vertical mode decomposition, and performed multi-processor simulations on super
computers.<br>
<br>
<br>
<br>
The SOSE includes two academic divisions, Marine Science, and Coastal Sciences, and several R&D centers including: Hydrographic Science Research Center, Marine Research Center, and Thad Cochran Marine Aquaculture Center. The Division of Marine Science is based
at the NASA Stennis Space Center. John C. Stennis Space Center is a ?federal city? that boasts the world?s largest concentration of oceanographers and hydrographers. Marine Science faculty benefit from close working relationships with a number of on-site federal
agencies, including the Naval Research Laboratory-Stennis Space Center, the Naval Oceanographic Office, the Naval Meteorology and Oceanography Command, the USGS and NOAA, including the National Data Buoy Center, and the National Centers for Environmental Information.<br>
<br>
<br>
<br>
As an Affirmative Action/Equal Employment Opportunity employer/Americans with Disabilities Act institution, The University of Southern Mississippi encourages minorities, women, veterans and persons with disabilities to apply. For more information about Southern
Miss, visit our website at <a href="http://www.usm.edu">http://www.usm.edu</a>.<br>
<br>
<br>
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>
Maarten Buijsman (maarten.buijsman@usm.edu)<br>
Assistant Professor Physical Oceanography<br>
University of Southern Mississippi<br>
Division of Marine Science<br>
1020 Balch Blvd<br>
Stennis Space Center, MS 39529<br>
Phone: 228.688.2385 , Fax: 228.688.1121<br>
<a href="https://sites.google.com/site/maartenbuijsman/">https://sites.google.com/site/maartenbuijsman/</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200107/166ddcc2/attachment-0001.html">http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200107/166ddcc2/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 8 Jan 2020 02:30:25 +0000<br>
From: pan fei <feiamanda@hotmail.com><br>
To: "mitgcm-support-request@mitgcm.org"<br>
<mitgcm-support-request@mitgcm.org>, "mitgcm-support@mitgcm.org"<br>
<mitgcm-support@mitgcm.org><br>
Subject: [MITgcm-support] ??: MITgcm-support Digest, Vol 199, Issue 5<br>
Message-ID:<br>
<PS2PR06MB3464E864A3228082C39409B9A33E0@PS2PR06MB3464.apcprd06.prod.outlook.com><br>
<br>
Content-Type: text/plain; charset="gb2312"<br>
<br>
Hi Matt,<br>
<br>
I have checked the diagnosed wind stress, the wind stress pattern and value are OK.<br>
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).<br>
<br>
Here is my data.exf and EXF_OPTIONS.h.<br>
<br>
# *********************<br>
# External Forcing Data<br>
# *********************<br>
&EXF_NML_01<br>
exf_iprec = 32,<br>
exf_debugLev = 3,<br>
&<br>
<br>
# *********************<br>
&EXF_NML_02<br>
#hfluxstartdate1 = 19920101,<br>
#hfluxstartdate2 = 000000,<br>
hfluxperiod = -12.,<br>
#<br>
#sfluxstartdate1 = 19920101,<br>
#sfluxstartdate2 = 000000,<br>
sfluxperiod = -12.,<br>
#<br>
#ustressstartdate1 = 19920101,<br>
#ustressstartdate2 = 000000,<br>
ustressperiod = -12.,<br>
#<br>
#vstressstartdate1 = 19920101,<br>
#vstressstartdate2 = 000000,<br>
vstressperiod = -12.,<br>
#<br>
#climsststartdate1 = 19920101,<br>
#climsststartdate2 = 000000,<br>
climsstperiod = -12.,<br>
# 2 months restoring timescale for temperature<br>
climsstTauRelax = 5184000.,<br>
#<br>
#climsssstartdate1 = 19920101,<br>
#climsssstartdate2 = 000000,<br>
climsssperiod = -12.,<br>
# 6 months restoring timescale for salinity<br>
climsssTauRelax = 15552000.,<br>
#<br>
hfluxfile = 'qnet.bin',<br>
sfluxfile = 'emp.bin',<br>
# here we deviate from the general pattern and prescribe a forcing fiels<br>
# with (i,j)=(1,1) at the NW (.flipped) instead of the SW (.bin) corner<br>
# to test the correpsonding code, see also below<br>
ustressfile = 'taux.bin',<br>
vstressfile = 'tauy.bin',<br>
climsstfile = 'SST.bin',<br>
climsssfile = 'SSS.bin',<br>
#<br>
&<br>
<br>
# *********************<br>
&EXF_NML_03<br>
&<br>
<br>
# *********************<br>
&EXF_NML_04<br>
hflux_lon0 = 0.,<br>
hflux_lon_inc = 1.,<br>
hflux_lat0 = -80.,<br>
hflux_lat_inc = 159*1.,<br>
hflux_nlon = 360,<br>
hflux_nlat = 160,<br>
#<br>
sflux_lon0 = 0.,<br>
sflux_lon_inc = 1.,<br>
sflux_lat0 = -80.,<br>
sflux_lat_inc = 159*1.,<br>
sflux_nlon = 360,<br>
sflux_nlat = 160,<br>
#<br>
ustress_lon0 = 0.,<br>
ustress_lon_inc = 1.,<br>
ustress_lat0 = -80.,<br>
ustress_lat_inc = 159*1.,<br>
ustress_nlon = 360,<br>
ustress_nlat = 160,<br>
#<br>
vstress_lon0 = 0.,<br>
vstress_lon_inc = 1.,<br>
vstress_lat0 = -80.,<br>
vstress_lat_inc = 159*1.,<br>
vstress_nlon = 360,<br>
vstress_nlat = 160,<br>
#<br>
climsst_lon0 = 0.,<br>
climsst_lon_inc = 1.,<br>
climsst_lat0 = -80.,<br>
climsst_lat_inc = 159*1.,<br>
climsst_nlon = 360,<br>
climsst_nlat = 160,<br>
#<br>
climsss_lon0 = 0.,<br>
climsss_lon_inc = 1.,<br>
climsss_lat0 = -80.,<br>
climsss_lat_inc = 159*1.,<br>
climsss_nlon = 360,<br>
climsss_nlat = 160,<br>
#<br>
&<br>
<br>
<br>
<br>
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
<br>
<br>
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>
C $Name: checkpoint65o $<br>
<br>
CBOP<br>
C !ROUTINE: EXF_OPTIONS.h<br>
C !INTERFACE:<br>
C #include "EXF_OPTIONS.h"<br>
<br>
C !DESCRIPTION:<br>
C *==================================================================*<br>
C | CPP options file for EXternal Forcing (EXF) package:<br>
C | Control which optional features to compile in this package code.<br>
C *==================================================================*<br>
CEOP<br>
<br>
#ifndef EXF_OPTIONS_H<br>
#define EXF_OPTIONS_H<br>
#include "PACKAGES_CONFIG.h"<br>
#include "CPP_OPTIONS.h"<br>
<br>
#ifdef ALLOW_EXF<br>
#ifdef ECCO_CPPOPTIONS_H<br>
<br>
C-- When multi-package option-file ECCO_CPPOPTIONS.h is used (directly included<br>
C in CPP_OPTIONS.h), this option file is left empty since all options that<br>
C are specific to this package are assumed to be set in ECCO_CPPOPTIONS.h<br>
<br>
#else /* ndef ECCO_CPPOPTIONS_H */<br>
<br>
C-- Package-specific Options & Macros go here<br>
<br>
C pkg/exf CPP options:<br>
C --------------------<br>
C<br>
C > ( EXF_VERBOSE ) < replaced with run-time, logical parameter "exf_verbose".<br>
C<br>
C >>> ALLOW_ATM_WIND <<<<br>
C If defined, 10-m wind fields can be read-in from files.<br>
C<br>
C >>> ALLOW_ATM_TEMP <<<<br>
C If defined, atmospheric temperature and specific<br>
C humidity fields can be read-in from files.<br>
C<br>
C >>> ALLOW_DOWNWARD_RADIATION <<<<br>
C If defined, downward long-wave and short-wave radiation<br>
C can be read-in form files or computed from lwflux and swflux.<br>
C<br>
C >>> ALLOW_ZENITHANGLE <<<<br>
C If defined, ocean albedo varies with the zenith angle, and<br>
C incoming fluxes at the top of the atmosphere are computed<br>
C<br>
C >>> ALLOW_BULKFORMULAE <<<<br>
C Allows the use of bulk formulae in order to estimate<br>
C turbulent and radiative fluxes at the ocean surface.<br>
C<br>
C >>> EXF_READ_EVAP <<<<br>
C If defined, evaporation fields are read-in, rather than<br>
C computed from atmospheric state.<br>
C<br>
C >>> ALLOW_RUNOFF <<<<br>
C If defined, river and glacier runoff can be read-in from files.<br>
C<br>
C >>> ALLOW_RUNOFTEMP <<<<br>
C If defined, river and glacier runoff temperature<br>
C can be read-in from files.<br>
C<br>
C >>> ATMOSPHERIC_LOADING <<<<br>
C If defined, atmospheric pressure can be read-in from files.<br>
C WARNING: this flag is set (define/undef) in CPP_OPTIONS.h<br>
C and cannot be changed here (in EXF_OPTIONS)<br>
C<br>
C >>> EXF_SEAICE_FRACTION <<<<br>
C If defined, seaice fraction can be read-in from files (areaMaskFile)<br>
C<br>
C >>> ALLOW_CLIMSST_RELAXATION <<<<br>
C Allow the relaxation to a monthly climatology of sea surface<br>
C temperature, e.g. the Reynolds climatology.<br>
C<br>
C >>> ALLOW_CLIMSSS_RELAXATION <<<<br>
C Allow the relaxation to a monthly climatology of sea surface<br>
C salinity, e.g. the Levitus climatology.<br>
C<br>
C >>> USE_EXF_INTERPOLATION <<<<br>
C Allows specification of arbitrary Cartesian input grids.<br>
C<br>
C ====================================================================<br>
C<br>
C The following CPP options:<br>
C<br>
C ALLOW_ATM_WIND (WIND)<br>
C ALLOW_ATM_TEMP (TEMP)<br>
C ALLOW_DOWNWARD_RADIATION (DOWN)<br>
C ALLOW_BULKFORMULAE (BULK)<br>
C EXF_READ_EVAP (EVAP)<br>
C<br>
C permit the ocean-model forcing configurations listed in the<br>
C table below. The first configuration is the default,<br>
C flux-forced, ocean model. The next four are stand-alone<br>
C configurations that use pkg/exf, open-water bulk formulae to<br>
C compute the missing surface fluxes from atmospheric variables.<br>
C The last four configurations can be used in conjunction with<br>
C pkg/seaice to model ice-covered regions. The forcing fields<br>
C in the rightmost column are defined in exf_fields.<br>
C<br>
C<br>
C WIND |TEMP |DOWN |BULK |EVAP | actions<br>
C -----|-----|-----|-----|-----|-------------------------------------<br>
C | | | | |<br>
C - | - | - | - | - | Read-in ustress, vstress, hflux,<br>
C | | | | | swflux, and sflux.<br>
C | | | | |<br>
C def | def | def | def | - | Read-in uwind, vwind, atemp, aqh,<br>
C | | | | | swdown, lwdown, precip, and runoff.<br>
C | | | | | Compute ustress, vstress, hflux,<br>
C | | | | | swflux, and sflux.<br>
C | | | | |<br>
C def | def | - | def | - | Read-in uwind, vwind, atemp, aqh,<br>
C | | | | | swflux, lwflux, precip, and runoff.<br>
C | | | | | Compute ustress, vstress, hflux,<br>
C | | | | | and sflux.<br>
C | | | | |<br>
C def | - | - | def | - | Read-in uwind, vwind, hflux,<br>
C | | | | | swflux, and sflux.<br>
C | | | | | Compute ustress and vstress.<br>
C | | | | |<br>
C - | def | - | def | - | Read-in ustress, vstress, atemp,<br>
C | | | | | aqh, swflux, lwflux, precip, and<br>
C | | | | | runoff. Compute hflux and sflux.<br>
C | | | | |<br>
C def | def | - | - | def | Read-in uwind, vwind, atemp, aqh,<br>
C | | | | | swflux, lwflux, precip, runoff,<br>
C | | | | | and evap.<br>
C | | | | |<br>
C def | def | - | def | - | Read-in uwind, vwind, atemp, aqh,<br>
C | | | | | swflux, lwflux, precip, and runoff.<br>
C | | | | | Compute open-water ustress, vstress,<br>
C | | | | | hflux, swflux, and evap.<br>
C | | | | |<br>
C def | def | def | - | def | Read-in uwind, vwind, atemp, aqh,<br>
C | | | | | swdown, lwdown, precip, runoff,<br>
C | | | | | and evap.<br>
C | | | | |<br>
C def | def | def | def | - | Read-in uwind, vwind, atemp, aqh,<br>
C | | | | | swdown, lwdown, precip, and runoff.<br>
C | | | | | Compute open-water ustress, vstress,<br>
C | | | | | hflux, swflux, and evap.<br>
C<br>
C ====================================================================<br>
<br>
C Bulk formulae related flags.<br>
#undef ALLOW_ATM_TEMP<br>
#undef ALLOW_ATM_WIND<br>
#undef ALLOW_DOWNWARD_RADIATION<br>
#undef ALLOW_RUNOFF<br>
#undef ALLOW_RUNOFTEMP<br>
#if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))<br>
# define ALLOW_BULKFORMULAE<br>
# undef ALLOW_BULK_LARGEYEAGER04<br>
#endif<br>
<br>
C Zenith Angle/Albedo related flags.<br>
#ifdef ALLOW_DOWNWARD_RADIATION<br>
# undef ALLOW_ZENITHANGLE<br>
#endif<br>
<br>
C Use ocean_emissivity*lwdown in lwFlux. This flag should be defined<br>
C unless to reproduce old results (obtained with inconsistent old code)<br>
#ifdef ALLOW_DOWNWARD_RADIATION<br>
# define EXF_LWDOWN_WITH_EMISSIVITY<br>
#endif<br>
<br>
C Relaxation to monthly climatologies.<br>
#define ALLOW_CLIMSST_RELAXATION<br>
#define ALLOW_CLIMSSS_RELAXATION<br>
<br>
C Allows to read-in seaice fraction from files (areaMaskFile)<br>
#undef EXF_SEAICE_FRACTION<br>
<br>
C Use spatial interpolation to interpolate<br>
C forcing files from input grid to model grid.<br>
#define USE_EXF_INTERPOLATION<br>
C for interpolated vector fields, rotate towards model-grid axis<br>
C using old rotation formulae (instead of grid-angles)<br>
#undef EXF_USE_OLD_VEC_ROTATION<br>
C for interpolation around N & S pole, use the old formulation<br>
C (no pole symmetry, single vector-comp interp, reset to 0 zonal-comp @ N.pole)<br>
#undef EXF_USE_OLD_INTERP_POLE<br>
<br>
#undef EXF_INTERP_USE_DYNALLOC<br>
#if ( defined (EXF_INTERP_USE_DYNALLOC) && defined (USING_THREADS) )<br>
# define EXF_IREAD_USE_GLOBAL_POINTER<br>
#endif<br>
<br>
#endif /* ndef ECCO_CPPOPTIONS_H */<br>
#endif /* ALLOW_EXF */<br>
#endif /* EXF_OPTIONS_H */<br>
<br>
<br>
Thank you so much.<br>
<br>
Best,<br>
<br>
Fei<br>
<br>
________________________________<br>
???: MITgcm-support <mitgcm-support-bounces@mitgcm.org> ?? mitgcm-support-request@mitgcm.org <mitgcm-support-request@mitgcm.org><br>
????: 2020?1?8? 1:00<br>
???: mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
??: MITgcm-support Digest, Vol 199, Issue 5<br>
<br>
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>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200108/c22db2da/attachment.html">http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200108/c22db2da/attachment.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 6<br>
**********************************************<br>
</div>
</span></font></div>
<span>_______________________________________________</span><br><span>MITgcm-support mailing list</span><br><span>MITgcm-support@mitgcm.org</span><br><span>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</span><br></div></blockquote></div></body></html>