Changeset 2270 for LMDZ5/trunk/libf/dyn3d/dynetat0.F

Timestamp:

May 7, 2015, 5:45:04 PM (10 years ago)

Author:

crisi

Message:

Adding isotopes in the dynamics and more generally tracers of tracers.
CRisi

File:

• LMDZ5/trunk/libf/dyn3d/dynetat0.F (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/dynetat0.F

@@ -297 +297 @@
            write(lunout,*)"          Il est donc initialise a zero"
            q(:,:,:,iq)=0.
+
+           ! CRisi: pour les isotopes, on peut faire init théorique
+           ! distill de Rayleigh très simplifiée
+           if (ok_isotopes) then
+              if ((iso_num(iq).gt.0).and.(zone_num(iq).eq.0)) then
+                q(:,:,:,iq)=q(:,:,:,iqpere(iq))                         
+     &                   *tnat(iso_num(iq))                             
+     &                   *(q(:,:,:,iqpere(iq))/30.e-3)                  
+     &                   **(alpha_ideal(iso_num(iq))-1)
+              endif
+              if ((iso_num(iq).gt.0).and.(zone_num(iq).eq.1)) then
+                  q(:,:,:,iq)=q(:,:,:,iqiso(iso_indnum(iq),
+     &                   phase_num(iq)))
+              endif  
+           endif !if (ok_isotopes) then 
         ELSE
            ierr = NF90_GET_VAR(nid, nvarid, q(:,:,:,iq))