Changeset 3714 for trunk/LMDZ.VENUS

Timestamp:

Apr 4, 2025, 4:13:19 PM (3 months ago)

Author:

slebonnois

Message:

SL: small adjustment in radlwsw.F

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• phytrac_emiss.F (modified) (9 diffs)
• radlwsw.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/phytrac_emiss.F ¶

@@ -35 +35 @@
       use dimphy
       USE geometry_mod, only: cell_area
-      USE chemparam_mod,only:M_tr
+      USE chemparam_mod,only:M_tr,type_tr
       USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
       IMPLICIT none
@@ -70 +70 @@
       real :: deltatr(klon,klev,nqtot)
 
-
-
-
-      integer,parameter :: nblat=3,nblon=3,nbaire=3,nbflux=2,maxcell=16 
+c several identical tracers emitted at different places, with different
+c fluxes, at the same time : 
+c same_tracer=.true. and various number of areas and fluxes
+c      logical,parameter :: same_tracer=.true.
+c      integer,parameter :: nblat=3,nblon=3,nbaire=3,nbflux=2,maxcell=16 
+c Otherwise, emission of one of the tracers only, in one area only
+      logical,parameter :: same_tracer=.false.
+      integer,parameter :: nblat=1,nblon=1,nbaire=1,nbflux=1
+      integer,save :: iq_co2,iq_n2
+
+      integer,parameter :: maxcell=10000
       integer,parameter :: Nheight=3 ! layer emission (150m)
       integer,save :: Ncell(nbaire) 
@@ -81 +88 @@
       integer,save :: ig_zone(nblat,nblon,nbaire,nbflux,maxcell) 
       integer,save :: numcell(nblat,nblon,nbaire,nbflux)
-        
 
       INTEGER i, k, it
@@ -102 +108 @@
 
         ALLOCATE(M_tr(nqtot))
-        M_tr(:)=44.                ! CO2
+        ALLOCATE(type_tr(nqtot))
+        do i=1,nqtot
+           if (tname(i)(1:3).eq.'co2') then   ! CO2
+               M_tr(i)=44.
+               iq_co2 =i
+               type_tr(i)=1
+           endif
+           if (tname(i)(1:2).eq.'n2')  then   ! N2
+               M_tr(i)=28.
+               iq_n2 =i
+               type_tr(i)=1
+           endif
+        enddo
 
 C=========================================================================
 c Caracteristiques des traceurs emis:
 C=========================================================================
+
 
 c nombre total de traceur
@@ -115 +134 @@
          endif 
                  
-
-
-
 c flux de CO2 (kg/s/m2)
-         flux_surface_co2(1) = 5.*10.**-9.
-         flux_surface_co2(2) = 5.*10.**-15.
+c 4.5e-7 correspond a 1% (en masse) d une couche de 2E6 Pa en 500 jV
+c donc avec 5E-7, on vide le N2 en 1000 jV
+c         flux_surface_co2(1) = 5.E-5
+         flux_surface_co2(1) = 0.
 
 c nombre de cellule pour le cote du carre d'aire 
-         Ncell(1)= 2
-         Ncell(2)= 3
-         Ncell(3)= 4
-      
-
+         Ncell(1)= 9999
 
 c localisation zone emission
-         lat_zone(1) = 08.
-         lat_zone(2) = -50.
-         lat_zone(3) = 35.
-         lon_zone(1) = -172.
-         lon_zone(2) = -20.
-         lon_zone(3) = 70.
-
+         lat_zone(1) = 0.
+         lon_zone(1) = 0.
  
          if ((nbp_lon*nbp_lat)==1) then ! running a 1D simulation
@@ -153 +162 @@
             do iaire=1,nbaire
 
-             if ((xlat(i).ge.lat_zone(ilat))
+            if (  
+     &       (Ncell(iaire).ne.9999)
+c emission on some areas only
+     &      .and.((xlat(i).ge.lat_zone(ilat))
      &      .and.((xlat(i)-Ncell(iaire)*deltalat)
      &      .lt.lat_zone(ilat))
      &      .and.(xlon(i).le.lon_zone(ilon))
      &      .and.((xlon(i)+Ncell(iaire)*deltalon)
-     &      .gt.lon_zone(ilon))) then
+     &      .gt.lon_zone(ilon))) ) then
 
               do iflux=1,nbflux
@@ -178 +190 @@
           tr_seri(:,:,:)=0.
 
-
            do i=1,klon
             do k=1,klev
