Changeset 3521 for trunk

Timestamp:

Nov 15, 2024, 6:57:21 PM (6 days ago)

Author:

afalco

Message:

Pluto: start2archive add thermal inertia as param.
AF

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F (modified) (25 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F

@@ -9 +9 @@
 c
 c   Objet:   Passage des  fichiers netcdf d'etat initial "start" et
-c   -----    "startfi" a un fichier netcdf unique "start_archive" 
+c   -----    "startfi" a un fichier netcdf unique "start_archive"
 c
 c  "start_archive" est une banque d'etats initiaux:
 c  On peut stocker plusieurs etats initiaux dans un meme fichier "start_archive"
 c    (Veiller dans ce cas avoir un day_ini different pour chacun des start)
-c 
+c
 c
 c
@@ -21 +21 @@
       use infotrac, only: infotrac_init, nqtot, tname
       USE comsoil_h
-      
+
 !      use slab_ice_h, only: noceanmx
 c      USE ocean_slab_mod, ONLY: nslay
@@ -42 +42 @@
 
       include "dimensions.h"
-      integer, parameter :: ngridmx = (2+(jjm-1)*iim - 1/jjm) 
+      integer, parameter :: ngridmx = (2+(jjm-1)*iim - 1/jjm)
       include "paramet.h"
       include "comdissip.h"
@@ -55 +55 @@
 c -----------------------------
       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm) ! vents covariants
-      REAL teta(ip1jmp1,llm)                    ! temperature potentielle 
+      REAL teta(ip1jmp1,llm)                    ! temperature potentielle
       REAL,ALLOCATABLE :: q(:,:,:)   ! champs advectes
       REAL pks(ip1jmp1)                      ! exner (f pour filtre)
@@ -64 +64 @@
       REAL ps(ip1jmp1)                       ! pression au sol
       REAL p3d(iip1, jjp1, llm+1)            ! pression aux interfaces
-      
+
 c Variable Physiques (grille physique)
 c ------------------------------------
@@ -104 +104 @@
 
       CHARACTER*2 str2
-      CHARACTER*80 fichier 
+      CHARACTER*80 fichier
       data  fichier /'startfi'/
 
@@ -115 +115 @@
 c Netcdf
 c-------
-      integer varid,dimid,timelen 
+      integer varid,dimid,timelen
       INTEGER nid,nid1
 
 c-----------------------------------------------------------------------
-c   Initialisations 
+c   Initialisations
 c-----------------------------------------------------------------------
 
@@ -152 +152 @@
         CALL ABORT
        ENDIF
-                                                
+
       ierr = NF_INQ_VARID (nid1, "controle", varid)
       IF (ierr .NE. NF_NOERR) THEN
@@ -188 +188 @@
       ENDIF
       ierr = NF_CLOSE(nid1)
-      
+
       ! allocate arrays of nsoilmx size
       allocate(tsoil(ngridmx,nsoilmx))
@@ -195 +195 @@
 
 c-----------------------------------------------------------------------
-c   Initialisations 
+c   Initialisations
 c-----------------------------------------------------------------------
 
@@ -216 +216 @@
 ! Allocate saved arrays (as in firstcall of physiq)
       call phys_state_var_init(nqtot)
-      
+
 ! Initialize tracer names, indexes and properties
       CALL initracer(ngridmx,nqtot)
@@ -222 +222 @@
       CALL phyetat0(.true.,ngridmx,llm,fichnom,0,Lmodif,nsoilmx,nqtot,
      .      day_ini_fi,timefi,
-     .      tsurf,tsoil,emis,q2,qsurf) 
+     .      tsurf,tsoil,emis,q2,qsurf,inertiedat)
 !       change FF 05/2011
 c     .       cloudfrac,totalcloudfrac,hice,
@@ -238 +238 @@
         CALL ABORT
        ENDIF
-                                                
+
       ierr = NF_INQ_VARID (nid1, "controle", varid)
       IF (ierr .NE. NF_NOERR) THEN
@@ -260 +260 @@
 c Controle de la synchro
 c-----------------------------------------------------------------------
-!mars a voir      if ((day_ini_fi.ne.day_ini).or.(abs(timefi-timedyn).gt.1.e-10)) 
-      if ((day_ini_fi.ne.day_ini)) 
+!mars a voir      if ((day_ini_fi.ne.day_ini).or.(abs(timefi-timedyn).gt.1.e-10))
+      if ((day_ini_fi.ne.day_ini))
      &  stop ' Probleme de Synchro entre start et startfi !!!'
 
@@ -284 +284 @@
 
 c-----------------------------------------------------------------------
-c Vent   (depend de la resolution horizontale) 
+c Vent   (depend de la resolution horizontale)
 c-----------------------------------------------------------------------
 c
@@ -292 +292 @@
 c-----------------------------------------------------------------------
 
-      call covnat(llm,ucov, vcov, un, vn) 
-      call wind_scal(un,vn,us,vs) 
+      call covnat(llm,ucov, vcov, un, vn)
+      call wind_scal(un,vn,us,vs)
 
 c-----------------------------------------------------------------------
@@ -310 +310 @@
 
 c-----------------------------------------------------------------------
-c Variable physique 
+c Variable physique
 c-----------------------------------------------------------------------
 c
@@ -353 +353 @@
          DO i=1,iim
            ptotal=ptotal+aire(i+(iim+1)*(j-1))*ps(i+(iim+1)*(j-1))/g
-!           n2icetotal = n2icetotal + 
+!           n2icetotal = n2icetotal +
 !     &            n2iceS(i+(iim+1)*(j-1))*aire(i+(iim+1)*(j-1))
          ENDDO
@@ -372 +372 @@
 
 c-----------------------------------------------------------------------
-c Ouverture de "start_archive" 
+c Ouverture de "start_archive"
 c-----------------------------------------------------------------------
 
       ierr = NF_OPEN ('start_archive.nc', NF_WRITE,nid)
- 
+
 c-----------------------------------------------------------------------
 c  si "start_archive" n'existe pas:
@@ -383 +383 @@
 c-----------------------------------------------------------------------
 c ini_archive:
-c On met dans l'entete le tab_cntrl dynamique (1 a 16) 
+c On met dans l'entete le tab_cntrl dynamique (1 a 16)
 c  On y ajoute les valeurs du tab_cntrl_fi (a partir de 51)
 c  En plus les deux valeurs ptotal et n2icetotal (99 et 100)
@@ -426 +426 @@
 c    On passe donc une des couches de q2 a part
 c    comme une variable 2D (la couche au sol: "q2surf")
-c    Les lmm autres couches sont nommees "q2atm" (3D) 
+c    Les lmm autres couches sont nommees "q2atm" (3D)
 c-----------------------------------------------------------------------
 
@@ -463 +463 @@
 c-----------------------------------------------------------------------
 c "tsoil" Temperature au sol definie dans 10 couches dans le sol
-c   Les 10 couches sont lues comme 10 champs 
+c   Les 10 couches sont lues comme 10 champs
 c  nommees Tg[1,10]
 
@@ -523 +523 @@
 c     &            ,'tsea_ice','',2,tsea_iceS)
 c        endif !ok_slab_ocean
-      
+
       ENDIF ! of IF(ierr.EQ.0)
 
@@ -540 +540 @@
 
 c-----------------------------------------------------------------------
-c Fin 
+c Fin
 c-----------------------------------------------------------------------
       ierr=NF_CLOSE(nid)
@@ -546 +546 @@
       write(*,*) "start2archive: All is well that ends well."
 
-      end 
+      end