Changeset 4124 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

Timestamp:

Apr 8, 2022, 4:47:04 PM (3 years ago)

Author:

dcugnet

Message:

Remove solsym, ok_isotopes (=niso>0), ok_isotrac (=nzone>0)

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -40 +40 @@
     USE indice_sol_mod
     USE infotrac, ONLY: iso_num, iso_indnum
-    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac, nqCO2, ok_isotopes, indnum_fn_num
+    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac, nqCO2, indnum_fn_num
     USE readTracFiles_mod, ONLY: addPhase
     USE strings_mod,  ONLY: strIdx, strStack, int2str
@@ -126 +126 @@
 
 #ifdef ISO
-    USE infotrac_phy, ONLY:  &
-        iqiso,ok_isotrac,niso, ntraciso
-     USE isotopes_mod, ONLY: iso_eau,iso_HDO,iso_O18,iso_O17,iso_HTO, &
+    USE infotrac_phy, ONLY: iqiso,niso, ntraciso, nzone
+    USE isotopes_mod, ONLY: iso_eau,iso_HDO,iso_O18,iso_O17,iso_HTO, &
         & bidouille_anti_divergence,ok_bidouille_wake, &
         & modif_ratqs,essai_convergence,iso_init,ridicule_rain,tnat, &
         & ridicule,ridicule_snow
-     USE isotopes_routines_mod, ONLY: iso_tritium
+    USE isotopes_routines_mod, ONLY: iso_tritium
 #ifdef ISOVERIF
     USE isotopes_verif_mod, ONLY: errmax,errmaxrel, &
@@ -1595 +1594 @@
 ! dyn3dmem et physiq
 #ifdef ISO
-    write(*,*) 'physiq 1846a: ok_isotopes,ntraciso,niso=',ok_isotopes,ntraciso,niso
-    if (.not.ok_isotopes) then
-      CALL abort_gcm('physiq 1756','options iso incompatibles',1)
-    endif
+    write(*,*) 'physiq 1846a: ok_isotopes,ntraciso,niso=',niso>0,ntraciso,niso
+    IF(niso  <= 0) CALL abort_gcm('physiq 1756','options iso incompatibles',1)
 #ifdef ISOTRAC
-    if (.not.ok_isotrac) then
-      CALL abort_gcm('physiq 1758','options isotrac incompatibles',1)
-    endif    
+    IF(nzone <= 0) CALL abort_gcm('physiq 1758','options isotrac incompatibles',1)
 #else
-! #ifdef ISOTRAC
-    if (ok_isotrac) then
-      CALL abort_gcm('physiq 1762','options isotrac incompatibles',1)
-    endif
-#endif
-!! #ifdef ISOTRAC
-! -> on supprime opion ISOTRAC, tout passe par ok_isotrac
+    IF(nzone  > 0) CALL abort_gcm('physiq 1762','options isotrac incompatibles',1)
+#endif
 #else
-! #ifdef ISO
-    if (ok_isotopes) then
-      CALL abort_gcm('physiq 1772','options iso incompatibles',1)
-    endif
-#endif
-! #ifdef ISO
+    if(niso   > 0) CALL abort_gcm('physiq 1772','options iso incompatibles',1)
+#endif
 
 #ifdef ISO
@@ -1623 +1609 @@
 #ifdef ISOVERIF
            write(*,*) 'physiq 1366: call iso_init'
-           write(*,*) 'ok_isotopes=',ok_isotopes
-#endif
-        if (ok_isotopes) then
-           call iso_init()
-        endif
+           write(*,*) 'ok_isotopes=',niso > 0
+#endif
+        if (niso > 0) call iso_init()
 #ifdef ISOTRAC
-if (ok_isotrac) then
+IF(nzone > 0) then
         write(*,*) 'physiq 1416: call iso_traceurs_init'
         call iso_traceurs_init()
@@ -2424 +2408 @@
 #ifdef ISOVERIF
 !    write(*,*) 'physiq 1847: qx(1,1,:)=',qx(1,1,:)
-    write(*,*) 'physiq 1846b: ok_isotopes,ntraciso,niso=',ok_isotopes,ntraciso,niso
+    write(*,*) 'physiq 1846b: ok_isotopes,ntraciso,niso=',niso>0,ntraciso,niso
 #endif
     do ixt=1,ntraciso
@@ -2451 +2435 @@
 #endif
 ! #ifdef ISO
-
     !
     !--OB mass fixer 
@@ -2521 +2504 @@
     ! Diagnostiquer la tendance dynamique
 #ifdef ISOVERIF
-      !write(*,*) 'physiq tmp 2010: ancien_ok=',ancien_ok       
       if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then  
         do i=1,klon 
@@ -6821 +6803 @@
             endif ! if (iso_eau.gt.0) then
 #ifdef ISOTRAC
-            if (ok_isotrac) then      
-            call iso_verif_traceur(xt_ancien(1,i,k),'physiq 4802')
-            endif !if (ok_isotrac) then 
+            IF(nzone > 0) CALL iso_verif_traceur(xt_ancien(1,i,k),'physiq 4802')
 #endif          
           enddo