Changeset 3040 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Sep 12, 2023, 11:41:22 AM (16 months ago)

Author:

jbclement

Message:

Mars PCM:
The variable 'timeperi' (defined in "planete_h.F90" and computed in "iniorbit.F") is renamed into 'lsperi' and thus slightly changed to be coherent to the solar longitude of perihelion in radian. It can now be used out of the box by other subroutines/programs like the PEM.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• orbit_param_criterion_mod.F90 (modified) (2 diffs)
• recomp_orb_param_mod.F90 (modified) (4 diffs)

trunk/LMDZ.COMMON/libf/evolution/orbit_param_criterion_mod.F90

@@ -22 +22 @@
 #endif
 #ifndef CPP_STD
-    use planete_h,   only: e_elips, obliquit, timeperi
+    use planete_h,   only: e_elips, obliquit, lsperi
 #else
-    use planete_mod, only: e_elips, obliquit, timeperi
+    use planete_mod, only: e_elips, obliquit, lsperi
 #endif
   use time_evol_mod,  only: year_bp_ini, var_obl, var_ecc, var_lsp, convert_years
@@ -57 +57 @@
 ! 0. Initializations
 ! **********************************************************************
-Year = year_bp_ini + i_myear        ! We compute the current year
-Lsp = modulo(-timeperi*180./pi,360.) ! We convert in degree
+Year = year_bp_ini + i_myear ! We compute the current year
+Lsp = lsperi*180./pi         ! We convert in degree
 
 call ini_lask_param_mod ! Allocation of variables

trunk/LMDZ.COMMON/libf/evolution/recomp_orb_param_mod.F90

@@ -15 +15 @@
 #ifndef CPP_STD      
     use comconst_mod, only: pi
-    use planete_h,    only: e_elips, obliquit, timeperi, periheli, aphelie, p_elips, peri_day, year_day
+    use planete_h,    only: e_elips, obliquit, lsperi, periheli, aphelie, p_elips, peri_day, year_day
 #else
-    use planete_mod,  only: e_elips, obliquit, timeperi, periastr, apoastr, p_elips, peri_day, year_day
+    use planete_mod,  only: e_elips, obliquit, lsperi, periastr, apoastr, p_elips, peri_day, year_day
     use comcstfi_mod, only: pi
 #endif
@@ -56 +56 @@
 #endif
 
-Year = year_bp_ini + i_myear        ! We compute the current year
-Lsp = modulo(-timeperi*180./pi,360.) ! We convert in degree
+Year = year_bp_ini + i_myear ! We compute the current year
+Lsp = lsperi*180./pi         ! We convert in degree
 
 write(*,*) 'recomp_orb_param, Old year =', Year - year_iter
@@ -100 +100 @@
 enddo
 halfaxe = 227.94
-timeperi = Lsp*pi/180.
+Lsp = Lsp*pi/180.
 periheli = halfaxe*(1 - e_elips)
 aphelie =  halfaxe*(1 + e_elips)
@@ -107 +107 @@
 
 #ifndef CPP_STD  
-    call call_dayperi(timeperi,e_elips,peri_day,year_day)
+    call call_dayperi(Lsp,e_elips,peri_day,year_day)
 #endif