Changeset 2425 for trunk

Timestamp:

Nov 2, 2020, 4:56:49 PM (4 years ago)

Author:

emillour

Message:

Mars GCM:
Adaptation for when using dust injection and/or dustscaling_mode==2. The
dust scenarios are then meant to just hold one value per sol deemed to
be valid for a given local time (t_scenario_sol in dust_param_mod). Thus
to avoid unwanted temporal interpolation in read_dust_scenario, the time
axis there should simply contain integers.
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/read_dust_scenario.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -3186 +3186 @@
 - 2D fields "dust_rad_adjust_prev" and "dust_rad_adjust_next" required to
   compute coefficient "dust_rad_adjust" need to be stored in (re)startfi files
+
+==02/11/2020 == EM
+Adaptation for when using dust injection and/or dustscaling_mode==2. The
+dust scenarios are then meant to just hold one value per sol deemed to
+be valid for a given local time (t_scenario_sol in dust_param_mod). Thus
+to avoid unwanted temporal interpolation in read_dust_scenario, the time
+axis there should simply contain integers.

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

@@ -6 +6 @@
 use geometry_mod, only: latitude, longitude ! in radians
 use datafile_mod, only: datadir
-use dust_param_mod, only: odpref
+use dust_param_mod, only: odpref, dustscaling_mode
+use planete_h, only: year_day
 implicit none
 
@@ -158 +159 @@
    ENDIF
 
+   ! Adapt time axis values if using dustinjection scheme and/or
+   ! dustscaling_mode == 2
+   IF ((dustinjection>0).or.(dustscaling_mode==2)) THEN
+     ! Sanity check: we expect to have one map per sol
+     if (timelen/=year_day) then
+       write(*,*)"read_dust_scenario error: timelen=",timelen
+       write(*,*)" whereas it was exepceted to be year_day=",year_day
+       call abort_physic("read_dust_scenario","timelen/=year_day",1)
+     endif
+     ! fill time() with intergers from 0 to timelen-1
+     write(*,*)"read_dust_scenario: adapting time axis to integer values"
+     do iloop=1,timelen
+       time(iloop)=iloop-1
+     enddo
+   ENDIF ! of IF ((dustinjection>0).or.(dustscaling_mode==2))
+
    ierr=nf90_inq_varid(nid,"latitude",nvarid)
    ierr=nf90_get_var(nid,nvarid,lat)