Changeset 3623 for trunk/LMDZ.PLUTO/libf

Timestamp:

Feb 12, 2025, 7:09:35 PM (22 hours ago)

Author:

afalco

Message:

Pluto: better message for change in "temperature" written by diagfi (replaced the deprecated "temp" name).
zrecast also reads "temperature" if "temp" fails.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/writediagfi.F (modified) (21 diffs)

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

@@ -1 +1 @@
       subroutine writediagfi(ngrid,nom,titre,unite,dim,px)
 
-!  Ecriture de variables diagnostiques au choix dans la physique 
+!  Ecriture de variables diagnostiques au choix dans la physique
 !  dans un fichier NetCDF nomme  'diagfi'. Ces variables peuvent etre
 !  3d (ex : temperature), 2d (ex : temperature de surface), ou
@@ -7 +7 @@
 !  solaire)
 !  (ou encore 1d, dans le cas de testphys1d, pour sortir une colonne)
-!  La periode d'ecriture est donnee par 
+!  La periode d'ecriture est donnee par
 !  "ecritphy " regle dans le fichier de controle de run :  run.def
 !
@@ -16 +16 @@
 ! WARNING : les variables dynamique (u,v,t,q,ps)
 !  sauvees par writediagfi avec une
-! date donnee sont legerement differentes que dans le fichier histoire car 
+! date donnee sont legerement differentes que dans le fichier histoire car
 ! on ne leur a pas encore ajoute de la dissipation et de la physique !!!
 ! IL est  RECOMMANDE d'ajouter les tendance physique a ces variables
 ! avant l'ecriture dans diagfi (cf. physiq.F)
-!  
+!
 ! Modifs: Aug.2010 Ehouarn: enforce outputs to be real*4
 !         Oct 2011 Francois: enable having a 'diagfi.def' file to select
@@ -30 +30 @@
 !                (ngrid = 2+(jjm-1)*iim - 1/jjm)
 !                 (= nlon ou klon dans la physique terrestre)
-!      
+!
 !      unit : unite logique du fichier de sortie (toujours la meme)
 !      nom  : nom de la variable a sortir (chaine de caracteres)
@@ -99 +99 @@
       logical,save :: firstcall=.true.
 !$OMP THREADPRIVATE(firstcall)  !diagfi_def,n_nom_def,nom_def read in diagfi.def
-      
+
 #ifdef CPP_PARA
 ! Added to work in parallel mode
@@ -145 +145 @@
               read(99,fmt='(a)',end=88) nom_def(n)
               write(*,*) 'Output in diagfi: ', trim(nom_def(n))
-            end do 
+              if (nom_def(n).eq.'temp') then
+                 write(*,*) 'WARNING: "temp" is outdated.'
+                 write(*,*) 'Do you really have a "temp" variable?'
+                 write(*,*) 'Otherwise use "temperature" instead.'
+              end if
+            end do
  88         continue
             if (n.ge.n_nom_def_max) then
@@ -151 +156 @@
                call abort_physic("writediagfi",
      &             "n_nom_def_max too small",1)
-            end if 
+            end if
             n_nom_def=n-1
             close(99)
@@ -185 +190 @@
            call abort_physic("writediagfi","firstnom too short",1)
          endif
-         
+
          zitau = -1 ! initialize zitau
 
@@ -242 +247 @@
           enddo
          ENDIF
-         
+
          ! write "header" of file (longitudes, latitudes, geopotential, ...)
          IF (klon_glo>1) THEN ! general 3D case
@@ -299 +304 @@
               write(*,*) "***** PUT_VAR matter in writediagfi_nc"
               write(*,*) "***** with time"
-              write(*,*) 'ierr=', ierr,": ",NF_STRERROR(ierr)  
+              write(*,*) 'ierr=', ierr,": ",NF_STRERROR(ierr)
 c             call abort
            endif
@@ -372 +377 @@
      &             trim(nom)//" already exists",1)
              endif
-           endif 
+           endif
 
            corner(1)=1
@@ -407 +412 @@
               call abort_physic("writediagfi",
      &             "failed writing "//trim(nom),1)
-           endif 
+           endif
 
           endif !of if (is_master)
@@ -491 +496 @@
 
 !#ifdef NC_DOUBLE
-!           ierr = NF_PUT_VARA_DOUBLE (nid,varid,corner,edges,dx2) 
-!#else         
+!           ierr = NF_PUT_VARA_DOUBLE (nid,varid,corner,edges,dx2)
+!#else
            IF (klon_glo>1) THEN ! General case
              ierr= NF_PUT_VARA_REAL(nid,varid,corner,edges,dx2)
@@ -498 +503 @@
              ierr= NF_PUT_VARA_REAL(nid,varid,corner,edges,dx2_1d)
            ENDIF
-!#endif     
+!#endif
 
            if (ierr.ne.NF_NOERR) then
@@ -506 +511 @@
               call abort_physic("writediagfi",
      &             "failed writing "//trim(nom),1)
-           endif 
+           endif
 
           endif !of if (is_master)
@@ -525 +530 @@
             dx1(l)=px(1,l)
           enddo
-          
+
           ierr= NF_INQ_VARID(nid,nom,varid)
            if (ierr /= NF_NOERR) then
@@ -538 +543 @@
 
               call def_var(nid,nom,titre,unite,2,id,varid,ierr)
-              
+
            else
              if (ntime==0) then
@@ -548 +553 @@
              endif
            endif
-           
+
            corner(1)=1
            corner(2)=ntime
-           
+
            edges(1)=nbp_lev
            edges(2)=1
@@ -566 +571 @@
               call abort_physic("writediagfi",
      &             "failed writing "//trim(nom),1)
-           endif 
+           endif
 
 !Case of a 0D variable (ie: a time-dependent scalar)
@@ -603 +608 @@
 
 !#ifdef NC_DOUBLE
-!           ierr = NF_PUT_VARA_DOUBLE (nid,varid,corner,edges,dx0)  
+!           ierr = NF_PUT_VARA_DOUBLE (nid,varid,corner,edges,dx0)
 !#else
            ierr= NF_PUT_VARA_REAL(nid,varid,corner,edges,dx0)
@@ -613 +618 @@
               call abort_physic("writediagfi",
      &             "failed writing "//trim(nom),1)
-           endif 
+           endif
 
           endif !of if (is_master)