Changeset 2298 for LMDZ5/branches/testing/libf/dyn3d/dynetat0.F

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/dyn3d/dynetat0.F (modified) (1 diff)

LMDZ5/branches/testing/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))