Changeset 2117 for trunk

Timestamp:

Mar 18, 2019, 4:27:59 PM (6 years ago)

Author:

jvatant

Message:

Follow-up of r2116, fix errors for runs mufi+chem in tracer gestion
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (9 diffs)

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -522 +522 @@
             ctimestep = ptimestep*REAL(ichim)
 
-            zdqchi(:,:,:)  = 0.0
-
          endif
 
@@ -1094 +1092 @@
             ! Microphysics condensation for 2D fields to sent non-saturated fields to photochem
             if ( callclouds .and. moyzon_ch .and. mod(icount-1,ichim).eq.0 ) then
-              zdqfibar(:,:,:) = 0.0 ! We work in zonal average -> forget processes other than condensation
+              zdqfibar(:,:,:) = 0.D0 ! We work in zonal average -> forget processes other than condensation
               call calmufi(ptimestep,zplevbar,zzlevbar,zplaybar,zzlaybar, &
                            gzlat,ztfibar,zqfibar,zdqfibar,zdqmufibar)
@@ -1135 +1133 @@
             call calc_ysat(ngrid,nlayer,pplay/100.0,pt,ysat) ! Compute saturation profiles for every grid point
 
-            dyccond(:,:,:) = 0.0 ! Default value -> no condensation
+            dyccond(:,:,:) = 0.D0 ! Default value -> no condensation
 
             do iq=1,nkim
@@ -1142 +1140 @@
             enddo
 
-            if ( callclouds ) dyccond(:,:,ices_indx) = 0.0 ! Condensation have been calculated in the cloud microphysics 
+            if ( callclouds ) dyccond(:,:,ices_indx) = 0.D0 ! Condensation have been calculated in the cloud microphysics 
 
             do iq=1,nkim
@@ -1157 +1155 @@
               call calc_ysat(ngrid,nlayer,zplaybar/100.0,ztfibar,ysat) ! Compute saturation profiles for every grid point for the zon-ave fields
 
-              dyccondbar(:,:,:) = 0.0 ! Default value -> no condensation
+              dyccondbar(:,:,:) = 0.D0 ! Default value -> no condensation
               
               do iq = 1,nkim
@@ -1164 +1162 @@
               enddo
 
-              if ( callclouds ) dyccondbar(:,:,ices_indx) = 0.0 ! Condensation have been calculated in the cloud microphysics 
+              if ( callclouds ) dyccondbar(:,:,ices_indx) = 0.D0 ! Condensation have been calculated in the cloud microphysics 
 
               do iq=1,nkim
@@ -1172 +1170 @@
               ! Pseudo-evap ( forcing constant surface humidity )
               do ig=1,ngrid
-                 if ( ychimbar(ig,1,7+nmicro) .lt. botCH4 ) ychimbar(ig,1,7+nmicro) = botCH4
+                 if ( ychimbar(ig,1,7) .lt. botCH4 ) ychimbar(ig,1,7) = botCH4
               enddo
 
@@ -1198 +1196 @@
             ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
             do ig=1,ngrid
-               if ( (ychim(ig,1,7+nmicro)+dycchi(ig,1,7+nmicro)*ptimestep) .lt. botCH4 ) then ! +dycchi because ychim not yet updated
-                  dycevapCH4(ig) = ( -ychim(ig,1,7+nmicro)+botCH4 ) / ptimestep - dycchi(ig,1,7+nmicro)
+               if ( (ychim(ig,1,7)+dycchi(ig,1,7)*ptimestep) .lt. botCH4 ) then ! +dycchi because ychim not yet updated
+                  dycevapCH4(ig) = ( -ychim(ig,1,7)+botCH4 ) / ptimestep - dycchi(ig,1,7)
                else
-                  dycevapCH4(ig) = 0.0
+                  dycevapCH4(ig) = 0.D0
                endif
             enddo
@@ -1209 +1207 @@
             ! v. Updates and positivity check
             ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+            zdqchi(:,:,:)  = 0.D0 ! -> for the microphysical tracers we need to update to 0 each time ;-)
+
             do iq=1,nkim
               zdqchi(:,:,iq+nmicro) = dycchi(:,:,iq)*rat_mmol(iq+nmicro) ! convert tendencies to mass mixing ratio
 
               where( (pq(:,:,iq+nmicro) + ( pdq(:,:,iq+nmicro)+zdqchi(:,:,iq+nmicro) )*ptimestep ) .LT. 0.)  & ! When using zonal means we set the same tendency
-                      zdqchi(:,:,iq+nmicro) = 1.e-30 - pdq(:,:,iq+nmicro) - pq(:,:,iq+nmicro)/ptimestep        ! everywhere in longitude -> could lead to negs !
+                      zdqchi(:,:,iq+nmicro) = 1.D-30 - pdq(:,:,iq+nmicro) - pq(:,:,iq+nmicro)/ptimestep        ! everywhere in longitude -> could lead to negs !
             enddo