Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/ini_undefSTD.f90

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/ini_undefSTD.f90 (modified) (3 diffs)

LMDZ6/branches/contrails/libf/phylmd/ini_undefSTD.f90

@@ -3 +3 @@
 
 
-SUBROUTINE ini_undefstd(itap, itapm1)
-  USE clesphys_mod_h
-  USE dimphy
-  USE phys_state_var_mod ! Variables sauvegardees de la physique
-  USE phys_cal_mod, ONLY: mth_len
-  IMPLICIT NONE
+SUBROUTINE ini_undefstd(itap)
+USE clesphys_mod_h, ONLY : freq_outnmc
+USE dimphy, ONLY : klon
+USE phys_state_var_mod, ONLY : tnondef,tsumstd,usumstd,vsumstd,wsumstd,phisumstd,qsumstd,rhsumstd
+USE phys_state_var_mod, ONLY : uvsumstd,vqsumstd,vtsumstd,wqsumstd,vphisumstd,wtsumstd,u2sumstd
+USE phys_state_var_mod, ONLY : v2sumstd,t2sumstd,o3sumstd,o3daysumstd,nlevSTD,nout,phys_tstep
+USE phys_state_var_mod, ONLY : missing_val_nf90
+USE phys_cal_mod, ONLY: mth_len
+USE phys_output_var_mod, ONLY : clef_files
+
+IMPLICIT NONE
 
   ! ====================================================================
@@ -34 +39 @@
   ! variables Input/Output
   ! INTEGER nlevSTD, klevSTD, itap
-  INTEGER itap, itapm1
+  INTEGER itap
   ! PARAMETER(klevSTD=17)
   ! REAL dtime
@@ -54 +59 @@
     ! initialisation variables en debut du mois, de la journee ou des 6h
 
-    IF (n==1 .AND. itap-itapm1==1 .OR. n>1 .AND. mod(itap,nint( &
-        freq_outnmc(n)/phys_tstep))==1) THEN
-      ! print*,'ini_undefSTD n itap',n,itap
+    !IF ( clef_files(n) .and. mod(itap,nint(freq_outnmc(n)/phys_tstep))==1) THEN
+    IF ( mod(itap,nint(freq_outnmc(n)/phys_tstep))==1) THEN
       DO k = 1, nlevstd
         DO i = 1, klon
           tnondef(i, k, n) = 0.
-          tsumstd(i, k, n) = 0.
-          usumstd(i, k, n) = 0.
-          vsumstd(i, k, n) = 0.
-          wsumstd(i, k, n) = 0.
-          phisumstd(i, k, n) = 0.
-          qsumstd(i, k, n) = 0.
-          rhsumstd(i, k, n) = 0.
-          uvsumstd(i, k, n) = 0.
-          vqsumstd(i, k, n) = 0.
-          vtsumstd(i, k, n) = 0.
-          wqsumstd(i, k, n) = 0.
-          vphisumstd(i, k, n) = 0.
-          wtsumstd(i, k, n) = 0.
-          u2sumstd(i, k, n) = 0.
-          v2sumstd(i, k, n) = 0.
-          t2sumstd(i, k, n) = 0.
-          o3sumstd(i, k, n) = 0.
-          o3daysumstd(i, k, n) = 0.
         END DO !i
       END DO !k
+      IF ( clef_files(n) ) THEN
+        ! print*,'ini_undefSTD n itap',n,itap
+        DO k = 1, nlevstd
+          DO i = 1, klon
+            tnondef(i, k, n) = 0.
+            tsumstd(i, k, n) = 0.
+            usumstd(i, k, n) = 0.
+            vsumstd(i, k, n) = 0.
+            wsumstd(i, k, n) = 0.
+            phisumstd(i, k, n) = 0.
+            qsumstd(i, k, n) = 0.
+            rhsumstd(i, k, n) = 0.
+            uvsumstd(i, k, n) = 0.
+            vqsumstd(i, k, n) = 0.
+            vtsumstd(i, k, n) = 0.
+            wqsumstd(i, k, n) = 0.
+            vphisumstd(i, k, n) = 0.
+            wtsumstd(i, k, n) = 0.
+            u2sumstd(i, k, n) = 0.
+            v2sumstd(i, k, n) = 0.
+            t2sumstd(i, k, n) = 0.
+            o3sumstd(i, k, n) = 0.
+            o3daysumstd(i, k, n) = 0.
+          END DO !i
+        END DO !k
+      END IF !
     END IF !
   END DO !n