[MITgcm-support] EXF files from ncep R1
Santha Akella
santha.akella at gmail.com
Sat Jan 9 20:09:45 EST 2010
Dear MITgcm Users,
I am trying to use the configuration in
verification/natl_box_adjoint
I am trying to figure out the correct files to go in the data.exf using ncep
R1 surface fluxes from:
http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.surfaceflux.html
1. atempfile:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=Air+Temperature
2. aqhfile:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=3&tid=25330&vid=29
(converted to Specific Humidity in Kg/Kg)
3. uwind :
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=u-wind&Level=10
4. vwind:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=v-wind&Level=10
5. precipfile:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=Precipitation+rate
(converted to m/s)
6. lwdownfile:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=Downward+longwave+radiation+flux
7. swdownfile:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=Downward+solar+radiation+flux
Also, I have attached CPP_OPTIONS.h and EXF_OPTIONS.h that I am using.
Did I get the correct files ? Or alternatively, if I was to force with
hflux, sflux, atemp, aqh, uwind, vwind, precip, lwflux & swflux, which files
would I need to calculate hflux & sflux? I guess for
--lwflux:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=4&tid=25331&vid=1261
--swflux:
http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?Dataset=NCEP+Reanalysis+Surface+Flux&Variable=Upward+solar+radiation+flux
For the above, I tried to get some hints from
verification/tutorial_global_oce_latlon/diags_matlab/ncep2global_ocean.m
but I don't know how to get ncep .cdf files used.
Thanks a lot for the help!
Santha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100109/84cf1950/attachment.htm>
-------------- next part --------------
C $Header: /u/gcmpack/MITgcm/verification/natl_box_adjoint/code_bulk/CPP_OPTIONS.h,v 1.2 2003/11/06 22:30:59 heimbach Exp $
C
C
C CPP flags controlling which code in included in the files that
C will be compiled.
#ifndef CPP_OPTIONS_H
#define CPP_OPTIONS_H
#include "PACKAGES_CONFIG.h"
C o Shortwave heating as extra term in external_forcing.F
#ifdef ALLOW_KPP
#define SHORTWAVE_HEATING
#endif
C o Include/exclude phi_hyd calculation code
#define INCLUDE_PHIHYD_CALCULATION_CODE
C o Include/exclude call to S/R FIND_RHO
#define INCLUDE_FIND_RHO_CALL
C o Include/exclude call to S/R CONVECT
#define INCLUDE_CONVECT_CALL
C o Include/exclude call to S/R CALC_ISOSLOPES
#define INCLUDE_CALC_ISOSLOPES_CALL
C o Include/exclude call to S/R CALC_DIFFUSIVITY
#define INCLUDE_CALC_DIFFUSIVITY_CALL
C o Allow nonHydrostatic code
#undef ALLOW_NONHYDROSTATIC
C o Use "natural" boundary conditions for salinity
C instead of the "virtual salt flux"
#undef USE_NATURAL_BCS
C o Include/exclude exact volume conservaion
#define EXACT_CONSERV
C o Execution environment support options
#include "CPP_EEOPTIONS.h"
C o Include/exclude code specific to the ECCO/SEALION version.
#ifdef ALLOW_ECCO
# include "ECCO_CPPOPTIONS.h"
#endif
#endif /* CPP_OPTIONS_H */
-------------- next part --------------
C $Header: /u/gcmpack/MITgcm/verification/lab_sea/code/EXF_OPTIONS.h,v 1.2 2007/05/10 22:31:12 jmc Exp $
C $Name: $
#ifndef EXF_OPTIONS_H
#define EXF_OPTIONS_H
#include "PACKAGES_CONFIG.h"
#ifdef ALLOW_EXF
#include "CPP_OPTIONS.h"
#ifdef ALLOW_AUTODIFF_TAMC
C When compile for AD mode,
C this package is currently configured by ECCO_CPPOPTIONS.h
C which is directly included in CPP_OPTIONS.h
#else /* ndef ALLOW_AUTODIFF_TAMC */
C CPP flags controlling which code is included in the files that
C will be compiled.
C
c pkg/exf CPP options:
c --------------------
c
c >>> INCLUDE_EXTERNAL_FORCING_PACKAGE <<<
c Include this package into the setup.
c
c >>> EXF_VERBOSE <<<
c Do a bit more printout for the log file than usual.
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_BULKFORMULAE <<<
c Allows the use of bulk formulae in order to estimate
c turbulent and radiative fluxes at the ocean's 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 >>> 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 >>> 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 Do more printout for the protocol file than usual.
#undef EXF_VERBOSE
C Bulk formulae related flags:
C those options are required to use pkg/exf with pkg/seaice.
#define ALLOW_ATM_TEMP
#define ALLOW_ATM_WIND
#define ALLOW_DOWNWARD_RADIATION
#define ALLOW_RUNOFF
#define ALLOW_BULKFORMULAE
#undef ALLOW_BULK_LARGEYEAGER04
C Relaxation to monthly climatologies.
#define ALLOW_CLIMSST_RELAXATION
#define ALLOW_CLIMSSS_RELAXATION
C Use spatial interpolation to interpolate
C forcing files from input grid to model grid.
#undef USE_EXF_INTERPOLATION
#define EXF_INTERP_USE_DYNALLOC
#if ( defined (EXF_INTERP_USE_DYNALLOC) & defined (USING_THREADS) )
# define EXF_IREAD_USE_GLOBAL_POINTER
#endif
#endif /* ndef ALLOW_AUTODIFF_TAMC */
#endif /* ALLOW_EXF */
#endif /* EXF_OPTIONS_H */
More information about the MITgcm-support
mailing list