Changeset 3860 for trunk


Ignore:
Timestamp:
Jul 18, 2025, 6:36:32 PM (3 days ago)
Author:
jbclement
Message:

PEM:

  • Correction of "visu_evol_layering.py" to compute Lsp variations.
  • Few cleanings to open files throughout the PEM.

JBC

Location:
trunk/LMDZ.COMMON/libf/evolution
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/evolution/changelog.txt

    r3859 r3860  
    752752== 18/07/2025 == JBC
    753753Setting up an automatic selection of the XIOS output file for the PEM according to the number of slopes.
     754
     755== 18/07/2025 == JBC
     756- Correction of "visu_evol_layering.py" to compute Lsp variations.
     757- Few cleanings to open files throughout the PEM.

  • trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90

    r3821 r3860  
    4343
    4444! #---------- Martian years parameters from launching script
    45 open(100,file = 'info_PEM.txt',status = 'old',form = 'formatted',iostat = ierr)
     45open(100,file = 'info_PEM.txt',status = 'old',form = 'formatted',action = 'read',iostat = ierr)
    4646if (ierr /= 0) then
    4747    write(*,*) 'Cannot find required file "info_PEM.txt"!'

  • trunk/LMDZ.COMMON/libf/evolution/deftank/visu_evol_layering.py

    r3851 r3860  
    969969
    970970    eccentricity  = (aphelion - perihelion) / (aphelion + perihelion)
    971     ls_perihelion = mars_ls(
    972         date_peri_day,
    973         date_peri_day,
    974         eccentricity,
    975         year_day
    976     )
     971    ls_perihelion = mars_ls(date_peri_day,0.,eccentricity,year_day)
    977972
    978973    return dates_yr, obliquity, eccentricity, ls_perihelion

  • trunk/LMDZ.COMMON/libf/evolution/lask_param_mod.F90

    r3571 r3860  
    2929
    3030nlask = 0
    31 open(1,file = 'obl_ecc_lsp.asc')
     31open(1,file = 'obl_ecc_lsp.asc',status = 'old',action = 'read')
    3232do
    3333    read(1,*,iostat = ierr)

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

    r3791 r3860  
    7575
    7676! We read the file
    77 open(73,file = 'obl_ecc_lsp.asc')
     77open(73,file = 'obl_ecc_lsp.asc',status = 'old',action = 'read')
    7878last_ilask = 0
    7979do ilask = 1,size(yearlask,1)

  • trunk/LMDZ.COMMON/libf/evolution/writediagpem.F90

    r3603 r3860  
    130130!$OMP MASTER
    131131    ! Open diagpem.def definition file if there is one:
    132         open(99,file="diagpem.def",status='old',form='formatted',iostat=ierr2)
     132        open(99,file="diagpem.def",status='old',form='formatted',action = 'read',iostat=ierr2)
    133133
    134134        if (ierr2 == 0) then
Note: See TracChangeset for help on using the changeset viewer.