Changeset 3772 for trunk/LMDZ.PLUTO/libf/phypluto/iniwrite.F

Timestamp:

May 21, 2025, 6:37:40 PM (2 months ago)

Author:

afalco

Message:

Pluto: allow to write controle in XIOS output.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/iniwrite.F (modified) (11 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/iniwrite.F

@@ -2 +2 @@
 
       use comsoil_h, only: mlayer, nsoilmx
-      USE comcstfi_mod, only: g, mugaz, omeg, rad, rcp, pi 
+      USE comcstfi_mod, only: g, mugaz, omeg, rad, rcp, pi
       USE vertical_layers_mod, ONLY: ap,bp,aps,bps,pseudoalt
 !      USE logic_mod, ONLY: fxyhypb,ysinus
@@ -10 +10 @@
       USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
       USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev
+      use tabfi_mod, only: tab_cntrl_mod
       IMPLICIT NONE
 
@@ -41 +42 @@
 c   Local:
 c   ------
-      INTEGER length,l
-      parameter (length = 100)
-      REAL tab_cntrl(length) ! run parameters are stored in this array
+    !   INTEGER length,l
+    !   parameter (length = 100)
+    !   REAL tab_cntrl(length) ! run parameters are stored in this array
       INTEGER ierr
       REAl,ALLOCATABLE :: lon_reg_ext(:) ! extended longitudes
@@ -50 +51 @@
       integer :: idim_rlatu,idim_llmp1,idim_llm
       integer :: idim_nsoilmx ! "subsurface_layers" dimension ID #
-      integer, dimension(2) :: id  
+      integer, dimension(2) :: id
 c-----------------------------------------------------------------------
 
@@ -61 +62 @@
       ENDIF
 
-      DO l=1,length
-         tab_cntrl(l)=0.
-      ENDDO
-      tab_cntrl(1)  = real(nbp_lon)
-      tab_cntrl(2)  = real(nbp_lat-1)
-      tab_cntrl(3)  = real(nbp_lev)
-      tab_cntrl(4)  = real(idayref)
-      tab_cntrl(5)  = rad
-      tab_cntrl(6)  = omeg
-      tab_cntrl(7)  = g
-      tab_cntrl(8)  = mugaz
-      tab_cntrl(9)  = rcp
-      tab_cntrl(10) = daysec
-      tab_cntrl(11) = dtphys
+    !   DO l=1,length
+    !      tab_cntrl(l)=0.
+    !   ENDDO
+    !   tab_cntrl(1)  = real(nbp_lon)
+    !   tab_cntrl(2)  = real(nbp_lat-1)
+    !   tab_cntrl(3)  = real(nbp_lev)
+    !   tab_cntrl(4)  = real(idayref)
+    !   tab_cntrl(5)  = rad
+    !   tab_cntrl(6)  = omeg
+    !   tab_cntrl(7)  = g
+    !   tab_cntrl(8)  = mugaz
+    !   tab_cntrl(9)  = rcp
+    !   tab_cntrl(10) = daysec
+    !   tab_cntrl(11) = dtphys
 !      tab_cntrl(12) = etot0
 !      tab_cntrl(13) = ptot0
@@ -82 +83 @@
 c
 c    ..........    P.Le Van  ( ajout le 8/04/96 )    .........
-c         .....        parametres  pour le zoom          ......   
+c         .....        parametres  pour le zoom          ......
 !      tab_cntrl(17)  = clon
 !      tab_cntrl(18)  = clat
@@ -105 +106 @@
 
 ! Define dimensions
-    
-      ierr = NF_REDEF (nid)
-
-      ierr = NF_DEF_DIM (nid, "index", length, idim_index)
+
+      ierr = NF_REDEF (nid)
+
+      ierr = NF_DEF_DIM (nid, "index", SIZE(tab_cntrl_mod), idim_index)
 !      ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
       ierr = NF_DEF_DIM (nid, "latitude", nbp_lat, idim_rlatu)
@@ -126 +127 @@
       ierr = NF_REDEF (nid)
 #ifdef NC_DOUBLE
-      ierr = NF_DEF_VAR (nid, "controle", NF_DOUBLE, 1, 
+      ierr = NF_DEF_VAR (nid, "controle", NF_DOUBLE, 1,
      .       idim_index,nvarid)
 #else
-      ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1, 
+      ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1,
      .       idim_index,nvarid)
 #endif
@@ -136 +137 @@
       ierr = NF_ENDDEF(nid)
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
-#else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl_mod)
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl_mod)
 #endif
 
@@ -177 +178 @@
 c
 c --------------------------
-      
+
       lon_reg_ext(1:nbp_lon)=lon_reg(1:nbp_lon)
       IF (nbp_lon*nbp_lat/=1) THEN
@@ -204 +205 @@
       ierr = NF_REDEF (nid)
 #ifdef NC_DOUBLE
-      ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1, 
+      ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1,
      .       idim_llm,nvarid)
 #else
-      ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1, 
+      ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1,
      .       idim_llm,nvarid)
 #endif