Changeset 2217 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Jan 14, 2020, 8:50:50 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in newstart "composition" option; added option "q=factor" in
newstart to multiply a tracer mixing ratio by a constant factor.
EM

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -500 +500 @@
       write(*,*) 'qname        : change tracer name'
       write(*,*) 'q=0          : ALL tracer =zero'
+      write(*,*) 'q=factor     : change tracer value by a multiplicative
+     & factor'
       write(*,*) 'q=x          : give a specific uniform value to one
      $ tracer'
@@ -770 +772 @@
              ENDDO
            ENDDO
+
+c       q=factor : change value of tracer by a multiplicative factor
+c       ------------------------------------------------------------
+        else if (trim(modif) .eq. 'q=factor') then
+             write(*,*) 'Which tracer do you want to modify ?'
+             do iq=1,nqtot
+               write(*,*)iq,' : ',trim(tname(iq))
+             enddo
+             write(*,*) '(choose between 1 and ',nqtot,')'
+             read(*,*) iq 
+             if ((iq.lt.1).or.(iq.gt.nqtot)) then
+               ! wrong value for iq, go back to menu
+               write(*,*) "wrong input value:",iq
+               cycle
+             endif
+             write(*,*)"factor to multiply current mixing ratio by?"
+             read(*,*) val
+             
+             q(1:iip1,1:jjp1,1:llm,iq)=q(1:iip1,1:jjp1,1:llm,iq)*val
+             qsurf(1:ngridmx,iq)=qsurf(1:ngridmx,iq)*val
 
 c       q=x : initialise tracer manually 
@@ -1019 +1041 @@
        !   molar mass of non-co2 atmosphere measured by MSL at Ls ~184
        
-           write(*,*)
+          else
+            ! Assume co2/n2/ar/o2/co are available
+            Mair_old=(q(iref,jref,1,igcm_co2)*mmol(igcm_co2)
+     &               +q(iref,jref,1,igcm_n2)*mmol(igcm_n2)
+     &               +q(iref,jref,1,igcm_ar)*mmol(igcm_ar)
+     &               +q(iref,jref,1,igcm_o2)*mmol(igcm_o2)
+     &               +q(iref,jref,1,igcm_co)*mmol(igcm_co))/Smmr
+          end if
+
+          write(*,*)
      &      "Air molar mass (g/mol) at reference site= ",Mair_old
-          end if
 
         ! Ask for new volume mixing ratio at reference site
@@ -1114 +1144 @@
            enddo
           enddo
-
-
-!          do l=1,llm
-!           do j=1,jjp1
-!            do i=1,iip1
-!              Smmr_old = 0.
-!              Smmr_new = 0.
-!              do iq=1,nqtot  
-!                if ((iq.eq.igcm_co2).or.(iq.eq.igcm_n2)
-!     &          .or.(iq.eq.igcm_ar)
-!     &          .or. (iq.eq.igcm_o2).or.(iq.eq.igcm_co)) then
-!                   Smmr_old = Smmr_old + q(i,j,l,iq) ! sum of old mmr 
-!                   q(i,j,l,iq)=q(i,j,l,iq)*coefvmr(iq)*Mair_old/Mair_new
-!                   Smmr_new = Smmr_new + q(i,j,l,iq) ! sum of new mmr
-!                end if 
-!              enddo
-!              iloc = maxloc(q(i,j,l,:))
-!              iqmax = iloc(1)
-!              q(i,j,l,iqmax) = q(i,j,l,iqmax) + Smmr_old - Smmr_new
-!            enddo
-!           enddo
-!          enddo
 
           write(*,*)