Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/calcul_STDlev.h

Timestamp:

Jul 28, 2025, 7:23:15 PM (6 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/calcul_STDlev.h (modified) (3 diffs)

LMDZ6/branches/contrails/libf/phylmd/calcul_STDlev.h

@@ -2 +2 @@
 ! $Header$
 !
+!-------------------------------------------------------------------
+! Interpolation sur des niveaux de pression.
+!
+! Enchainement des operations en prenant pour exemple la temperature
+!
+! 1. CALL ini_undefSTD :
+!  On met les champs cumules à 0 si 
+!  if mod(itap,nint(freq_outnmc(n)/phys_tstep))==1 
+!          tnondef(i, k, n) = 0.
+!          tsumstd(i, k, n) = 0.
+!
+! 2. Ensuite : 
+!   on interpole les champs sur les niveaux STD de pression 
+!   a chaque pas de temps de la physique avec plevel_new
+!   t_seri -> tlevSTD
+!
+! 3. CALL undefSTD : 
+!   On ajoute les champs à la moyenne tous les frec_calnmc
+!      pas de temps
+!   IF (mod(itap,nint(freq_calnmc(n)/phys_tstep))==0) THEN
+!         IF (tlevstd(i,k)==missing_val) THEN
+!           tnondef(i, k, n) = tnondef(i, k, n) + 1.
+!         ELSE IF (tlevstd(i,k)/=missing_val) THEN
+!           tsumstd(i, k, n) = tsumstd(i, k, n) + tlevstd(i, k)
+!
+! 4. CALL moy_undefSTD
+!    On normalise à freq_outnmc, avant de reinitialiser plus haut dans ini_undef
+!    IF mod(itap,nint(freq_outnmc(n)/phys_tstep))==0
+!         tsumstd(i, k, n) = tsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k,n))
+!
+! 5. Intrpolation des flux (pourquoi pas avec 2 ?)
+!
+! 6. Choix entre instantanné et moyenne en temps
+!     twriteSTD(:,:,1)=tsumSTD(:,:,1)
+!     twriteSTD(:,:,3)=tlevSTD(:,:)
+! 
+!  Premiere série de modifications (2025/05/07, FH) :
+!   - Fait en sorte que freq_outnmc=ecrit_files
+!   - Mis tous les USE avec ONLY dans *undefSTD.f90
+!   - Utilie ecrit_files comme frequence de sortie freq_outnmc
+!   - Nettoye la lecture des freq_outnmc et freq_calnmc dans conf_phys
+!   - Passe les dimensionnement de certaines variables de 3 à 10
+!     nombre de fichiers par défaut. En gros nout=nfiles
+!  
+!  Proposition de modifications :
+!   - Ne plus calculer l'interpolation à tous les pas de temps
+!      mais seulement à la fréquence freq_calcnmc
+!   - Mettre le call ini_undefSTD juste avant undef_STD et moy_undefSTD
+!   - Enlever ecrit_files(7) = freq_outNMC(1) dans phys_output_mod
+!   - Enlever les vieux controles de fichiers dans conf_phys
+!   - enlever les dimensionnement à 10 en dur
+!   - Enlever toute chaine de caractère NMC dans le modèle.
+!
+!-------------------------------------------------------------------
+
 !IM on initialise les variables 
 !
 !       missing_val=nf90_fill_real
 !
+
+! On peut à terme utiliser directement ecrit_files (le nom est moins
+! bien choisi ...
+
+       !print*,'STDSTD ecrit_files ',ecrit_files(1:4)
+       freq_outnmc(:)=ecrit_files(:)
+       freq_calnmc(:)=pdtphys
+
        IF (.not. ok_all_xml) then
-        CALL ini_undefSTD(itap,itapm1)
+        CALL ini_undefSTD(itap)
        ENDIF
-!
-!IM on interpole les champs sur les niveaux STD de pression 
-!IM a chaque pas de temps de la physique
 !
 !-------------------------------------------------------c
@@ -134 +194 @@
        IF (.not. ok_all_xml) then
         CALL undefSTD(itap, read_climoz)
-        CALL moy_undefSTD(itap,itapm1)
+        CALL moy_undefSTD(itap)
        ENDIF
 !
@@ -160 +220 @@
            lwup,LWup200)
 !
