Changeset 1806

Timestamp:

Oct 20, 2017, 11:22:19 AM (7 years ago)

Author:

bclmd

Message:

inichim_newstart with Mars case by default

File:

• trunk/LMDZ.GENERIC/libf/aeronostd/inichim_newstart.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/aeronostd/inichim_newstart.F90

@@ -448 +448 @@
          endif
          flagnitro = .true.
-         nspe = 35
+         nspe = 18
       endif
 
@@ -489 +489 @@
 
       niq(1) = igcm_co2
-!      niq(2) = igcm_ar
-      niq(2) = igcm_n2
-      niq(3) = igcm_o2
-      niq(4) = igcm_co
-      niq(5) = igcm_o
-      niq(6) = igcm_h2
+      niq(2) = igcm_ar
+      niq(3) = igcm_n2
+      niq(4) = igcm_o2
+      niq(5) = igcm_co
+      niq(6) = igcm_o
+      niq(7) = igcm_h2
 
 ! minor species:
 
-      niq(7)  = igcm_h
-      niq(8)  = igcm_oh
-      niq(9) = igcm_ho2
-      niq(10) = igcm_h2o_vap
-      niq(11) = igcm_h2o2
-      niq(12) = igcm_o1d 
-      niq(13) = igcm_o3
+      niq(8)  = igcm_h
+      niq(9)  = igcm_oh
+      niq(10) = igcm_ho2
+      niq(11) = igcm_h2o_vap
+      niq(12) = igcm_h2o2
+      niq(13) = igcm_o1d 
+      niq(14) = igcm_o3
 
 ! nitrogen species:
 
       if (flagnitro) then
-         niq(14) = igcm_n
-         niq(15) = igcm_no
-         niq(16) = igcm_no2
-         niq(17) = igcm_n2d         
+         niq(15) = igcm_n
+         niq(16) = igcm_no
+         niq(17) = igcm_no2
+         niq(18) = igcm_n2d         
       end if
 
@@ -519 +519 @@
 
 ! carbon species:
-      niq(18) = igcm_ch4
-      niq(19) = igcm_ch3
-      niq(20) = igcm_ch
-      niq(21) = igcm_1ch2
-      niq(22) = igcm_3ch2
-      niq(23) = igcm_cho
-      niq(24) = igcm_ch2o
-      niq(25) = igcm_ch3o
-      niq(26) = igcm_c
-      niq(27) = igcm_c2
-      niq(28) = igcm_c2h
-      niq(29) = igcm_c2h2
-      niq(30) = igcm_c2h3
-      niq(31) = igcm_c2h4
-      niq(32) = igcm_c2h6
-      niq(33) = igcm_ch2co
-      niq(34) = igcm_ch3co
-      niq(35) = igcm_hcaer
+!      niq(18) = igcm_ch4
+!      niq(19) = igcm_ch3
+!      niq(20) = igcm_ch
+!      niq(21) = igcm_1ch2
+!      niq(22) = igcm_3ch2
+!      niq(23) = igcm_cho
+!      niq(24) = igcm_ch2o
+!      niq(25) = igcm_ch3o
+!      niq(26) = igcm_c
+!      niq(27) = igcm_c2
+!      niq(28) = igcm_c2h
+!      niq(29) = igcm_c2h2
+!      niq(30) = igcm_c2h3
+!      niq(31) = igcm_c2h4
+!      niq(32) = igcm_c2h6
+!      niq(33) = igcm_ch2co
+!      niq(34) = igcm_ch3co
+!      niq(35) = igcm_hcaer
 
 
@@ -610 +610 @@
       end if
 
-
-      if (1.eq.1) then
+! initialization for the early eath
+      if (1.eq.0) then
         do l = 1,nalt
          vmrinit(l,:)=0.0