Changeset 1910 for LMDZ5/branches/testing/libf/phylmd/calcul_STDlev.h

Timestamp:

Nov 29, 2013, 9:40:25 AM (11 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r1860:1909 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (2 props)
• libf (modified) (1 prop)
• libf/phylmd (modified) (1 prop)
• libf/phylmd/calcul_STDlev.h (modified) (15 diffs, 1 prop)

LMDZ5/branches/testing/libf/phylmd/calcul_STDlev.h

@@ -1 @@
-c
-c $Header$
-c
-cIM on initialise les variables 
-c
+!
+! $Header$
+!
+!IM on initialise les variables 
+!
         missing_val=nf90_fill_real
-c
-cIM freq_moyNMC = frequences auxquelles on moyenne les champs accumules 
-cIM               sur les niveaux de pression standard du NMC
+!
+!IM freq_moyNMC = frequences auxquelles on moyenne les champs accumules 
+!IM               sur les niveaux de pression standard du NMC
       DO n=1, nout
        freq_moyNMC(n)=freq_outNMC(n)/freq_calNMC(n)
       ENDDO
-c
+!
         CALL ini_undefSTD(itap,freq_outNMC)
-c
-cIM on interpole les champs sur les niveaux STD de pression 
-cIM a chaque pas de temps de la physique
-c
-c-------------------------------------------------------c
-c positionnement de l'argument logique a .false.        c
-c pour ne pas recalculer deux fois la meme chose !      c
-c a cet effet un appel a plevel_new a ete deplace       c
-c a la fin de la serie d'appels                         c
-c la boucle 'DO k=1, nlevSTD' a ete internalisee        c
-c dans plevel_new, d'ou la creation de cette routine... c
-c-------------------------------------------------------c
-c
-        CALL plevel_new(klon,klev,nlevSTD,.true.,pplay,rlevSTD,
-     &              t_seri,tlevSTD)
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             u_seri,ulevSTD)
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             v_seri,vlevSTD)
-c
-
-c
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zphi/RG,philevSTD)
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             qx(:,:,ivap),qlevSTD)
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_rh*100.,rhlevSTD)
-c
+!
+!IM on interpole les champs sur les niveaux STD de pression 
+!IM a chaque pas de temps de la physique
+!
+!-------------------------------------------------------c
+! positionnement de l'argument logique a .false.        c
+! pour ne pas recalculer deux fois la meme chose !      c
+! a cet effet un appel a plevel_new a ete deplace       c
+! a la fin de la serie d'appels                         c
+! la boucle 'DO k=1, nlevSTD' a ete internalisee        c
+! dans plevel_new, d'ou la creation de cette routine... c
+!-------------------------------------------------------c
+!
+        CALL plevel_new(klon,klev,nlevSTD,.true.,pplay,rlevSTD, &
+                    t_seri,tlevSTD)
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   u_seri,ulevSTD)
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   v_seri,vlevSTD)
+!
+
+!
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zphi/RG,philevSTD)
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   qx(:,:,ivap),qlevSTD)
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_rh*100.,rhlevSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -47 +47 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,uvSTD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,uvSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -55 +55 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,vqSTD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,vqSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -63 +63 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,vTSTD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,vTSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -71 +71 @@
          ENDDO !i 
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,wqSTD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,wqSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -79 +79 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,vphiSTD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,vphiSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -87 +87 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,wTSTD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,wTSTD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -95 +95 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,u2STD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,u2STD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -103 +103 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,v2STD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,v2STD)
+!
         DO l=1, klev
          DO i=1, klon
@@ -111 +111 @@
          ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,T2STD)
-
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,T2STD)
+
+!
       zx_tmp_fi3d(:,:)=wo(:,:,1) * dobson_u * 1e3 / zmasse / rmo3 * rmd
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,O3STD)
-c
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,O3STD)
+!
       if (read_climoz == 2) THEN
       zx_tmp_fi3d(:,:)=wo(:,:,2) * dobson_u * 1e3 / zmasse / rmo3 * rmd
