noaa-gfdl / am4 Goto Github PK

The input tarball has been truncated to about 2GiB and is missing quite a few files. It's possible to download the individual files from the AM4_run directory but it would be nicer to have the compress file available.

ftp://nomads.gfdl.noaa.gov/users/Ming.Zhao/AM4Documentation/GFDL-AM4.0/inputData

Hello,

I'm having trouble running a model.
It's an error about land_frac,
Can I know how to solve it?
Below is the error message.

FATAL from PE 35: diag_manager_mod::register_static_field: module/output_field land/land_frac AREA measures field NOT found in diag_table. Contact the model liaison.n
[cli_35]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 35
FATAL from PE 26: diag_manager_mod::register_static_field: module/output_field land/land_frac AREA measures field NOT found in diag_table. Contact the model liaison.n
….
(skip)

Thank you.

Adding additional text to the README.md to better explain how to obtain the input data files would be useful, and my remove some user confusion.

The make templates cannot determine how to find the HDF5 Fortran module, nor the libraries. It would be good if the templates could attempt to guess these locations, or give the user the ability to add the required flags easier.

In order to work on the fv_mapz test case, we would need to get AM4 running with GNU rather than intel.

There is a problem in obtaining the input data. I can not download the data even though I had closed the firewall。
[liguanghe@localhost ~]$ wget ftp://nomads.gfdl.noaa.gov/users/Ming.Zhao/AM4Documentation/GFDL-AM4.0/inputData/AM4_run.tar.gz
--2022-03-19 18:46:31-- ftp://nomads.gfdl.noaa.gov/users/Ming.Zhao/AM4Documentation/GFDL-AM4.0/inputData/AM4_run.tar.gz
=> 'AM4_run.tar.gz'
Resolving nomads.gfdl.noaa.gov (nomads.gfdl.noaa.gov)... 2610:20:8020:6c01::31:31, 140.208.31.31
Connecting to nomads.gfdl.noaa.gov (nomads.gfdl.noaa.gov)|2610:20:8020:6c01::31:31|:21... connected.
Logging in as anonymous ... Logged in!
==> SYST ... done. ==> PWD ... done.
==> TYPE I ... done. ==> CWD (1) /users/Ming.Zhao/AM4Documentation/GFDL-AM4.0/inputData ... done.
==> SIZE AM4_run.tar.gz ... 22268830166
==> EPSV ... ==> LPSV ...
Cannot initiate PASV transfer.
==> EPRT ...

Hi there -

The variable n_sponge is set to 0 by default, which would suggest there is no sponge layer. This doesn't seem correct: there should be some sponge layer, but I'm not sure where it's defined. Should I assume the sponge layer extends down to 8 hPa like in Zhao et al. 2018?

To introduce small perturbations to our initial conditions, I used the variable add_noise. I was wondering at which vertical levels is thermal noise added? Is it only to the lowest levels, or is it to all levels?

Thanks for any help,
Marshall

The run instructions need to be updated to include copying the updated ascii files into the run directory, per #28

The FMScoupler project has a dependency on missing land_model_mod.mod. See NOAA-GFDL/FMScoupler#15.

So I am trying to build land_model_mod.mod so that I can get the autotools build for FMScoupler working.

However, when I build AM4 according to your home rolled build system (a.k.a. your "workflow") it does not build land_model_mod.mod.

find ~/AM4/ -name 'land_model*'
/home/Edward.Hartnett/AM4/src/land_lad2/land_model.F90

How do I get AM4, using your legacy build system, to build land_model_mod.mod?

Thanks!

Hello,

I am trying to do an ensemble simulation using the AM4.
When conducting ensemble simulations, how do you change the initial conditions?
Is it correct to reconstruct initial condition variables such as U, V, T, etc., in the form of adding or subtracting small random values for each point.

In coupler/Makefile, the src file of atmos_ocean_dep_fluxes_calc.o is $(SRCROOT)FMScoupler/full/atmos_ocean_dep_fluxes_calc.f90, but the src file's postfix is F90, not f90. It would throw errors when I built with Dockerfile.gnu.
It is the same with atmos_ocean_fluxes_calc.o.

I was running the 2021.03 version but the run returned with a fatal error as it cannot save some variables in the diag table.

For example, there are "frac_soil" and "area_soil" in diag table,
"land", "frac_soil", "frac_soil", "land_static", "all", .false., "none", 2

but there are no such variables in the code
( that is, grep -r "frac_soil" . returns nothing)

