Changeset 2298 for LMDZ5/branches/testing/libf/dyn3dmem/iniacademic_loc.F90

Timestamp:

Jun 14, 2015, 9:13:32 PM (9 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes -r2237:2291 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/dyn3dmem/iniacademic_loc.F90 (modified) (3 diffs)

LMDZ5/branches/testing/libf/dyn3dmem/iniacademic_loc.F90

@@ -7 +7 @@
   use exner_hyb_m, only: exner_hyb
   use exner_milieu_m, only: exner_milieu
-  USE infotrac, ONLY : nqtot
+  USE infotrac, ONLY: nqtot,niso_possibles,ok_isotopes,iqpere,ok_iso_verif,tnat,alpha_ideal, &
+        & iqiso,phase_num,iso_indnum,iso_num,zone_num
   USE control_mod, ONLY: day_step,planet_type
   USE parallel_lmdz, ONLY: ijb_u, ije_u, ijb_v, ije_v
@@ -110 +111 @@
   ztot0      = 0.
   stot0      = 0.
-  ang0       = 0.
+  ang0       = 0.      
 
   if (llm == 1) then
@@ -269 +270 @@
         if (planet_type=="earth") then
            ! Earth: first two tracers will be water
+
            do i=1,nqtot
               if (i == 1) q(ijb_u:ije_u,:,i)=1.e-10
               if (i == 2) q(ijb_u:ije_u,:,i)=1.e-15
               if (i.gt.2) q(ijb_u:ije_u,:,i)=0.
+
+              ! CRisi: init des isotopes
+              ! distill de Rayleigh très simplifiée
+              if (ok_isotopes) then
+                if ((iso_num(i).gt.0).and.(zone_num(i).eq.0)) then          
+                   q(ijb_u:ije_u,:,i)=q(ijb_u:ije_u,:,iqpere(i))       &
+      &                  *tnat(iso_num(i))                             &
+      &                  *(q(ijb_u:ije_u,:,iqpere(i))/30.e-3)                              &
+     &                   **(alpha_ideal(iso_num(i))-1)
+                endif                
+                if ((iso_num(i).gt.0).and.(zone_num(i).eq.1)) then
+                  q(ijb_u:ije_u,:,i)=q(ijb_u:ije_u,:,iqiso(iso_indnum(i),phase_num(i)))
+                endif
+              endif !if (ok_isotopes) then
+
            enddo
         else
            q(ijb_u:ije_u,:,:)=0
         endif ! of if (planet_type=="earth")
+
+        if (ok_iso_verif) then
+           call check_isotopes(q,ijb_u,ije_u,'iniacademic_loc')
+        endif !if (ok_iso_verif) then
 
         ! add random perturbation to temperature