Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/qminimum_loc.f90

Timestamp:

Jul 24, 2024, 1:17:08 PM (2 months ago)

Author:

abarral

Message:

Rename modules in misc from *_mod > lmdz_*
Put cbrt.f90, ch*.f90, pch*.f90 in new lmdz_libmath_pch.f90

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/qminimum_loc.f90 (modified) (10 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/qminimum_loc.f90

@@ -55 +55 @@
 !$OMP BARRIER
 
-  ! !write(lunout,*) 'qminimum 52: entree'
+  !write(lunout,*) 'qminimum 52: entree'
   IF(first) THEN
      iq_vap = strIdx(tracers(:)%name, addPhase('H2O', 'g'))
@@ -83 +83 @@
   ENDDO
 
-  ! !write(lunout,*) 'qminimum 57'
+  !write(lunout,*) 'qminimum 57'
   DO k = 1, llm
 !$OMP DO SCHEDULE(STATIC)
@@ -103 +103 @@
   ! le defaut en prennant de l'eau vapeur de la couche au-dessous.
   !
-  ! !write(lunout,*) 'qminimum 81'
+  !write(lunout,*) 'qminimum 81'
   DO k = llm, 2, -1
   !cc      zx_abc = dpres(k) / dpres(k-1)
@@ -127 +127 @@
   ! doit imprimer un message d'avertissement (saturation possible).
   !
-  ! !write(lunout,*) 'qminimum 106'
+  !write(lunout,*) 'qminimum 106'
   nb_pump=0
 !$OMP DO SCHEDULE(STATIC)
@@ -149 +149 @@
   ENDIF
 
-  ! !write(lunout,*) 'qminimum 128'
+  !write(lunout,*) 'qminimum 128'
   if (niso > 0) then
-          ! !write(lunout,*) 'qminimum 140'
-  ! ! CRisi: traiter de même les traceurs d'eau
-  ! ! Mais il faut les prendre à l'envers pour essayer de conserver la
-  ! ! masse.
-  ! ! 1) pompage dans le sol
-  ! ! On suppose que ce pompage se fait sans isotopes -> on ne modifie
-  ! ! rien ici et on croise les doigts pour que ça ne soit pas trop
-  ! ! génant
-  ! ! en fait, si, c'est genant quand les isotopes doivent eux même transporter des
-  ! ! traceurs -> apporter aussi un peu d'isotopes... Combien?
-  ! ! Essayer tnat/2 = -500 permil? C'est déjà mieux que -1000
-  ! ! permil...
-  ! ! pb: que faire pour les traceurs?
+          !write(lunout,*) 'qminimum 140'
+  ! CRisi: traiter de même les traceurs d'eau
+  ! Mais il faut les prendre à l'envers pour essayer de conserver la
+  ! masse.
+  ! 1) pompage dans le sol
+  ! On suppose que ce pompage se fait sans isotopes -> on ne modifie
+  ! rien ici et on croise les doigts pour que ça ne soit pas trop
+  ! génant
+  ! en fait, si, c'est genant quand les isotopes doivent eux même transporter des
+  ! traceurs -> apporter aussi un peu d'isotopes... Combien?
+  ! Essayer tnat/2 = -500 permil? C'est déjà mieux que -1000
+  ! permil...
+  ! pb: que faire pour les traceurs?
 !$OMP DO SCHEDULE(STATIC)
   DO i = ijb, ije
@@ -172 +172 @@
 !$OMP END DO NOWAIT
 
-  ! ! 2) transfert de vap vers les couches plus hautes
-  ! !write(lunout,*) 'qminimum 158'
+  ! 2) transfert de vap vers les couches plus hautes
+  !write(lunout,*) 'qminimum 158'
   do k=2,llm
 !$OMP DO SCHEDULE(STATIC)
     DO i = ijb, ije
       if (zx_defau_diag(i,k,1)>0.0) then
-          ! ! on ajoute la vapeur en k
+          ! on ajoute la vapeur en k
           !  write(lunout,*) 'i,k,q_follow(i,k-1,ivap)=',
   ! :                 i,k,q_follow(i,k-1,1)
@@ -218 +218 @@
           endif
 
-          ! ! et on la retranche en k-1
+          ! et on la retranche en k-1
            q(i,k-1,iqIsoPha(ixt,iq_vap)) = &
                  q(i,k-1,iqIsoPha(ixt,iq_vap)) &
@@ -257 +257 @@
 
 
-    ! ! 3) transfert d'eau de la vapeur au liquide
-    ! !write(*,*) 'qminimum 164'
+    ! 3) transfert d'eau de la vapeur au liquide
+    !write(*,*) 'qminimum 164'
     do k=1,llm
 !$OMP DO SCHEDULE(STATIC)
@@ -264 +264 @@
       if (zx_defau_diag(i,k,2)>0.0) then
 
-          ! ! on ajoute eau liquide en k en k
+          ! on ajoute eau liquide en k en k
           do ixt=1,ntiso
            q(i,k,iqIsoPha(ixt,iq_liq))=q(i,k,iqIsoPha(ixt,iq_liq)) &
                  +zx_defau_diag(i,k,2) &
                  *q(i,k,iqIsoPha(ixt,iq_vap))/q_follow(i,k,1)
-          ! ! et on la retranche à la vapeur en k
+          ! et on la retranche à la vapeur en k
            q(i,k,iqIsoPha(ixt,iq_vap))=q(i,k,iqIsoPha(ixt,iq_vap)) &
                  -zx_defau_diag(i,k,2) &
@@ -286 +286 @@
 
   endif !if (niso > 0) then
-  ! !write(*,*) 'qminimum 188'
+  !write(*,*) 'qminimum 188'
 !$OMP BARRIER