Changeset 1992 for LMDZ5/trunk/libf/phylmd/undefSTD.F90

Timestamp:

Mar 5, 2014, 2:19:12 PM (11 years ago)

Author:

lguez

Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:

• LMDZ5/trunk/libf/phylmd/undefSTD.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/undefSTD.F) (1 diff)

LMDZ5/trunk/libf/phylmd/undefSTD.F90

@@ -1 @@
-!
+
 ! $Id$
-!
-      SUBROUTINE undefSTD(itap, read_climoz)
-      USE netcdf
-      USE dimphy
-      USE phys_state_var_mod ! Variables sauvegardees de la physique
-      IMPLICIT none
-      include "clesphys.h"
-c
-c====================================================================
-c
-c I. Musat : 09.2004
-c
-c Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef)) 
-c          ou la variable tlevSTD est bien definie (.NE.missing_val), 
-c et 
-c        * de la somme de tlevSTD => tsumSTD
-c
-c nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
-c        ! de la physique
-c nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
-c        ! de la physique
-c nout=3 !var. mensuelle "NMC" moyennee toutes les ecrit_hf
-c
-c
-c NB: mettre "inst(X)" dans le write_hist*NMC.h !
-c====================================================================
-c
-cym#include "dimensions.h"
-cym      integer jjmp1
-cym      parameter (jjmp1=jjm+1-1/jjm)
-cym#include "dimphy.h"
-c variables Input
-c
-c     INTEGER nlevSTD, klevSTD, itap
-c     PARAMETER(klevSTD=17)
-      INTEGER itap
-c     REAL dtime
-c
-c variables locales
-c     INTEGER i, k, nout, n
-c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
-      INTEGER i, k, n
-c     REAL freq_calNMC(nout)
-      INTEGER read_climoz
-c
-c variables Output
-c     REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
-c     LOGICAL oknondef(klon,klevSTD,nout)
-c     REAL tnondef(klon,klevSTD,nout)
-c
-      REAL missing_val
-c
-      missing_val=nf90_fill_real
-c
-      DO n=1, nout
-c
-c
-c calcul variables tous les freq_calNMC(n)/dtime pas de temps 
-c de la physique
-c
-       IF(MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0) THEN
-        DO k=1, nlevSTD
-         DO i=1, klon
-          IF(tlevSTD(i,k).EQ.missing_val) THEN
-c          IF(oknondef(i,k,n)) THEN          
-            tnondef(i,k,n)=tnondef(i,k,n)+1.
-c          ENDIF !oknondef(i,k)
-c
-          ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
-           tsumSTD(i,k,n)=tsumSTD(i,k,n)+tlevSTD(i,k)
-           usumSTD(i,k,n)=usumSTD(i,k,n)+ulevSTD(i,k)
-           vsumSTD(i,k,n)=vsumSTD(i,k,n)+vlevSTD(i,k)
-           wsumSTD(i,k,n)=wsumSTD(i,k,n)+wlevSTD(i,k)
-           phisumSTD(i,k,n)=phisumSTD(i,k,n)+philevSTD(i,k)
-           qsumSTD(i,k,n)=qsumSTD(i,k,n)+qlevSTD(i,k)
-           rhsumSTD(i,k,n)=rhsumSTD(i,k,n)+rhlevSTD(i,k)
-           uvsumSTD(i,k,n)=uvsumSTD(i,k,n)+uvSTD(i,k)
-           vqsumSTD(i,k,n)=vqsumSTD(i,k,n)+vqSTD(i,k)
-           vTsumSTD(i,k,n)=vTsumSTD(i,k,n)+vTSTD(i,k)
-           wqsumSTD(i,k,n)=wqsumSTD(i,k,n)+wqSTD(i,k)
-           vphisumSTD(i,k,n)=vphisumSTD(i,k,n)+vphiSTD(i,k)
-           wTsumSTD(i,k,n)=wTsumSTD(i,k,n)+wTSTD(i,k)
-           u2sumSTD(i,k,n)=u2sumSTD(i,k,n)+u2STD(i,k)
-           v2sumSTD(i,k,n)=v2sumSTD(i,k,n)+v2STD(i,k)
-           T2sumSTD(i,k,n)=T2sumSTD(i,k,n)+T2STD(i,k)
-           O3sumSTD(i,k,n)=O3sumSTD(i,k,n)+O3STD(i,k)
-           IF (read_climoz==2) 
-     &          O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)+O3daySTD(i,k)
 
