Changeset 2689 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
Jun 18, 2022, 7:50:16 PM (2 years ago)
Author:
emillour
Message:

Mars GCM:
Enable reading MY35 E107 forcing.
EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/aeronomars/param_read_e107.F

    r2610 r2689  
    252252
    253253      !Tabulated values for E10.7
    254       if(solvaryear.eq.23) then
    255          filename="e107_MY23.dat"
    256       else if(solvaryear.eq.24) then
    257          filename="e107_MY24.dat"
    258       else if(solvaryear.eq.25) then
    259          filename="e107_MY25.dat"
    260       else if (solvaryear.eq.26) then
    261          filename="e107_MY26.dat"
    262       else if (solvaryear.eq.27) then
    263          filename="e107_MY27.dat"
    264       else if (solvaryear.eq.28) then
    265          filename="e107_MY28.dat"
    266       else if (solvaryear.eq.29) then
    267          filename="e107_MY29.dat"
    268       else if(solvaryear.eq.30) then
    269          filename="e107_MY30.dat"
    270       else if(solvaryear.eq.31) then
    271          filename="e107_MY31.dat"
    272       else if(solvaryear.eq.32) then
    273          filename="e107_MY32.dat"
    274       else if(solvaryear.eq.33) then
    275          filename="e107_MY33.dat"
    276       else if(solvaryear.eq.34) then
    277          filename="e107_MY34.dat"
     254      if((solvaryear.ge.23).and.(solvaryear.le.35)) then
     255        write(filename,fmt="(a,i2,a)")"e107_MY",solvaryear,".dat"
    278256      else
    279257         write(*,*)"param_read_e107: "
    280258         write(*,*)"bad value for solvaryear in callphys.def"
    281          write(*,*)"solvaryear must be between 24 and 33"
     259         write(*,*)"solvaryear must be between 24 and 35"
    282260         stop
    283261      endif
Note: See TracChangeset for help on using the changeset viewer.