[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:54 EDT 2003


HyperNews notification of new message.  See: 
  http://escher.JPL.NASA.GOV:2000/HyperNews/get/forums/model/19/1.html

2-deg North Atlantic Test Case
==============================

The North Atlantic test case is now configured and running
for three flavors of the MIT GCM: MIT's c31, JPL's mitgcmuv,
and JPL's ecco.  All three codes use exactly the same KPP
package and give qualitatively similar results, although
differences between c31 and JPL, order 10^-3 deg C at surface
after a month of integration, are larger than differences
between the two JPL codes (order 10^-6 deg C).

Files and instructions for configuring and running the tests
and for comparing the various results are in
triton.jpl.nasa.gov:/dm1/dimitri/mitgcm/NATLANtest

Directory is also permanently stored in
nebula.jpl.nasa.gov:/silo/dmenem/triton.mitgcm/NATLANtest/NATLANtest_3oct00.tar

Instructions for configuring tests and comparing results are
in NATLANtest/NAtlExp/README.natl which is appended below.

>>>>>>>>>>>>> NATLANtest/NAtlExp/README.natl <<<<<<<<<<<<<<<<

North Atlantic test
===================

This test configuration is intended as a cross-platform, cross-version test of
KPP mixing scheme.  Test domain consists of a closed 20x16x23 box.  Forcing
files are daily for a period of 1 month.  Description of input and output
files for various model versions and instructions for setting up and executing
test integrations follow.

File "NAtlExp/matlab/comp_jpl_c31.m" is an example routine for reading and
comparing c31 and JPL output.

File "NAtlExp/matlab/comp_jpl.m" is an example routine for reading and
comparing JPL mitgcmuv and JPL ecco 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];

JPL c19_gm13
============

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

jpl_topog      Bathymetry (m)
   Format is real*4 FORTRAN unformatted sequential binary 20*16

jpl_saltLev    Levitus salt climatology (PSU)
jpl_thetaLev   Levitus temperature climatology (deg C)
   Format is real*4 FORTRAN unformatted sequential binary 20*16*23

jpl_windx      Zonal wind stress; >0 from West to East (N/m^2)
jpl_windy      Meridional wind stress; >0 from South to North (N/m^2)
jpl_qnet       Heat flux less shortwave; >0 = cooling (W/m^2)
jpl_sw         Shortwave radiation; >0 = cooling (W/m^2)
jpl_ep         Surface salt flux; >0 = evaporation (kg/m^2/s)
jpl_sst        Sea-surface temperature (deg C)
jpl_dqdt       Thermal relaxation coefficient (W/m^2/deg)
jpl_sss        Sea-surface salinity (PSU)
   Format is real*4 JPL direct access (20*16*1+1)*32; t=0:1:31 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_00720_240          SSH                      (m)
KPPhbl_01_01.00000_00720_240     KPP OBL depth            (m)
   read using: [H t]=readgcm('H_01_01.00000_00720_240',0,20,16,1);

KPPvisc_01_01.00000_00720_240    KPP vertical viscosity   (m^2/s)
KPPdiff_01_01.00000_00720_240    KPP vertical diffusivity (m^2/s)
KPPghat_01_01.00000_00720_240    KPP non-local transport  (s/m^2)
S_01_01.00000_00720_240          salinity                 (PSU) 
T_01_01.00000_00720_240          temperature              (deg C)
U_01_01.00000_00720_240          zonal velocity           (m/s)
V_01_01.00000_00720_240          meridional velocity      (m/s)
   read using: [V t]=readgcm('V_01_01.00000_00720_240',0,20,16,23);

Instructions:
-------------

Instructions for running JPL mitgcmuv test follow.

1) Expand tgz files

   gunzip NAtlExp.tar.gz
   tar xvf NAtlExp.tar
   gunzip c19_gm13.tar.gz
   tar xvf c19_gm13.tar

   Note: c19_gm13.tar can also be obtained from JPL CVS tree
   cvs export -r c19_gm13 mitgcmuv

