Changeset 2826 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Nov 23, 2022, 11:11:26 AM (2 years ago)

Author:

romain.vande

Message:

Mars GCM:
The variable co2ice is deleted. All the co2 ice on surface is now in qsurf(:,igcm_co2).
CO2 tracer is now mandatory. diagfi output is unchanged.
RV

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -9 +9 @@
 
 subroutine phyetat0 (fichnom,tab0,Lmodif,nsoil,ngrid,nlay,nq, &
-                     day_ini,time0,tsurf,tsoil,albedo,emis,q2,qsurf,co2ice, &
+                     day_ini,time0,tsurf,tsoil,albedo,emis,q2,qsurf, &
                      tauscaling,totcloudfrac,wstar,watercap)
 
@@ -63 +63 @@
   real,intent(out) :: q2(ngrid,nlay+1) ! 
   real,intent(out) :: qsurf(ngrid,nq) ! tracers on surface
-  real,intent(out) :: co2ice(ngrid) ! co2 ice cover
   real,intent(out) :: tauscaling(ngrid) ! dust conversion factor
   real,intent(out) :: totcloudfrac(ngrid) ! total cloud fraction
@@ -348 +347 @@
 endif ! if (startphy_file)
 
-! CO2 ice cover
-if (startphy_file) then
-   call get_field("co2ice",co2ice,found,indextime)
-   if (.not.found) then
-     call abort_physic(modname, &
-                "phyetat0: Failed loading <co2ice>",1)
-   endif
-else
-   co2ice(:)=0.
-endif !if (startphy_file)
-write(*,*) "phyetat0: CO2 ice cover <co2ice> range:", &
-            minval(co2ice), maxval(co2ice)
-
 ! Dust conversion factor
 if (startphy_file) then
@@ -561 +547 @@
 
     if (startphy_file) then
-       call get_field(txt,qsurf(:,iq),found,indextime)
-       if (.not.found) then
-         write(*,*) "phyetat0: Failed loading <",trim(txt),">"
-         write(*,*) "         ",trim(txt)," is set to zero"
-         qsurf(:,iq)=0.
-       endif
-    else
+      if (txt.eq."co2") then
+        ! We first check if co2ice exist in the startfi file (old way)
+        ! CO2 ice cover
+        call get_field("co2ice",qsurf(:,iq),found,indextime)
+        ! If not, we load the corresponding tracer in qsurf
+        if (.not.found) then
+          call get_field(txt,qsurf(:,iq),found,indextime)
+          if (.not.found) then
+            call abort_physic(modname, &
+                "phyetat0: Failed loading co2ice. there is neither the variable co2ice nor qsurf",1)
+          endif
+        endif
+      else ! (not the co2 tracer)
+        call get_field(txt,qsurf(:,iq),found,indextime)
+        if (.not.found) then
+          write(*,*) "phyetat0: Failed loading <",trim(txt),">"
+          write(*,*) "         ",trim(txt)," is set to zero"
+          qsurf(:,iq)=0.
+        endif
+      endif !endif co2
+    else !(not startphy_file)
       qsurf(:,iq)=0.
     endif ! of if (startphy_file)
@@ -573 +573 @@
                  minval(qsurf(:,iq)), maxval(qsurf(:,iq))
   enddo ! of do iq=1,nq
+
+    if (txt.eq."co2") then
+      ! We first check if co2ice exist in the startfi file (old way)
+      ! CO2 ice cover
+      if (startphy_file) then
+        call get_field("co2ice",qsurf(:,iq),found,indextime)
+      ! If not, we load the corresponding tracer in qsurf
+        if (.not.found) then
+          call get_field(txt,qsurf(:,iq),found,indextime)
+          if (.not.found) then
+            call abort_physic(modname, &
+                "phyetat0: Failed loading co2ice",1)
+          endif
+        endif
+      else
+       ! If we run without startfile, co2ice is set to 0
+        qsurf(:,iq)=0.
+      endif !if (startphy_file)
+        write(*,*) "phyetat0: CO2 ice cover <co2ice> range:", &
+            minval(qsurf(:,iq)), maxval(qsurf(:,iq))
+    endif
+
 endif ! of if (nq.ge.1)