Changeset 2913 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/ini_archive.F

Timestamp:

Mar 14, 2023, 10:07:33 AM (2 years ago)

Author:

romain.vande

Message:

Mars PCM:
Adapt start2archive.F to the subslope parametrisation.
Small correction for some dimensions of variables.
RV

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/ini_archive.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/ini_archive.F

@@ -1 +1 @@
 c=======================================================================
-      subroutine ini_archive(nid,idayref,phis,ith,tab_cntrl_fi)
+      subroutine ini_archive(nid,idayref,phis,ith,tab_cntrl_fi
+     &                      ,def_slope,subslope_dist)
 c=======================================================================
 c
@@ -40 +41 @@
       USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy
       USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
+      use comslope_mod, ONLY: nslope
       implicit none
 
@@ -66 +68 @@
       REAL              phis(ip1jmp1)
       real ith(ip1jmp1,nsoilmx)
+      real subslope_dist(ip1jmp1,nslope)
+      real def_slope(nslope+1)
       REAL              tab_cntrl_fi(length)
 
@@ -77 +81 @@
       INTEGER idim_tim
       INTEGER idim_nsoilmx ! "subsurface_layers" dimension ID #
+      INTEGER idim_nslope, idim_nslope_p1
       INTEGER nid,nvarid
       real sig_s(llm),s(llm)
@@ -163 +168 @@
       ierr = NF_DEF_DIM (nid,"interlayer", llmp1, idim_llmp1)
       ierr = NF_DEF_DIM (nid,"Time", NF_UNLIMITED, idim_tim)
+      ierr = NF_DEF_DIM (nid,"nslope", nslope, idim_nslope)
+      ierr = NF_DEF_DIM (nid,"nslope_plus_1",nslope+1,idim_nslope_p1)
 
 c
@@ -513 +520 @@
 #endif
 
+c Put subslope dist
+      dims3(1)=idim_rlonv
+      dims3(2)=idim_rlatu
+      dims3(3)=idim_nslope
+      ierr = NF_REDEF (nid)
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "subslope_dist", NF_DOUBLE, 3,
+     .       dims3,nvarid)
+#else
+      ierr = NF_DEF_VAR (nid, "subslope_dist", NF_FLOAT, 3,
+     .       dims3,nvarid)
+#endif
+      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",13,
+     .    "subslope_dist")
+
+      ierr = NF_ENDDEF(nid)
+#ifdef NC_DOUBLE
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,subslope_dist)
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,subslope_dist)
+#endif
+
+c Put def_slope
+
+      ierr = NF_REDEF (nid)
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "def_slope", NF_DOUBLE, 1,
+     .       [idim_nslope_p1],nvarid)
+#else
+      ierr = NF_DEF_VAR (nid, "def_slope", NF_FLOAT, 1,
+     .       [idim_nslope_p1],nvarid)
+#endif
+      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",7,"def_slope")
+      ierr = NF_ENDDEF(nid)
+#ifdef NC_DOUBLE
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,def_slope)
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,def_slope)
+#endif
+
       PRINT*,'iim,jjm,llm,idayref',iim,jjm,llm,idayref
       PRINT*,'rad,omeg,g,mugaz,kappa',