Changeset 3173 for trunk/LMDZ.COMMON/libf/evolution
- Timestamp:
- Jan 9, 2024, 10:55:18 AM (12 months ago)
- Location:
- trunk/LMDZ.COMMON/libf/evolution
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/changelog.txt
r3172 r3173 196 196 Fixing bug when rewritting the startfi.nc for the PCM: fluxgeo, read in the run_PEM.def, was not written in the startfi. 197 197 198 == 09/01/2023 == JBC 199 - Correction of script "launch_orb_1Dchained.sh" which read orbital parameters missing one interval of years in two. 200 - Addition of the Martian date in "info_PEM.txt" for post-processing. -
trunk/LMDZ.COMMON/libf/evolution/deftank/launch_orb_1Dchained.sh
r3149 r3173 85 85 yearlask=$(echo "$eyears_bp_ini + $i_myear*$convert_years" | bc -l) 86 86 found=false 87 while read -r y1 obl1 ecc1 lsp1 && read -r y2 obl2 ecc2 lsp2; do 88 y1=$(echo "$y1*1000" | bc -l) 89 y2=$(echo "$y2*1000" | bc -l) 87 read -r y1 obl1 ecc1 lsp1 < "$orb_data" 88 while read -r y2 obl2 ecc2 lsp2; do 89 if [ "$first_line" = true ]; then 90 first_line=false 91 continue # Skip the first line as it's already read outside the loop 92 fi 90 93 if [ "$(echo "$y1 >= $yearlask && $yearlask > $y2" | bc)" -eq 1 ]; then 91 94 found=true … … 107 110 break 108 111 fi 112 y1=$y2 113 obl1=$obl2 114 ecc1=$ecc2 115 lsp1=$lsp2 109 116 done < "$orb_data" 110 117 if [ "$found" = false ]; then -
trunk/LMDZ.COMMON/libf/evolution/info_PEM_mod.F90
r3149 r3173 17 17 !======================================================================= 18 18 19 use time_evol_mod, only: convert_years 19 use time_evol_mod, only: convert_years, year_bp_ini 20 20 21 21 implicit none … … 43 43 endif 44 44 open(1,file = 'info_PEM.txt',status = "old",position = "append",action = "write") 45 write(1,*) year_iter, i_myear, stopPEM 45 ! Martian date, Number of Martians years done by the PEM run, Number of Martians years done by the chainded simulation, Code of the stopping criterion 46 ! The conversion ratio from Planetary years to Earth years is given in the header of the file 47 write(1,*) year_bp_ini + i_myear, year_iter, i_myear, stopPEM 46 48 close(1) 47 49 else
Note: See TracChangeset
for help on using the changeset viewer.