when I remove these terms in diag table, it returned with another error
FATAL from PE 548: diag_manager_mod::init_field_cell_measure: AREA field not in diag_table for field soil/soil_liq

I wonder if the diag table is outdated.

Thanks,

About the output of GFDL-AM4 (eg.
19790101.atmos_daily.tile1.nc.0000
19790101.atmos_daily.tile1.nc.0001
19790101.atmos_daily.tile1.nc.0002
19790101.atmos_daily.tile1.nc.0003
19790101.atmos_daily.tile2.nc.0001...... How to merge and did not find latitude and longitude information. Is it using FRE-NCtool? How to operate exactly? Thank you very much for your help

I download from "ftp://nomads.gfdl.noaa.gov/users/Ming.Zhao/AM4Documentation/GFDL-AM4.0/inputData/AM4_run.tar.gz" AM4_run met below but in running Error message, I hope you can help me

FATAL from PE 0: The checksum in the file:INPUT/fv_core.res.tile1.nc and variable:u does not match the checksum calculated from the data. file:1C000927D95F8925 from data:EAFE4EF3FAAE9405

Can you provide the user guide or description of the namelist parameter? Thank you

&coupler_nml
months = 6,
days = 0,
current_date = 2015,9,1,0,0,0,
calendar = 'julian'
dt_atmos = 1800,
dt_cpld = 7200,
use_lag_fluxes = .true.
concurrent = .false.
do_ocean = .false.
ocean_npes = 0
atmos_npes = 24
atmos_nthreads = 1
radiation_nthreads=1
do_concurrent_radiation=.true.
use_hyper_thread = .false.

NOTE from PE 0: radiation_driver_mod:: concurrent radiation restart active

FATAL from PE 0: The domain associated with the file: does not have an io_domain.

Dear all,
I tried to compile the AM4 model with the intel ifort, mpiifort based on the configuration of the intel.mk template.
I encountered an error in the file AM4/src/FMS/mpp/include/mpp_type_mpi.h, which is related to the definition of the rank
The error message is as follows:
make[1]: warning: Clock skew detected. Your build may be incomplete.
make[1]: Leaving directory /disk01/liuy/models/AM4/exec/coupler' mpiifort coupler/libcoupler.a atmos_dyn/libatmos_dyn.a ice_sis/libice_sis.a atmos_cubed_sphere/libatmos_cubed_sphere.a atmos_phys/libatmos_phys.a mom6/libmom6.a land_lad2/libland_lad2.a fms/libfms.a -L/disk01/liuy/apps/netcdf/lib -lnetcdf -lnetcdff -L/opt/intel/impi/5.0.2.044/intel64/lib -lmpi -L/disk01/liuy/apps/hdf5/lib -lhdf5 -lhdf5_fortran -lhdf5_hl -lhdf5hl_fortran -L/opt/intel/mkl/lib/intel64 -lmkl_blas95_lp64 -lmkl_lapack95_lp64 -lmkl_intel_lp64 -lmkl_core -lmkl_sequential -o am4_xanadu_2021.02.x ld: warning: libimf.so, needed by /disk01/liuy/apps/netcdf/lib/libnetcdf.so, not found (try using -rpath or -rpath-link) ld: warning: libsvml.so, needed by /disk01/liuy/apps/netcdf/lib/libnetcdf.so, not found (try using -rpath or -rpath-link) ld: warning: libirng.so, needed by /disk01/liuy/apps/netcdf/lib/libnetcdf.so, not found (try using -rpath or -rpath-link) ld: warning: libintlc.so.5, needed by /disk01/liuy/apps/netcdf/lib/libnetcdf.so, not found (try using -rpath or -rpath-link) ld: warning: libifport.so.5, needed by /disk01/liuy/apps/netcdf/lib/libnetcdff.so, not found (try using -rpath or -rpath-link) fms/libfms.a(mpp.o): In function mpp_type_create_int4':
/disk01/liuy/models/AM4/src/FMS/mpp/include/mpp_type_mpi.h:46: undefined reference to rank_' /disk01/liuy/models/AM4/src/FMS/mpp/include/mpp_type_mpi.h:60: undefined reference to rank_'
/disk01/liuy/models/AM4/src/FMS/mpp/include/mpp_type_mpi.h:75: undefined reference to rank_' /disk01/liuy/models/AM4/src/FMS/mpp/include/mpp_type_mpi.h:76: undefined reference to rank_'
/disk01/liuy/models/AM4/src/FMS/mpp/include/mpp_type_mpi.h:77: undefined reference to rank_' fms/libfms.a(mpp.o):/disk01/liuy/models/AM4/src/FMS/mpp/include/mpp_type_mpi.h:81: more undefined references to rank_' follow
ld: am4_xanadu_2021.02.x: hidden symbol `__intel_cpu_features_init_x' in /opt/intel/composer_xe_2015.1.133/compiler/lib/intel64/libirc.a(cpu_feature_disp.o) is referenced by DSO
ld: final link failed: Nonrepresentable section on output
make: *** [am4_xanadu_2021.02.x] Error 1

How do can fix this issue.
Thanks so much

I would like to ask if AM4 has a setting to increase the heat source (eg.forceing_nml) just like( https://github.com/lukelbd/gfdl-fms/blob/master/exp/input.nml ) gfdl-fms
&forcing_nml
no_forcing = .false., ! turn off forcing altogether (same as setting each 'k' to zero)
locked_heating = .false., ! use input file for diabatic heating
teq_mode = 'hs', ! can be hs, pk, or pkmod (i.e. my custom stuff)
damp_mode = 'hs', ! can be hs, pk, or pkmod (i.e. my custom stuff)
conserve_energy = .true., ! employ frictional warming
surf_schneider = .true., ! conserve mean surface temp when delh is changed?
strat_vtx = .true., ! (pk/pkmod) turn on and off the polar vortex
strat_sponge = .true., ! (pk/pkmod) turn on sponge
strat_damp = 'constant', ! (pkmod) 'linear': transition from kstrat at tropopause to kmeso at TOA linearly; 'constant': use kstrat everywhere in
t_mean = 300.0, ! mean surface temp
t_zero = 315.0, ! temp at equator surface
t_strat = 200.0, ! using pk stratosphere, will offset US standard temperature by minus 216.65
delh = 60.0, ! equator-pole temp contrast
delv = 10.0, ! lapse rate scaling, probably should be left alone
eps = 0.0, ! adds seasons; delh=delh+eps in north hemisphere, delh=delh-eps in south hemisphere
exp_h = 0, ! positive number shifts baroclinic zone poleward from hs, negative equatorward
exp_b = 4, ! cosine latitude exponent for boundary layer damping equation; Shneider uses 8, Held-Suarez 4
sigma_b = 0.7, ! boundary level; below which friction and elevated thermal damping start
z_ozone = 20.0, ! (pk) level at which stratospheric warming begins for U.S. standard atmosphere (km)
z_kdepth = 5.0, ! (pkmod) depth of transition region for troposphere-to-stratosphere damping
lat_ref = 0.0, ! (pkmod) latitude at which we use exact US standard atmosphere; everywhere else, it is warped
p_sponge = 0.5, ! (pk/pkmod) level above which sponge on winds is applied (hPa)
p_logeval = 200.0, ! (pkmod) pressure at which we evaluate the log(p/p0) when inverting Held-Suarez
vtx_edge = 50.0, ! (pk/pkmod) edge of the vortex
vtx_width = 20.0, ! (pk/pkmod) width of the vortext-US standard atmosphere transition region
vtx_gamma = 2.0, ! (pk/pkmod) lapse rate in polar vortex region (K/km), paper uses 2 and 4
ksponge = -0.5, ! (pk/pkmod) sponge damping rate
kfric = -1.0, ! friction damping rate
kbl = -4.0, ! maximum boundary layer thermal damping, realized at equatorial surface
ktrop = -40.0, ! troposphere damping (if hs, also applies in stratosphere)
kstrat = -40.0, ! (pk/pkmod) stratosphere damping
kmeso = -4.0, ! (pk/pkmod) damping rate at upper boundary, Holton-Mass style
trsink = -4.0, ! damping timescale for tracer, if tracers are registered
trflux = 1.e-5, ! surface flux for tracer, if tracers are registered
q0_tropical = 0.0, ! amplitude for Butler upper-troposphere tropical heating
q0_vortex = 0.0, ! amplitude for Butler polar vortex heating
q0_arctic = 0.0, ! amplitude for Butler Arctic surface heating
q0_global = 0.0, ! globally uniform heating
q0_surface = 0.0, ! average strength of boundary layer heating, that decays linearly from surface to top of boundary layer
q0_realistic = 0.0, ! strength of 'realistically' vertically scaled heating at 1000hPa
q0_lsp = 0.0, ! intensity of lsp forcing
m_lsp = 1.0, ! next lines are LSP heating terms
p0_lsp = 800.0,
pt_lsp = 200.0,
lat0_lsp = 45.0,
slat_lsp = 10.0,