Changeset 2907 for trunk/LMDZ.MARS/libf/phymars/iostart.F90

Timestamp:

Feb 28, 2023, 2:39:58 PM (21 months ago)

Author:

romain.vande

Message:

Mars PCM:
Further implementation of subslope parametrisation.
Carefull ! This is a devolpment revision and it still need improvements and tests.
However, this commit should not change anything for nslope=1. Only nslope=1 is possible for now!

Sloped variables in the startfi are outputed along the nslope dimension.

Possibility to output variables on a specific subslope in diagfi using VAR_slopeXX in the diagfi.def,
with VAR the variable name (ex: tsurf, co2ice...) and XX the slope number (ex: 04, 07...).
Without any specific mention to slope, variable named VAR in the diagfi.def will correspond to the variable in the flat slope, this can change in the future.

Other code changes for nslope.gt.1 (sometimes the grid mesh average is used instead of the value of the subslope):

Changes for albedo :
albedo_mesh_avg is used in callradite

Changes for emis :
emis_mesh_avg is used in callradite

Changes for zdqsdif:
zdqsdif_mesh_avg is used for lifting in turbulent-resolving mode

Changes for qsurf:
qsurf_mesh_avg is used for callradite, rocketduststorm, topmons.
qsurf of the major slope is used for calchim and geticecover (tituscap)

Changes for tsurf is the next commit.

File:

• trunk/LMDZ.MARS/libf/phymars/iostart.F90 (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/iostart.F90

@@ -692 +692 @@
             ierr=NF90_REDEF(nid_restart)
 #ifdef NC_DOUBLE
+            if(field_name.eq. "tauscaling" .or. field_name.eq. "totcloudfrac" .or. field_name.eq. "wstar" ) then
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
                               (/idim2,idim7/),nvarid)
-#else
+            else
+            ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
+                              (/idim2,idim8,idim7/),nvarid)
+            endif
+#else
+            if(field_name.eq. "tauscaling" .or. field_name.eq. "totcloudfrac" .or. field_name.eq. "wstar" ) then
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
                               (/idim2,idim7/),nvarid)
+            else
+            ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
+                              (/idim2,idim8,idim7/),nvarid)
+            endif
 #endif
             if (ierr.ne.NF90_NOERR) then