Changeset 3801

Timestamp:

Jun 6, 2025, 5:52:05 PM (3 weeks ago)

Author:

mmaurice

Message:

Generic PCM:

Add qsurf flag in modern traceur.def (used in 1D)

MM

Location:

trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Files:

• inichim_1D.F90 (modified) (1 diff)
• initracer_1D.F90 (modified) (7 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/inichim_1D.F90

@@ -120 +120 @@
          else
             write(*,*) ' Parameter value (default)     : qxf=', qxf(iq)
+         end if
+         if (index(tracline,'qsurf='    ) /= 0) then
+             read(tracline(index(tracline,'qsurf=')+len('qsurf='):),*) qsurf(iq)
+             write(*,*) ' Parameter value (traceur.def) : qsurf=', qsurf(iq)
+         else
+             write(*,*) ' Parameter value (default)     : qsurf=', qsurf(iq)
          end if
       end do

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/initracer_1D.F90

@@ -3 +3 @@
 !implicit none
 
-    subroutine initracer_1D(nlayer, nq, pq, qsurf, play)
+    subroutine initracer_1D(nlayer, nq, pq, pqsurf, play)
 
     use tracer_h, only: noms, mmol
@@ -42 +42 @@
     ! outputs :
 
-    real,intent(out)   :: pq(nlayer,nq)  ! tracers (kg/kg_of_air)
-    real,intent(out)   :: qsurf(nq)      ! surface values (kg/m2) of tracers
+    real,intent(inout)   :: pq(nlayer,nq)  ! tracers (kg/kg_of_air)
+    real,intent(inout)   :: pqsurf(nq)     ! surface values (kg/m2) of tracers
 
     ! local :
@@ -49 +49 @@
     real, allocatable  :: pf(:)          ! pressure in vmr profile files set in traceur.def
     real, allocatable  :: qf(:)          ! vmr      in vmr profile files set in traceur.def
-
+    real    :: qsurf(nq)                 ! surface tracer value
     real    :: mmean(nlayer)             ! mean molecular mass (g)
     real    :: pqx(nlayer,nq)            ! tracers (vmr)
@@ -112 +112 @@
             write(*,*) ' Parameter value (default)     : qxf=', qxf(iq)
         end if
+        if (index(tracline,'qsurf='    ) /= 0) then
+            read(tracline(index(tracline,'qsurf=')+len('qsurf='):),*) qsurf(iq)
+            write(*,*) ' Parameter value (traceur.def) : qsurf=', qsurf(iq)
+        else
+            write(*,*) ' Parameter value (default)     : qsurf=', qsurf(iq)
+        end if
     end do
 
@@ -125 +131 @@
                 if (pq(ilay,iq) /= 0.) then
                     print*, "initracer_1D: Tracer "//trim(noms(iq))//" seems to be already initialized"
-                    print*, "              and you are now trying set it a homogeneous value (flag qx)"
+                    print*, "              and you are now trying set it to a homogeneous value (flag qx)"
                     call abort
                 end if
@@ -134 +140 @@
                 if (pq(ilay,iq) /= 0.) then
                     print*, "initracer_1D: Tracer "//trim(noms(iq))//" seems to be already initialized"
-                    print*, "              and you are now trying set it a homogeneous value (flag qx)"
+                    print*, "              and you are now trying set it to an input porfile (flag qxf)"
                     call abort
                 end if
@@ -211 +217 @@
     end do
 
+    ! surface tracers
+    do iq=1,nq
+        if (pqsurf(iq) == 0.) then
+            pqsurf(iq) = qsurf(iq)
+        else
+            write(*,*) 'initracer_1D: you provided qsurf for '//trim(noms(iq))//' defined as background'
+            write(*,*) '              but it seems to have already been initialized: keep previous value'
+        end if
+    end do
     ! 4. Hard coding
+
     ! Do whatever you want here to specify pq and qsurf
     ! Or use #ModernTrac-v1 and add another option section 2.