-      twriteSTD(:,:,1)=tsumSTD(:,:,1)
-      qwriteSTD(:,:,1)=qsumSTD(:,:,1)
-      rhwriteSTD(:,:,1)=rhsumSTD(:,:,1)
-      phiwriteSTD(:,:,1)=phisumSTD(:,:,1)
-      uwriteSTD(:,:,1)=usumSTD(:,:,1)
-      vwriteSTD(:,:,1)=vsumSTD(:,:,1)
-      wwriteSTD(:,:,1)=wsumSTD(:,:,1)
-
-      twriteSTD(:,:,2)=tsumSTD(:,:,2)
-      qwriteSTD(:,:,2)=qsumSTD(:,:,2)
-      rhwriteSTD(:,:,2)=rhsumSTD(:,:,2)
-      phiwriteSTD(:,:,2)=phisumSTD(:,:,2)
-      uwriteSTD(:,:,2)=usumSTD(:,:,2)
-      vwriteSTD(:,:,2)=vsumSTD(:,:,2)
-      wwriteSTD(:,:,2)=wsumSTD(:,:,2)
-
-      twriteSTD(:,:,3)=tlevSTD(:,:)
-      qwriteSTD(:,:,3)=qlevSTD(:,:)
-      rhwriteSTD(:,:,3)=rhlevSTD(:,:)
-      phiwriteSTD(:,:,3)=philevSTD(:,:)
-      uwriteSTD(:,:,3)=ulevSTD(:,:)
-      vwriteSTD(:,:,3)=vlevSTD(:,:)
-      wwriteSTD(:,:,3)=wlevSTD(:,:)
-
-      twriteSTD(:,:,4)=tlevSTD(:,:)
-      qwriteSTD(:,:,4)=qlevSTD(:,:)
-      rhwriteSTD(:,:,4)=rhlevSTD(:,:)
-      phiwriteSTD(:,:,4)=philevSTD(:,:)
-      uwriteSTD(:,:,4)=ulevSTD(:,:)
-      vwriteSTD(:,:,4)=vlevSTD(:,:)
-      wwriteSTD(:,:,4)=wlevSTD(:,:)
-!
-!IM initialisation 5eme fichier de sortie 
-      twriteSTD(:,:,5)=tlevSTD(:,:)
-      qwriteSTD(:,:,5)=qlevSTD(:,:)
-      rhwriteSTD(:,:,5)=rhlevSTD(:,:)
-      phiwriteSTD(:,:,5)=philevSTD(:,:)
-      uwriteSTD(:,:,5)=ulevSTD(:,:)
-      vwriteSTD(:,:,5)=vlevSTD(:,:)
-      wwriteSTD(:,:,5)=wlevSTD(:,:)
-!
-!IM initialisation 6eme fichier de sortie 
-      twriteSTD(:,:,6)=tlevSTD(:,:)
-      qwriteSTD(:,:,6)=qlevSTD(:,:)
-      rhwriteSTD(:,:,6)=rhlevSTD(:,:)
-      phiwriteSTD(:,:,6)=philevSTD(:,:)
-      uwriteSTD(:,:,6)=ulevSTD(:,:)
-      vwriteSTD(:,:,6)=vlevSTD(:,:)
-      wwriteSTD(:,:,6)=wlevSTD(:,:)
+      do n=1,nout
+         !print*,'STDSTD n nout',n,nout,clef_files(n)
+         if (clef_files(n)) then
+             if (type_ecri(n)=='ave(X)') then
+                 !print*,'STDSTD fichier ave(X), n=',n
+                 twriteSTD(:,:,n)=tsumSTD(:,:,n)
+                 qwriteSTD(:,:,n)=qsumSTD(:,:,n)
+                 rhwriteSTD(:,:,n)=rhsumSTD(:,:,n)
+                 phiwriteSTD(:,:,n)=phisumSTD(:,:,n)
+                 uwriteSTD(:,:,n)=usumSTD(:,:,n)
+                 vwriteSTD(:,:,n)=vsumSTD(:,:,n)
+                 wwriteSTD(:,:,n)=wsumSTD(:,:,n)
+             else
+                 !print*,'STDSTD fichier autre, n=',n,type_ecri(n)
+                 twriteSTD(:,:,n)=tlevSTD(:,:)
+                 qwriteSTD(:,:,n)=qlevSTD(:,:)
+                 rhwriteSTD(:,:,n)=rhlevSTD(:,:)
+                 phiwriteSTD(:,:,n)=philevSTD(:,:)
+                 uwriteSTD(:,:,n)=ulevSTD(:,:)
+                 vwriteSTD(:,:,n)=vlevSTD(:,:)
+                 wwriteSTD(:,:,n)=wlevSTD(:,:)
+             endif
+         endif
+      enddo
+
 !IM for NMC files
       DO n=1, nlevSTD3