Changeset 2228 for trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

Timestamp:

Jan 28, 2020, 4:35:21 PM (5 years ago)

Author:

jnaar

Message:

[GCM MARS 1D] Adding a new flag 'paleomars' (default false) in 1d GCM in order to use
eccentricity and solar longitude of perihelion as input parameters.
Added variables : halfaxe, excentric, Lsperi. New routine call_dayperi.F from FF to
compute perhelion date from prescribed eccentricity & Lsperi.

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -122 +122 @@
       character(len=44) :: txt
 
+c   New flag to compute paleo orbital configurations + few variables JN
+      REAL halfaxe, excentric, Lsperi
+      Logical paleomars
+
 c=======================================================================
 
@@ -151 +155 @@
       periheli = 206.66          ! minimum sun-mars distance (Mkm) ~206.66
       aphelie = 249.22           ! maximum sun-mars distance (Mkm) ~249.22
+      halfaxe = 227.94           ! demi-grand axe de l'ellipse
       peri_day =  485.           ! perihelion date (sols since N. Spring)
       obliquit = 25.2            ! Obliquity (deg) ~25.2         
+      excentric = 0.0934         ! Eccentricity (0.0934)         
  
 c     Planetary Boundary Layer and Turbulence parameters 
@@ -464 +470 @@
 c Orbital parameters
 c ------------------
-      print *,'Min. distance Sun-Mars (Mkm)?'
-      call getin("periheli",periheli)
-      write(*,*) " periheli = ",periheli
-
-      print *,'Max. distance Sun-Mars (Mkm)?'
-      call getin("aphelie",aphelie)
-      write(*,*) " aphelie = ",aphelie
-
-      print *,'Day of perihelion?'
-      call getin("periday",peri_day)
-      write(*,*) " periday = ",peri_day
-
-      print *,'Obliquity?'
-      call getin("obliquit",obliquit)
-      write(*,*) " obliquit = ",obliquit
+      paleomars=.false. ! Default: no water ice reservoir
+      call getin("paleomars",paleomars)
+      if (paleomars==.true.) then
+        write(*,*) "paleomars=", paleomars
+        write(*,*) "Orbital parameters from callphys.def"
+        write(*,*) "Enter eccentricity & Lsperi"
+        print *, 'Martian eccentricity (0<e<1) ?'
+        call getin('excentric ',excentric)
+        write(*,*)"excentric =",excentric 
+        print *, 'Solar longitude of perihelion (0<Ls<360) ?'
+        call getin('Lsperi',Lsperi )
+        write(*,*)"Lsperi=",Lsperi 
+        Lsperi = Lsperi*pi/180.0 ! Put it in rad for peri_day
+        periheli = halfaxe*(1-excentric)
+        aphelie = halfaxe*(1+excentric)
+        call call_dayperi(Lsperi,excentric,peri_day,year_day)
+        write(*,*) "Corresponding orbital params for GCM"
+        write(*,*) " periheli = ",periheli
+        write(*,*) " aphelie = ",aphelie
+        write(*,*) "date of perihelion (sol)",peri_day
+      else
+        write(*,*) "paleomars=", paleomars
+        write(*,*) "Default present-day orbital parameters"
+        write(*,*) "Unless specified otherwise"
+        print *,'Min. distance Sun-Mars (Mkm)?'
+        call getin("periheli",periheli)
+        write(*,*) " periheli = ",periheli
+
+        print *,'Max. distance Sun-Mars (Mkm)?'
+        call getin("aphelie",aphelie)
+        write(*,*) " aphelie = ",aphelie
+
+        print *,'Day of perihelion?'
+        call getin("periday",peri_day)
+        write(*,*) " periday = ",peri_day
+
+        print *,'Obliquity?'
+        call getin("obliquit",obliquit)
+        write(*,*) " obliquit = ",obliquit
+      end if
  
 c  latitude/longitude