Changeset 3080 for trunk/LMDZ.MARS/libf

Timestamp:

Oct 10, 2023, 4:02:50 PM (14 months ago)

Author:

emillour

Message:

Mars PCM: follow-up of r3078 which broke the 1D model.
Also added initialisation of non-oro GW tendencies stored in
startfi.nc when in 1D.
EM+JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -25 +25 @@
 use regular_lonlat_mod,       only: init_regular_lonlat
 use mod_interface_dyn_phys,   only: init_interface_dyn_phys
-use geometry_mod,             only: init_geometry
+use geometry_mod,             only: init_geometry, init_geometry_cell_area_for_outputs
 use dimphy,                   only: init_dimphy
 use comgeomfi_h,              only: sinlat, ini_fillgeom
@@ -38 +38 @@
 use phys_state_var_init_mod,  only: phys_state_var_init
 use turb_mod,                 only: q2
+use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd
 ! Mostly for XIOS outputs:
 use mod_const_mpi,            only: COMM_LMDZ
@@ -429 +430 @@
 call init_regular_lonlat(1,1,longitude,latitude,(/0.,0./),(/0.,0./))
 call init_geometry(1,longitude,latitude,(/0.,0.,0.,0./),(/0.,0.,0.,0./),cell_area)
+call init_geometry_cell_area_for_outputs(1,cell_area)
 ! Ehouarn: init_vertial_layers called later (because disvert not called yet)
 !      call init_vertical_layers(nlayer,preff,scaleheight,
@@ -503 +505 @@
 zgam(1) = 0.
 zthe(1) = 0.
+
+! For the non-orographic gravity waves scheme
+du_nonoro_gwd(1,:)=0
+dv_nonoro_gwd(1,:)=0
 
 ! For the slope wind scheme