[Mitgcm-support] 2-deg North Atlantic Test Case
mitgcm-support at dev.mitgcm.org
mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:56:31 EDT 2003
HyperNews notification of new message. See:
http://escher.JPL.NASA.GOV:2000/HyperNews/get/forums/model/19.html
I am appending the readme file "NAtlExp/README.natl" from
a North Atlantic test configuration. Differences between
our version of the code and Classic are being investigated.
North Atlantic test
===================
A North Atlantic test configuration has been prepared based on Patrick
Heimbach's North Atlantic MIT GCM Classic integrations. Test domain consists of
a 20x16x23 box, forcing files are from NCEP re-analysis, and integration
period is one year. This configuration is intended to fill the much-needed
gap of cross-platform cross-version testing. Of immediate interest are
updates to KPP and GM schemes. The test code presently resides in a tarball
on "triton.jpl.nasa.gov:/dm1/dimitri/mitgcm/NATLANtest/NAtlExp.tar.gz".
Description of input and output files for various configurations and
instructions for setting up and executing test integrations follow.
Subdirectory NAtlExp/matlab contains scripts to create forcing files for
various code versions and for comparing output.
Domain Configuration
====================
lon=322:2:360;
lat=13:2:43;
dZ=[10 10 15 20 20 25 35 50 75 100 150 200 275 ...
350 415 450 500 500 500 500 500 500 500];
Classic MIT GCM
===============
Code files:
-----------
cfg.closed Model parameter file
Input files:
------------
PH_pmask_020_016_23.dat Land mask (0=dry, 1=wet)
PH_LevCli.ptmp_020_016_23.dat Levitus temperature climatology (deg C)
PH_LevCli.salt_020_016_23.dat Levitus salt climatology (PSU)
Format is real*4 plain binary 20*16*23
PH_SST_fields_020_016_23.dat Sea-surface temperature (deg C)
Format is real*4 plain binary 20*16*12; t=0:30:330 days
PH_totflsfc_020_016_23.dat Heat flux less shortwave; >0 = cooling (W/m^2)
PH_swflsfc_020_016_23.dat Shortwave radiation; >0 = cooling (W/m^2)
PH_pratesfc_020_016_23.dat Surface salt flux; >0 = evaporation (PSU/s)
Format is real*4 plain binary 20*16*367; t=-.75:1:356.25 days
PH_ustress_020_016_23.dat Zonal wind stress (N/m^2)
PH_vstress_020_016_23.dat Meridional wind stress (N/m^2)
Format is real*4 plain binary 20*16*734; t=-.875:.5:365.625 days
Output files:
-------------
PH_hfield_000.dat SSH (m)
Format is real*4 plain binary 20*16*365; t=0:1:364 days
PH_kmfield_000.dat KPP vertical viscosity (m^2/s)
PH_ksfield_000.dat KPP vertical diffusivity (m^2/s)
PH_sfield_000.dat salinity (PSU)
PH_tfield_000.dat temperature (deg C)
PH_ufield_000.dat zonal velocity (m/s)
PH_vfield_000.dat meridional velocity (m/s)
Format is real*4 plain binary 20*16*23*365; t=0:1:364 days
JPL mitgcmuv
============
Code files:
-----------
CPP_OPTIONS.h Compiler options file
SIZE.h Domain definition file
Input files:
------------
data Model parameter file
data.kpp KPP parameter file
POLY3.COEFFS Equation of state coefficients
Forcing files are based on Classic input files and were created using matlab
script "NAtlExp/matlab/mk_input.m".
topog Bathymetry (m)
Format is real*4 FORTRAN unformatted sequential binary 20*16
saltLev Levitus salt climatology (PSU)
thetaLev Levitus temperature climatology (deg C)
Format is real*4 FORTRAN unformatted sequential binary 20*16*23
SST Sea-surface temperature (deg C)
dQdT Thermal relaxation coefficient (W/m^2/deg)
SSS Sea-surface salinity (PSU)
Format is real*4 JPL direct access (20*16*1+1)*14; t=0:30:390 days
qnet Heat flux less shortwave; >0 = cooling (W/m^2)
qsw Shortwave radiation; >0 = cooling (W/m^2)
EmPmR Surface salt flux; >0 = evaporation (PSU/s)
Format is real*4 JPL direct access (20*16*1+1)*367; t=-.75:1:356.25 days
windx Zonal wind stress (N/m^2)
windy Meridional wind stress (N/m^2)
Format is real*4 JPL direct access (20*16*1+1)*734; t=-.875:.5:365.625 days
Output files:
-------------
Output files are in JPL format. They can be read into 4-D arrays (lon, lat,
depth, time) using matlab function "NAtlExp/matlab/readgcm.m".
H_01_01.00000_08760_024 SSH (m)
KPPhbl_01_01.00000_08760_024 KPP OBL depth (m)
read using: [H t]=readgcm('H_01_01.00000_08760_024',0,20,16,1);
KPPvisc_01_01.00000_08760_024 KPP vertical viscosity (m^2/s)
KPPdiff_01_01.00000_08760_024 KPP vertical diffusivity (m^2/s)
KPPghat_01_01.00000_08760_024 KPP non-local transport (s/m^2)
S_01_01.00000_08760_024 salinity (PSU)
T_01_01.00000_08760_024 temperature (deg C)
U_01_01.00000_08760_024 zonal velocity (m/s)
V_01_01.00000_08760_024 meridional velocity (m/s)
read using: [V t]=readgcm('V_01_01.00000_08760_024',0,20,16,23);
Instructions:
-------------
Instructions for running JPL mitgcmuv test follow (assuming that NAtlExp.tar.gz
has been copied into mitgcmuv execution directory):
cd mitgcmuv/exe
gunzip NAtlExp.tar.gz
tar xvf NAtlExp.tar
cp NAtlExp/jpl/code/KPP_OPTIONS.h ../pkg/kpp
cp NAtlExp/jpl/code/CPP_OPTIONS.h ../model/inc
cp NAtlExp/jpl/code/SIZE.h ../model/inc
c ../bin
../tools/genmake -makefile
make depend
make
cd ../exe
ln -sf NAtlExp/jpl/input/POLY3.COEFFS .
ln -sf NAtlExp/jpl/input/data.kpp .
ln -sf NAtlExp/jpl/input/data .
ln -sf ../verification/exp5/input/eedata .
mitgcmuv > stdout
c30 mitgcmuv
============
Code files:
-----------
CPP_OPTIONS.h Compiler options file
genmake genmake file compatible with that of JPL code
Input files:
------------
Forcing files are same as those in Classic/Input except for bathymetry which
was created using "NAtlExp/matlab/mk_input.m".
topog Bathymetry (m)
Format is real*4 plain binary 20*16
Instructions:
-------------
Instructions for running c30 mitgcmuv test follow (assuming that
NAtlExp.tar.gz has been copied into mitgcmuv execution directory).
(This section is still under construction.)
cd mitgcmuv/exe
gunzip NAtlExp.tar.gz
tar xvf NAtlExp.tar
cp NAtlExp/c30/code/CPP_OPTIONS.h ../model/inc
cp NAtlExp/jpl/code/SIZE.h ../model/inc
c ../bin
ln -sf ../exe/NAtlExp/c30/code/genmake ../tools
../tools/genmake -makefile
make depend
make
cd ../exe
ln -sf NAtlExp/jpl/POLY3.COEFFS .
ln -sf NAtlExp/jpl/data.kpp .
ln -sf ../verification/exp2/inputNAtlExp/jpl/eedata .
ln -sf NAtlExp/c30/input/data.pkg .
ln -sf NAtlExp/c30/input/data .
mitgcmuv > stdout
More information about the MITgcm-support
mailing list