2) Set-up compile, link, and execute code:

   cd mitgcmuv
   mkdir exe
   mkdir bin
   cp ../NAtlExp/jpl/code/*.h model/inc
   cp ../NAtlExp/jpl/code/*.F model/src
   cd bin
   ../tools/genmake -makefile
   make depend
   make
   cd ../exe
   ln -sf ../../NAtlExp/jpl/input/POLY3.COEFFS .
   ln -sf ../../NAtlExp/jpl/input/data* .
   ln -sf ../verification/exp5/input/eedata .
   mitgcmuv > stdout

JPL ECCO kpp03
==============

Code files:
-----------

CPP_OPTIONS.h  Compiler options file
SIZE.h         Domain definition file
ini_forcing.F
load_external_fields.F

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_00720_240          SSH                      (m)
KPPhbl_01_01.00000_00720_240     KPP OBL depth            (m)
   read using: [H t]=readgcm('H_01_01.00000_00720_240',0,20,16,1);

KPPvisc_01_01.00000_00720_240    KPP vertical viscosity   (m^2/s)
KPPdiff_01_01.00000_00720_240    KPP vertical diffusivity (m^2/s)
KPPghat_01_01.00000_00720_240    KPP non-local transport  (s/m^2)
S_01_01.00000_00720_240          salinity                 (PSU) 
T_01_01.00000_00720_240          temperature              (deg C)
U_01_01.00000_00720_240          zonal velocity           (m/s)
V_01_01.00000_00720_240          meridional velocity      (m/s)
   read using: [V t]=readgcm('V_01_01.00000_00720_240',0,20,16,23);

Instructions:
-------------

Instructions for running JPL ecco test follow.

1) Expand tgz files and get ecco code from JPL repository

   gunzip NAtlExp.tar.gz
   tar xvf NAtlExp.tar
   cvs export -r kpp03 ecco

2) Set-up compile, link, and execute code:

   cd ecco
   cp verification/natlan/code/*.h src/includes
   cp verification/natlan/code/*.F src/model
   touch -r acinclude.m4 configure.in configure
   setenv CC cc
   setenv FC f90
   configure
   cd lapack
   make
   cd ..
   make
   cd src
   f90 -mp -mpio -o ecco -L../blas -L../lapack main.o \
   the_model_main.o model/libmodel.a model/find_alpha.o \
   solver/libsolver.a diags/libdiags.a wrapper/libwrapper.a \
   kpp/libkpp.a io/libio.a -llapack -lblas
   cd ..
   mkdir exe
   cd exe
   cp ../src/ecco .
   ln -sf ../../NAtlExp/jpl/input/POLY3.COEFFS .
   ln -sf ../../NAtlExp/jpl/input/data* .
   ln -sf ../verification/exp5/input/eedata .
   ecco > stdout

c31 mitgcmuv
============

Code files:
-----------

Input files:
------------

data           Model parameter file
data.pkg       Package parameter file
data.kpp       KPP parameter file
POLY3.COEFFS   Equation of state coefficients

c31_topog      Bathymetry (m)
   Format is real*4 plain binary 20*16

c31_saltLev    Levitus salt climatology (PSU)
c31_thetaLev   Levitus temperature climatology (deg C)
   Format is real*4 plain binary 20*16*23

c31_windx      Zonal wind stress; >0 from East to West (N/m^2)
c31_windy      Meridional wind stress; >0 from North to South (N/m^2)
c31_qnet       Heat flux less shortwave; >0 = cooling (W/m^2)
c31_sw         Shortwave radiation; >0 = cooling (W/m^2)
c31_ep         Surface salt flux; >0 = evaporation (m/s)
c31_sst        Sea-surface temperature (deg C)
c31_sss        Sea-surface salinity (PSU)
   Format is real*4 plain binary 20*16*32; t=0:1:31 days

Output files:
-------------

Output files are in direct access binary format.  Read to 4-D arrays
(lon, lat, depth, time) using matlab function "NAtlExp/matlab/readbin.m".

PS.001.001.data          SSH                      (m)
KPPhbl.001.001.data     KPP OBL depth            (m)
   read using: H=readbin('KPPhbl.001.001.data',[20 16 1 4],1);

KPPviscAz.001.001.data    KPP vertical viscosity   (m^2/s)
KPPdiffKzT.001.001.data    KPP vertical diffusivity (m^2/s)
KPPghat.001.001.data    KPP non-local transport  (s/m^2)
S.001.001.data          salinity                 (PSU) 
T.001.001.data          temperature              (deg C)
U.001.001.data          zonal velocity           (m/s)
V.001.001.data          meridional velocity      (m/s)
   read using: T=readbin('T.001.001.data',[20 16 23 4],1);

Instructions:
-------------

Instructions for running c31 test follow.

1) Expand tgz files

   gunzip NAtlExp.tar.gz
   tar xvf NAtlExp.tar
   gunzip c31_092600.tar.gz
   tar xvf c31_092600.tar

   Note: c31_092600.tar can also be obtained from MIT CVS tree
   mitcvs checkout -D "09/26/2000" models/MITgcmUV

2) Set-up compile, link, and execute code:

   cd NAtlExp/c31/code
   cp genmake ../../../MITgcmUV/tools
   cp kpp_do_diags.F ../../../MITgcmUV/pkg/kpp
   cp *.h ../../../MITgcmUV/model/inc
   cp e*.F i*.F s*.F w*.F ../../../MITgcmUV/model/src
   cd ../../../MITgcmUV/bin
   ../tools/genmake -makefile
   make depend
   make
   cd ../exe
   ln -sf ../../NAtlExp/c31/input/POLY3.COEFFS .
   ln -sf ../../NAtlExp/c31/input/*data* .
   mitgcmuv > stdout





More information about the MITgcm-support mailing list