-        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD,
-     &             zx_tmp_fi3d,O3daySTD)
+        CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, &
+                   zx_tmp_fi3d,O3daySTD)
       endif
-c
+!
         DO l=1, klev
         DO i=1, klon
@@ -130 +130 @@
         ENDDO !i
         ENDDO !l
-        CALL plevel_new(klon,klev,nlevSTD,.true.,zx_tmp_fi3d,rlevSTD,
-     &             omega,wlevSTD)
-c
-cIM on somme les valeurs toutes les freq_calNMC secondes
-c
+        CALL plevel_new(klon,klev,nlevSTD,.true.,zx_tmp_fi3d,rlevSTD, &
+                   omega,wlevSTD)
+!
+!IM on somme les valeurs toutes les freq_calNMC secondes
+!
        CALL undefSTD(itap,freq_calNMC, read_climoz)
-c
-cIM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes)
-c
+!
+!IM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes)
+!
        CALL moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
-c
-       CALL plevel(klon,klev,.true.,pplay,50000.,
-     &              zphi/RG,geo500)
-
-cIM on interpole a chaque pas de temps le SWup(clr) et SWdn(clr) a 200 hPa
-c
-      CALL plevel(klon,klevp1,.true.,paprs,20000.,
-     $     swdn0,SWdn200clr)
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     swdn,SWdn200)
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     swup0,SWup200clr)
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     swup,SWup200)
-c
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     lwdn0,LWdn200clr)
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     lwdn,LWdn200)
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     lwup0,LWup200clr)
-      CALL plevel(klon,klevp1,.false.,paprs,20000.,
-     $     lwup,LWup200)
-c
+!
+       CALL plevel(klon,klev,.true.,pplay,50000., &
+                    zphi/RG,geo500)
+
+!IM on interpole a chaque pas de temps le SWup(clr) et SWdn(clr) a 200 hPa
+!
+      CALL plevel(klon,klevp1,.true.,paprs,20000., &
+           swdn0,SWdn200clr)
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           swdn,SWdn200)
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           swup0,SWup200clr)
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           swup,SWup200)
+!
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           lwdn0,LWdn200clr)
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           lwdn,LWdn200)
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           lwup0,LWup200clr)
+      CALL plevel(klon,klevp1,.false.,paprs,20000., &
+           lwup,LWup200)
+!
       twriteSTD(:,:,1)=tsumSTD(:,:,1)
       qwriteSTD(:,:,1)=qsumSTD(:,:,1)
@@ -195 +195 @@
       vwriteSTD(:,:,4)=vlevSTD(:,:)
       wwriteSTD(:,:,4)=wlevSTD(:,:)
-c
-cIM initialisation 5eme fichier de sortie 
+!
+!IM initialisation 5eme fichier de sortie 
       twriteSTD(:,:,5)=tlevSTD(:,:)
       qwriteSTD(:,:,5)=qlevSTD(:,:)
@@ -204 +204 @@
       vwriteSTD(:,:,5)=vlevSTD(:,:)
       wwriteSTD(:,:,5)=wlevSTD(:,:)
-c
-cIM initialisation 6eme fichier de sortie 
+!
+!IM initialisation 6eme fichier de sortie 
       twriteSTD(:,:,6)=tlevSTD(:,:)
       qwriteSTD(:,:,6)=qlevSTD(:,:)
@@ -213 +213 @@
       vwriteSTD(:,:,6)=vlevSTD(:,:)
       wwriteSTD(:,:,6)=wlevSTD(:,:)
-cIM for NMC files
+!IM for NMC files
       DO n=1, nlevSTD3
        DO k=1, nlevSTD
@@ -227 +227 @@
        ENDDO 
       ENDDO 
-c
+!
       DO n=1, nlevSTD8
        DO k=1, nlevSTD