Changeset 669 for trunk/LMDZ.MARS/libf/phymars/tabfi.F

Timestamp:

May 24, 2012, 5:02:41 PM (13 years ago)

Author:

tnavarro

Message:

Ice effective radius for output is weighted by ice total surface instead of ice total mass in physiq.F

• New distribution for permanent ice reservoirs in surfini.F

icelocationmode = 1 allows for realistic ice distribution, no matter what resolution.
icelocationmode = 2 is predefined 32x24 or 64x48 and should be USED BY DEFAULT. It currently overestimates ice.
icelocationmode = 3 is good for quick logical definitions (paleoclimates,stability studies,etc...)

• Possibility to change perihelion date in solar longitude as well as in sol in tabfi.F for newstarts
• Improved latent heat release when ground ice sublimates: now a implicit scheme for ground temperature

File:

• trunk/LMDZ.MARS/libf/phymars/tabfi.F (modified) (5 diffs)

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

@@ -67 +67 @@
       INTEGER Lmodif
       INTEGER lmax
-      REAL p_rad,p_omeg,p_g,p_mugaz,p_daysec,time
+      REAL p_rad,p_omeg,p_g,p_mugaz,p_daysec,time,peri_ls
 
 c Variables
@@ -289 +289 @@
       write(*,*) '(18)        obliquit : planet obliquity (deg)'
       write(*,*) '(17)      peri_day : perihelion date (sol since Ls=0)'
+      write(*,*) '(  )      peri_ls : perihelion date (Ls since Ls=0)'
       write(*,*) '(15)      periheli : min. sun-mars dist (Mkm)'
       write(*,*) '(16)      aphelie  : max. sun-mars dist (Mkm)'
@@ -429 +430 @@
         else if (trim(modif) .eq. 'peri_day') then
           write(*,*) 'current value:',peri_day
-          write(*,*) 'peri_day should be 485 on current Mars'
+          write(*,*) 'peri_day should be 485 sols on current Mars'
           write(*,*) 'enter new value:'
  116      read(*,*,iostat=ierr) peri_day
@@ -435 +436 @@
           write(*,*) 
           write(*,*) ' peri_day (new value):',peri_day
+          
+        else if (trim(modif) .eq. 'peri_ls') then
+          write(*,*) 'peri_ls value is not stored in start files,'
+          write(*,*) 'but it should be 251 degrees on current Mars'
+          write(*,*) '(peri_day should be 485 sols on current Mars)'
+          write(*,*) 'enter new value:'
+ 1160     read(*,*,iostat=ierr) peri_ls
+          if(ierr.ne.0) goto 1160
+          write(*,*) 
+          write(*,*) ' peri_ls asked :',peri_ls
+          call lsp2solp(peri_ls,peri_day)
+          write(*,*) ' peri_day (new value):',peri_day
+
 
         else if (trim(modif) .eq. 'periheli') then
@@ -529 +543 @@
 c-----------------------------------------------------------------------
       end
+
+
+
+      
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! gives sol at perihelion for ls at perihelion (for precession cycles)
+      subroutine lsp2solp(lsp,solp)
+
+      implicit none
+!  Arguments:
+      real lsp     ! Input: ls at perihelion
+      real solp    ! Output: sol at perihelion
+      
+!  Local:
+      double precision zx0 ! eccentric anomaly at Ls=0
+      double precision year_day 
+      double precision e_elips
+      double precision pi,degrad 
+      
+      parameter (year_day=668.6d0) ! number of sols in a martian year
+      parameter (e_elips=0.0934d0)  ! eccentricity of orbit
+      parameter (pi=3.14159265358979d0)
+      parameter (degrad=57.2957795130823d0)
+      
+      zx0 = -2.0*datan(dtan(0.5*lsp/degrad)
+     .          *dsqrt((1.-e_elips)/(1.+e_elips)))
+      if (zx0 .le. 0.) zx0 = zx0 + 2.*pi
+      
+      solp  = year_day*(1.-(zx0-e_elips*dsin(zx0))/(2.*pi))
+
+
+      end subroutine lsp2solp
+
+
+