Changeset 2634

Timestamp:

Mar 3, 2022, 10:17:13 PM (3 years ago)

Author:

abierjon

Message:

Mars GCM:
Changes on dust_rad_adjust for GCM6:

• put an upper limit on dust_rad_adjust to avoid skyrocketing values when there are strong diurnal variations of opacity
• move the condition to compute dust_rad_adjust only once per time step before the whole call to compute_dust_rad_adjust, instead of repetitive local checks
• some cleaning

AB

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/aeropacity_mod.F (modified) (2 diffs)
• libf/phymars/callradite_mod.F (modified) (1 diff)
• libf/phymars/dust_rad_adjust_mod.F90 (modified) (4 diffs)
• libf/phymars/dust_scaling_mod.F90 (modified) (3 diffs)
• libf/phymars/rocketduststorm_mod.F90 (modified) (1 diff)
• libf/phymars/topmons_mod.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -3630 +3630 @@
 startfi.nc in the output file start_archive.nc, and newstart reads the variable
 albedo in start_archive.nc before outputing it in restartfi.nc
+
+== 03/03/2022 == AB
+Changes on dust_rad_adjust for GCM6:
+- put an upper limit on dust_rad_adjust to avoid skyrocketing values when there are strong diurnal
+  variations of opacity
+- move the condition to compute dust_rad_adjust only once per time step before the whole call to
+  compute_dust_rad_adjust, instead of repetitive local checks
+- some cleaning

trunk/LMDZ.MARS/libf/phymars/aeropacity_mod.F

@@ -109 +109 @@
                              ! false to compute RT in the topdust
       REAL,INTENT(OUT) :: tauscaling(ngrid) ! Scaling factor for qdust and Ndust
-      REAL,INTENT(OUT) :: dust_rad_adjust(ngrid) ! Radiative adjustment 
+      REAL,INTENT(INOUT) :: dust_rad_adjust(ngrid) ! Radiative adjustment 
                           ! factor for dust
       REAL,INTENT(IN) :: totcloudfrac(ngrid) ! total water ice cloud fraction
@@ -597 +597 @@
 c       and once in the part of the mesh without the sub-grid mountain (nohmons=true)
         aerosol(1:ngrid,1:nlayer,iaer) = 0.
-        IF (nohmons) THEN  ! considering part of the mesh without top flows
+        IF (nohmons) THEN  ! considering part of the mesh without top dust flows
           aerosol(1:ngrid,1:nlayer,iaer)=1.e-25
         ELSE  ! part of the mesh with concentrated dust flows

trunk/LMDZ.MARS/libf/phymars/callradite_mod.F

@@ -177 +177 @@
       REAL,INTENT(INOUT) :: tauscaling(ngrid) ! Conversion factor for 
                                ! qdust and Ndust
-      REAL,INTENT(OUT) :: dust_rad_adjust(ngrid) ! Radiative adjustment 
+      REAL,INTENT(INOUT) :: dust_rad_adjust(ngrid) ! Radiative adjustment 
                           ! factor for dust
       REAL,INTENT(IN) :: albedo(ngrid,2),emis(ngrid)

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

@@ -23 +23 @@
   integer,intent(in) :: ngrid ! number of atmospheric columns
   integer,intent(in) :: nlayer ! number of atmospheric levels
-  real,intent(in) :: zday ! tim (in sols and fraction thereof)
+  real,intent(in) :: zday ! time (in sols and fraction thereof)
   real,intent(in) :: pplev(ngrid,nlayer+1) ! pressure (Pa) at layer boundaries
   real,intent(in) :: taudust(ngrid) ! visible dust columns opacity in the GCM
@@ -31 +31 @@
   real,allocatable,save :: local_time(:) ! LT at current physics time step
   real,allocatable,save :: local_time_prevdt(:) ! LT at previous physics time step
-  real :: zday_prevdt !value of zday at previous physics time step
+  real :: zday_prevdt ! value of zday at previous physics time step
   real,save :: zday_scenario ! to fetch dod values from the scenario
   real,save :: zday_scenario_next ! to fetch dod values from the scenario the next day
-  logical,save :: firstcall=.true.
+  logical,save :: firstcall=.true. ! is it the first call of the run?
   integer :: ig
 !  real,allocatable,save :: tau_pref_scenario(:)
   real,allocatable,save :: tau_pref_scenario_next(:)
   real :: weight ! interpolation weight
-  real,save :: zday_prev_call=-666. ! stored value of zday from previous call
   
-!$OMP THREADPRIVATE(  local_time,local_time_prevdt,zday_scenario,zday_scenario_next)
-!$OMP THREADPRIVATE(firstcall,tau_pref_scenario_next,zday_prev_call)
+!$OMP THREADPRIVATE(local_time,local_time_prevdt,zday_scenario,zday_scenario_next)
+!$OMP THREADPRIVATE(firstcall,tau_pref_scenario_next)
 
   
@@ -60 +59 @@
   endif ! of if firstcall
   
-  ! 1. Compute local times (in sol fraction), if not already done
-  if (zday/=zday_prev_call) then
-    local_time(1:ngrid)=modulo(1.+(zday-INT(zday)) + &
-                         (longitude_deg(1:ngrid)/15)/24,1.)
-    zday_prevdt=zday-dtphys/daysec
-    local_time_prevdt(1:ngrid)=modulo(1.+(zday_prevdt-INT(zday_prevdt)) + &
-                         (longitude_deg(1:ngrid)/15)/24,1.)
-    
-    zday_scenario=zday-modulo(zday,1.) ! integer value of the day: the scenario
-    ! opacity is assumed to be measured at 2pm but stored at nidnight
-    zday_scenario_next=zday_scenario+1
-  endif ! of if (zday/=zday_prev_call)
+  ! 1. Compute local times (in sol fraction)
+  local_time(1:ngrid)=modulo(1.+(zday-INT(zday)) + &
+                       (longitude_deg(1:ngrid)/15)/24,1.)
+  zday_prevdt=zday-dtphys/daysec
+  local_time_prevdt(1:ngrid)=modulo(1.+(zday_prevdt-INT(zday_prevdt)) + &
+                       (longitude_deg(1:ngrid)/15)/24,1.)
+
+  zday_scenario=zday-modulo(zday,1.) ! integer value of the day: the scenario
+  ! opacity is assumed to be measured at 2pm but stored at nidnight
+  zday_scenario_next=zday_scenario+1
   
   ! 2. Load dust opacities for zday_scenario and zday_scenario_next
