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