-              tr_seri(i,k,:)=1.-28/43.44*max(min(0.035,
-     &           0.035*(1.-log(paprs(i,k)/6.e6)/log(9.e6/6.e6))),0.)
+c Profile for VeGa2
+c              tr_seri(i,k,iq_n2) = M_tr(i_n2)/43.44*max(min(0.035,
+c     &           0.035*(1.-log(paprs(i,k)/6.e6)/log(9.e6/6.e6))),0.)
+c              tr_seri(i,k,iq_co2)= 1. - tr_seri(i,k,iq_n2)
+
+c Uniform initialization, yN2=3.5%, regular mean mol mass = 43.44 g/mol
+              tr_seri(i,k,iq_n2) = 0.035*M_tr(iq_n2)/43.44
+              tr_seri(i,k,iq_co2)= 1. - tr_seri(i,k,iq_n2)
             end do
            end do
@@ -202 +219 @@
        flux(:,:)=0.
 
+        if (same_tracer) then
 c emet les traceurs qui sont presents sur la grille
-      do ilat  = 1,nblat
-       do ilon  = 1,nblon
+        do ilat  = 1,nblat
+        do ilon  = 1,nblon
         do iaire = 1,nbaire
-         do iflux = 1,nbflux
+        do iflux = 1,nbflux
      
               it=min( (ilat-1)*nblon*nbflux*nbaire+(iaire-1)*nbflux
      &         +(ilon-1)*nbaire*nbflux+iflux , nqtot )   
  
-
-c injection dans une seule cellule:
-c source en kg/kg/s
-c            deltatr(i,Nheight(iz),it) = so2_quantity/(86400.*Nemiss) ! kg/s
-c     $ *RG/( area_emiss(ilat,ilon)
-c     $      *(paprs(i,Nheight(iz))-paprs(i,Nheight(iz)+1)) )    ! /kg (masse cellule)
-     
-c            tr_seri(i,Nheight(iz),it) = tr_seri(i,Nheight(iz),it)
-c     $      + deltatr(i,Nheight(iz),it)*pdtphys
-
-c injection dans toute la colonne (a faire):
+         if (Ncell(iaire).ne.9999) then
+c column injection
           do ipos=1,maxcell
  
@@ -231 +240 @@
             do k=1,Nheight
              deltatr(i,k,it) = flux_surface_co2(iflux) ! kg/s/m2
-     $         *RG/(paprs(i,1)-paprs(i,Nheight+1))    ! /kg (masse colonne)
-     
-               tr_seri(i,k,it) = tr_seri(i,k,it)+deltatr(i,k,it)*pdtphys
+     $         *RG/(paprs(i,1)-paprs(i,Nheight+1))    ! / (kg/m2) (masse colonne)
+     
+             tr_seri(i,k,it) = tr_seri(i,k,it)+deltatr(i,k,it)*pdtphys
             end do
 
            end if 
           end do
- 
-        end do
-        end do
-        end do
-        end do
-
+         endif
+        end do
+        end do
+        end do
+        end do
+
+        else  ! same_tracer=.false.
+
+c column injection !! with constant mass !!
+          do i=1,nlon
+ 
+             flux(i,iq_co2)=flux_surface_co2(1)  
+              
+            do k=1,Nheight
+          deltatr(i,k,iq_co2) = min((1.-tr_seri(i,k,iq_co2))/pdtphys,
+     $                                flux(i,iq_co2)) ! kg/s/m2
+     $         *RG/(paprs(i,1)-paprs(i,Nheight+1))    ! / (kg/m2) (masse colonne)
+             deltatr(i,k,iq_n2)  = -deltatr(i,k,iq_co2)
+             tr_seri(i,k,iq_co2) = tr_seri(i,k,iq_co2)
+     $                            +deltatr(i,k,iq_co2)*pdtphys
+             tr_seri(i,k,iq_n2)  = tr_seri(i,k,iq_n2)
+     $                            +deltatr(i,k,iq_n2)*pdtphys
+            end do
+
+          end do
+
+        endif ! same_tracer
+ 
        
 c======================================================================
 c======================================================================
 
-   
-      
-
-     
 #ifdef CPP_XIOS      
-       do it=1,nqtot
-       CALL  send_xios_field("flux_"//tname(it),
-     &                     flux(:,it)) 
-      end do 
+!      do it=1,nqtot
+!      CALL  send_xios_field("flux_"//tname(it),
+!    &                     flux(:,it)) 
+!     end do 
 #endif
-
-
 
       RETURN

trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F ¶

@@ -62 +62 @@
       INTEGER k, kk, i, j, band
       integer,save :: i_sw
-      integer,parameter :: subloop=100
+      integer,save :: subloop
 
       REAL   PPB(klev+1)
@@ -101 +101 @@
 
       if (firstcall) then
+
+c s2: solar only 240 times per Vd if diurnal cycle
+        if (cycle_diurne) then
+            subloop = nbapp_rad/240
+        else
+            subloop = nbapp_rad  ! no diurnal cycle, once per Vd is enough
+        endif
 
 c ---------- ksive --------------