Ignore:
Timestamp:
Feb 20, 2025, 5:51:07 PM (4 months ago)
Author:
jbclement
Message:

Generic:
Small adjustments for the description of the 'controle' array in the "start.nc" file.
JBC

Location:
trunk/LMDZ.COMMON/libf
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d_common/dynredem.F90

    r3636 r3637  
    8282    start_file_type="planeto"
    8383    idecal = 4
    84     ldscrpt = 28
     84    ldscrpt = 29
    8585  else
    8686    if ( planet_type.eq."titan" ) then
    8787      ! Titan inherited Earth-like start files with idecal=5
    8888      write(lunout,*) trim(modname),' : Titan start file'
    89       ldscrpt = 32
    9089    else
    9190      write(lunout,*) trim(modname),' : Earth-like start file'
    92       ldscrpt = 32
    9391    endif
    9492    idecal = 5
     93    ldscrpt = 32
    9594  endif
    9695
    … …  
    134133  tab_cntrl(idecal+17)  = grossismx
    135134  tab_cntrl(idecal+18)  = grossismy
    136   dscrpt_tab_cntrl(idecal + 1:idecal + 18) = (/     &
    137       int2fmtstr(idecal + 1)//" Radius of the planet                 ", &
    138       int2fmtstr(idecal + 2)//" Rotation of the planet (rad/s)       ", &
    139       int2fmtstr(idecal + 3)//" Gravity (m/s2)                       ", &
    140       int2fmtstr(idecal + 4)//" Specific heat Cp (J.kg-1.K-1)        ", &
    141       int2fmtstr(idecal + 5)//" = r/Cp (=kappa)                      ", &
    142       int2fmtstr(idecal + 6)//" Lenght of a sol (s)                  ", &
    143       int2fmtstr(idecal + 7)//" Dynamical time step (s)              ", &
    144       int2fmtstr(idecal + 8)//" Total energy                         ", &
    145       int2fmtstr(idecal + 9)//" Total pressure                       ", &
     135  dscrpt_tab_cntrl(idecal + 1:idecal + 18) = (/                          &
     136      int2fmtstr(idecal + 1)// " Radius of the planet                 ", &
     137      int2fmtstr(idecal + 2)// " Rotation of the planet (rad/s)       ", &
     138      int2fmtstr(idecal + 3)// " Gravity (m/s2)                       ", &
     139      int2fmtstr(idecal + 4)// " Specific heat Cp (J.kg-1.K-1)        ", &
     140      int2fmtstr(idecal + 5)// " = r/Cp (=kappa)                      ", &
     141      int2fmtstr(idecal + 6)// " Lenght of a sol (s)                  ", &
     142      int2fmtstr(idecal + 7)// " Dynamical time step (s)              ", &
     143      int2fmtstr(idecal + 8)// " Total energy                         ", &
     144      int2fmtstr(idecal + 9)// " Total pressure                       ", &
    146145      int2fmtstr(idecal + 10)//" Total enstrophy                      ", &
    147146      int2fmtstr(idecal + 11)//" Total enthalpy                       ", &
    … …  
    161160    tab_cntrl(idecal+23) = taux
    162161    tab_cntrl(idecal+24) = tauy
    163     dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/  &
     162    dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/                       &
    164163      int2fmtstr(idecal + 19)//" 1.                                   ", &
    165164      int2fmtstr(idecal + 20)//" Extention (in longitude) of zoom     ", &
    … …  
    176175    tab_cntrl(idecal+24) = 0.
    177176    IF( ysinus )  tab_cntrl(idecal+22) = 1.
    178     dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/  &
     177    dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/                       &
    179178      int2fmtstr(idecal + 19)//" -                                    ", &
    180179      int2fmtstr(idecal + 20)//" Extention (in longitude) of zoom     ", &
    … …  
    199198    !tab_cntrl(29)=hour_ini
    200199    tab_cntrl(29)=0
     200    dscrpt_tab_cntrl(29) = "(29) 0                                    "
    201201  endif
    202202

  • trunk/LMDZ.COMMON/libf/dyn3dpar/dynredem_p.F90

    r3636 r3637  
    8888      ! Titan inherited Earth-like start files with idecal=5
    8989      write(lunout,*) trim(modname),' : Titan start file'
    90       ldscrpt = 32
    9190    else
    9291      write(lunout,*) trim(modname),' : Earth-like start file'
    93       ldscrpt = 32
    9492    endif
    9593    idecal = 5
     94    ldscrpt = 32
    9695  endif
    9796
    … …  
    135134  tab_cntrl(idecal+17)  = grossismx
    136135  tab_cntrl(idecal+18)  = grossismy
    137   dscrpt_tab_cntrl(idecal + 1:idecal + 18) = (/     &
    138       int2fmtstr(idecal + 1)//" Radius of the planet                 ", &
    139       int2fmtstr(idecal + 2)//" Rotation of the planet (rad/s)       ", &
    140       int2fmtstr(idecal + 3)//" Gravity (m/s2)                       ", &
    141       int2fmtstr(idecal + 4)//" Specific heat Cp (J.kg-1.K-1)        ", &
    142       int2fmtstr(idecal + 5)//" = r/Cp (=kappa)                      ", &
    143       int2fmtstr(idecal + 6)//" Lenght of a sol (s)                  ", &
    144       int2fmtstr(idecal + 7)//" Dynamical time step (s)              ", &
    145       int2fmtstr(idecal + 8)//" Total energy                         ", &
    146       int2fmtstr(idecal + 9)//" Total pressure                       ", &
     136  dscrpt_tab_cntrl(idecal + 1:idecal + 18) = (/                          &
     137      int2fmtstr(idecal + 1)// " Radius of the planet                 ", &
     138      int2fmtstr(idecal + 2)// " Rotation of the planet (rad/s)       ", &
     139      int2fmtstr(idecal + 3)// " Gravity (m/s2)                       ", &
     140      int2fmtstr(idecal + 4)// " Specific heat Cp (J.kg-1.K-1)        ", &
     141      int2fmtstr(idecal + 5)// " = r/Cp (=kappa)                      ", &
     142      int2fmtstr(idecal + 6)// " Lenght of a sol (s)                  ", &
     143      int2fmtstr(idecal + 7)// " Dynamical time step (s)              ", &
     144      int2fmtstr(idecal + 8)// " Total energy                         ", &
     145      int2fmtstr(idecal + 9)// " Total pressure                       ", &
    147146      int2fmtstr(idecal + 10)//" Total enstrophy                      ", &
    148147      int2fmtstr(idecal + 11)//" Total enthalpy                       ", &
    … …  
    162161    tab_cntrl(idecal+23) = taux
    163162    tab_cntrl(idecal+24) = tauy
    164     dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/  &
     163    dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/                       &
    165164      int2fmtstr(idecal + 19)//" 1.                                   ", &
    166165      int2fmtstr(idecal + 20)//" Extention (in longitude) of zoom     ", &
    … …  
    177176    tab_cntrl(idecal+24) = 0.
    178177    IF( ysinus )  tab_cntrl(idecal+22) = 1.
    179     dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/  &
     178    dscrpt_tab_cntrl(idecal + 19:idecal + 24) = (/                       &
    180179      int2fmtstr(idecal + 19)//" -                                    ", &
    181180      int2fmtstr(idecal + 20)//" Extention (in longitude) of zoom     ", &
    … …  
    200199!    tab_cntrl(29)=hour_ini
    201200    tab_cntrl(29)=0
     201    dscrpt_tab_cntrl(29) = "(29) 0                                    "
    202202  endif
    203203
Note: See TracChangeset for help on using the changeset viewer.