-          ENDIF 
-         ENDDO !i
-        ENDDO !k
-c
-       ENDIF !MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0
-c
-      ENDDO !n
-c
-      RETURN
-      END  
+SUBROUTINE undefstd(itap, read_climoz)
+  USE netcdf
+  USE dimphy
+  USE phys_state_var_mod ! Variables sauvegardees de la physique
+  IMPLICIT NONE
+  include "clesphys.h"
+
+  ! ====================================================================
+
+  ! I. Musat : 09.2004
+
+  ! Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef))
+  ! ou la variable tlevSTD est bien definie (.NE.missing_val),
+  ! et
+  ! * de la somme de tlevSTD => tsumSTD
+
+  ! nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
+  ! ! de la physique
+  ! nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
+  ! ! de la physique
+  ! nout=3 !var. mensuelle "NMC" moyennee toutes les ecrit_hf
+
+
+  ! NB: mettre "inst(X)" dans le write_hist*NMC.h !
+  ! ====================================================================
+
+  ! ym#include "dimensions.h"
+  ! ym      integer jjmp1
+  ! ym      parameter (jjmp1=jjm+1-1/jjm)
+  ! ym#include "dimphy.h"
+  ! variables Input
+
+  ! INTEGER nlevSTD, klevSTD, itap
+  ! PARAMETER(klevSTD=17)
+  INTEGER itap
+  ! REAL dtime
+
+  ! variables locales
+  ! INTEGER i, k, nout, n
+  ! PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
+  INTEGER i, k, n
+  ! REAL freq_calNMC(nout)
+  INTEGER read_climoz
+
+  ! variables Output
+  ! REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
+  ! LOGICAL oknondef(klon,klevSTD,nout)
+  ! REAL tnondef(klon,klevSTD,nout)
+
+  REAL missing_val
+
+  missing_val = nf90_fill_real
+
+  DO n = 1, nout
+
+
+    ! calcul variables tous les freq_calNMC(n)/dtime pas de temps
+    ! de la physique
+
+    IF (mod(itap,nint(freq_calnmc(n)/dtime))==0) THEN
+      DO k = 1, nlevstd
+        DO i = 1, klon
+          IF (tlevstd(i,k)==missing_val) THEN
+            ! IF(oknondef(i,k,n)) THEN
+            tnondef(i, k, n) = tnondef(i, k, n) + 1.
+            ! ENDIF !oknondef(i,k)
+
+          ELSE IF (tlevstd(i,k)/=missing_val) THEN
+            tsumstd(i, k, n) = tsumstd(i, k, n) + tlevstd(i, k)
+            usumstd(i, k, n) = usumstd(i, k, n) + ulevstd(i, k)
+            vsumstd(i, k, n) = vsumstd(i, k, n) + vlevstd(i, k)
+            wsumstd(i, k, n) = wsumstd(i, k, n) + wlevstd(i, k)
+            phisumstd(i, k, n) = phisumstd(i, k, n) + philevstd(i, k)
+            qsumstd(i, k, n) = qsumstd(i, k, n) + qlevstd(i, k)
+            rhsumstd(i, k, n) = rhsumstd(i, k, n) + rhlevstd(i, k)
+            uvsumstd(i, k, n) = uvsumstd(i, k, n) + uvstd(i, k)
+            vqsumstd(i, k, n) = vqsumstd(i, k, n) + vqstd(i, k)
+            vtsumstd(i, k, n) = vtsumstd(i, k, n) + vtstd(i, k)
+            wqsumstd(i, k, n) = wqsumstd(i, k, n) + wqstd(i, k)
+            vphisumstd(i, k, n) = vphisumstd(i, k, n) + vphistd(i, k)
+            wtsumstd(i, k, n) = wtsumstd(i, k, n) + wtstd(i, k)
+            u2sumstd(i, k, n) = u2sumstd(i, k, n) + u2std(i, k)
+            v2sumstd(i, k, n) = v2sumstd(i, k, n) + v2std(i, k)
+            t2sumstd(i, k, n) = t2sumstd(i, k, n) + t2std(i, k)
+            o3sumstd(i, k, n) = o3sumstd(i, k, n) + o3std(i, k)
+            IF (read_climoz==2) o3daysumstd(i, k, n) = o3daysumstd(i, k, n) + &
+              o3daystd(i, k)
+
+          END IF
+        END DO !i
+      END DO !k
+
+    END IF !MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0
+
+  END DO !n
+
+  RETURN
+END SUBROUTINE undefstd