-  !    if not already done
-  if (zday/=zday_prev_call) then
-!    call read_dust_scenario(ngrid,nlayer,zday_scenario,pplev,     &
-!                                         tau_pref_scenario)
-    call read_dust_scenario(ngrid,nlayer,zday_scenario_next,pplev, &
-                                         tau_pref_scenario_next)
-  endif ! of if (zday/=zday_prev_call)
+!  call read_dust_scenario(ngrid,nlayer,zday_scenario,pplev,     &
+!                                       tau_pref_scenario)
+  call read_dust_scenario(ngrid,nlayer,zday_scenario_next,pplev, &
+                                       tau_pref_scenario_next)
 
   ! 3. Update dust_rad_adjust_* for grid points which just reached 2pm
-  ! but only when this routine is called for the first time
-  ! during this time step
-  if (zday/=zday_prev_call) then
-   do ig=1,ngrid
-    if ((local_time(ig).ge.t_scenario_sol).and. &
-                 (local_time_prevdt(ig).lt.(t_scenario_sol))) then
-      ! store previous "next" as "prev" (NB we could also decide to recompute
-      ! it using the current taudust...)
-      dust_rad_adjust_prev(ig)=dust_rad_adjust_next(ig)
-      ! compute new target based on current dust opacity
-      dust_rad_adjust_next(ig)=tau_pref_scenario_next(ig)* &
-                               pplev(ig,1)/odpref/taudust(ig)
-    endif
-   enddo
-  endif ! of if (zday/=zday_prev_call)
+  do ig=1,ngrid
+   if ((local_time(ig).ge.t_scenario_sol).and. &
+                (local_time_prevdt(ig).lt.(t_scenario_sol))) then
+     ! store previous "next" as "prev" (NB we could also decide to recompute
+     ! it using the current taudust...)
+     dust_rad_adjust_prev(ig)=dust_rad_adjust_next(ig)
+     ! compute new target based on current dust opacity
+     dust_rad_adjust_next(ig)=tau_pref_scenario_next(ig)* &
+                              pplev(ig,1)/odpref/taudust(ig)
+                              
+     ! Upper limit for dust_rad_adjust
+     ! this is done to avoid skyrocketing taus given to the radiative transfer
+     ! when there is a huge diurnal variation of tau (especially in polar nights)
+      dust_rad_adjust_next(ig)=min(dust_rad_adjust_next(ig), 5.) !3.
+      
+   endif
+  enddo
   
   ! 4. Compute dust_rad_adjust using linear interpolation
@@ -116 +112 @@
                         (dust_rad_adjust_next(ig)-dust_rad_adjust_prev(ig))
   enddo! of do=ig=1,ngrid
-  
-  ! update zday_prev_call
-  zday_prev_call=zday
   
   end subroutine compute_dust_rad_adjust

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

@@ -25 +25 @@
                        ! opacity column at odpref reference pressure
     real,intent(out) :: tauscaling(ngrid) ! dust scaling factor
-    real,intent(out) :: dust_rad_adjust(ngrid) ! Radiative adjustment 
+    real,intent(inout) :: dust_rad_adjust(ngrid) ! Radiative adjustment 
                           ! factor for dust
     real,intent(inout) :: aerosol(ngrid,nlayer,naerkind) ! opacities
@@ -31 +31 @@
     integer :: ig, l , iaer
     real :: taudust(ngrid)
+    real,save :: zday_prev_call=-666. ! stored value of zday from previous call
+!$OMP THREADPRIVATE(zday_prev_call)
 
   ! 1. compute/set tauscaling
@@ -68 +70 @@
     elseif (dustscaling_mode==2) then
       ! GCM v6 style, compute dust_rad_adjust
-      call compute_dust_rad_adjust(ngrid,nlayer,zday,pplev, &
-                                   taudust,dust_rad_adjust)
+      ! only when this routine is called for the first time
+      ! during this time step
+      if (zday/=zday_prev_call) then
+        call compute_dust_rad_adjust(ngrid,nlayer,zday,pplev, &
+                                     taudust,dust_rad_adjust)
+       endif ! of if (zday/=zday_prev_call)
+                                     
+      ! update zday_prev_call
+      zday_prev_call=zday
     endif
 

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

@@ -79 +79 @@
       REAL, INTENT(IN) :: totstormfract(ngrid)
       REAL, INTENT(INOUT) :: tauscaling(ngrid)
-      REAL,INTENT(OUT) :: dust_rad_adjust(ngrid)
+      REAL,INTENT(INOUT) :: dust_rad_adjust(ngrid)
       
 !     subgrid scale water ice clouds

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

@@ -74 +74 @@
       INTEGER, INTENT(IN) :: igout
       REAL, INTENT(INOUT) :: tauscaling(ngrid)
-      REAL,INTENT(OUT) :: dust_rad_adjust(ngrid)
+      REAL,INTENT(INOUT) :: dust_rad_adjust(ngrid)
 !     input sub-grid scale rocket dust storm
       LOGICAL, INTENT(IN) :: clearatm