Ignore:
Timestamp:
12/18/18 18:25:28 (5 years ago)
Author:
dumas
Message:

format of the parameter file created by the code, possibility to use it in reading

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SOURCES/initial-phy-2.f90

    r124 r224  
    4747  read(num_param,runpar) 
    4848 
    49  
    50   ! ouverture du fichier num_rep_42 
    51   !------------------------------------- 
    52 428 format(A)             ! formats pour les ecritures dans 42 
    53  
    54   ! filin='../Param/param'//runname 
     49!  filin='../Param/param'//runname 
    5550  filin='param'//runname 
    5651  open(num_rep_42,file=filin) 
     
    6156  write(num_rep_42,*)'!  sources dans  ',trim(dirsource) 
    6257  write(num_rep_42,*) 
    63   write(num_rep_42,428)'!___________________________________________________________'  
    64   write(num_rep_42,428) '&runpar                  ! nom du bloc  parametres du run' 
     58  write(num_rep_42,'(A)')'!___________________________________________________________'  
     59  write(num_rep_42,'(A)') '&runpar                  ! nom du bloc  parametres du run' 
    6560  write(num_rep_42,*) 
    6661  write(num_rep_42,*) 'runname      =  "',runname,'"' 
    6762  write(num_rep_42,*) 'icompteur    = ', icompteur                               
    68   write(num_rep_42,*) 'iout    = ', iout  
     63  write(num_rep_42,*) 'iout         = ', iout  
    6964  write(num_rep_42,*) 'reprcptr     = "',trim(reprcptr),'"' 
     65  write(num_rep_42,*) 'itracebug    = ', itracebug 
    7066  write(num_rep_42,*) 'num_tracebug  = ', num_tracebug 
    71   write(num_rep_42,*) 'comment_run   ="',trim(comment_run),'"' 
     67  write(num_rep_42,*) 'comment_run   = "',trim(comment_run),'"' 
    7268  write(num_rep_42,*)'/'                             
    73   write(num_rep_42,428) '! runname    : nom de l experience (8 caracteres)  ' 
    74   write(num_rep_42,428) '! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hw' 
    75   write(num_rep_42,428) '! icompteur  : 3-> T seulement' 
    76   write(num_rep_42,428) '! iout   : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise ' 
    77   write(num_rep_42,428) '! reprcptr   : nom du fichier' 
    78   write(num_rep_42,428) '! itracebug: 1-> ecriture de traces au debut des routines' 
    79   write(num_rep_42,428) '! num_tracebug   numero d unite ecriture traces' 
    80   write(num_rep_42,428) '! comment_run: commentaire court sur le run' 
     69  write(num_rep_42,'(A)') '! runname    : nom de l experience (8 caracteres)  ' 
     70  write(num_rep_42,'(A)') '! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hw' 
     71  write(num_rep_42,'(A)') '! icompteur  : 3-> T seulement' 
     72  write(num_rep_42,'(A)') '! iout   : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise ' 
     73  write(num_rep_42,'(A)') '! reprcptr   : nom du fichier restart' 
     74  write(num_rep_42,'(A)') '! itracebug: 1-> ecriture de traces au debut des routines' 
     75  write(num_rep_42,'(A)') '! num_tracebug   numero d unite ecriture traces' 
     76  write(num_rep_42,'(A)') '! comment_run: commentaire court sur le run' 
    8177 
    8278  write(num_rep_42,*) 
     
    118114 
    119115 
    120   write(num_rep_42,428)'!___________________________________________________________'   
    121   write(num_rep_42,428) '&grdline                 ! bloc grounding line' 
     116  write(num_rep_42,'(A)')'!___________________________________________________________'   
     117  write(num_rep_42,'(A)') '&grdline                 ! bloc grounding line' 
    122118  write(num_rep_42,*) 
    123119  write(num_rep_42,*) 'igrdline     = ',igrdline 
     
    125121  write(num_rep_42,*) 'ibmelt_inv   = ',ibmelt_inv 
    126122  write(num_rep_42,*)'/'                            
    127   write(num_rep_42,428)'! igrdline   : 1 ligne d echouage fixée, sinon 0' 
    128   write(num_rep_42,428)'! Schoof     : 0 pas de Schoof, 1 flux de Schoof' 
    129   write(num_rep_42,428)'! ibmelt_inv : 0 cas std, 1 inversion du bmelt (avec igrdline=1)' 
     123  write(num_rep_42,'(A)')'! igrdline   : 1 ligne d echouage fixée, sinon 0' 
     124  write(num_rep_42,'(A)')'! Schoof     : 0 pas de Schoof, 1 flux de Schoof' 
     125  write(num_rep_42,'(A)')'! ibmelt_inv : 0 cas std, 1 inversion du bmelt (avec igrdline=1)' 
    130126  write(num_rep_42,*) 
    131127 
     
    145141 
    146142! dmr&aurel ##  namelist/timesteps/dtmin,dtmax,dtt,testdiag,tbegin,tend 
    147   write(num_rep_42,428)'!___________________________________________________________'   
    148   write(num_rep_42,428) '&timesteps                 ! bloc time steps' 
     143  write(num_rep_42,'(A)')'!___________________________________________________________'   
     144  write(num_rep_42,'(A)') '&timesteps                 ! bloc time steps' 
    149145  write(num_rep_42,*) 
    150146  write(num_rep_42,*) 'tend      = ',tend 
     
    157153 
    158154  write(num_rep_42,*)'/'                            
    159   write(num_rep_42,428)'! tous les temps en annees. tbegin et tend : debut et fin du run ' 
    160   write(num_rep_42,428)'! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax'       
    161   write(num_rep_42,428)'! dtt : pas de temps long' 
    162   write(num_rep_42,428)'! testdiag, pour gérer le pas de temps dynamique dt' 
    163   write(num_rep_42,428)'! ordres de grandeur (a moduler selon dx) : ' 
    164   write(num_rep_42,428)'! 40 km dtmin=2.e-2, dtmax=1., dtt=5., tesdiag=0.02' 
     155  write(num_rep_42,'(A)')'! tous les temps en annees. tbegin et tend : debut et fin du run ' 
     156  write(num_rep_42,'(A)')'! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax'       
     157  write(num_rep_42,'(A)')'! dtt : pas de temps long' 
     158  write(num_rep_42,'(A)')'! testdiag, pour gérer le pas de temps dynamique dt' 
     159  write(num_rep_42,'(A)')'! ordres de grandeur (a moduler selon dx) : ' 
     160  write(num_rep_42,'(A)')'! 40 km dtmin=2.e-2, dtmax=1., dtt=5., tesdiag=0.02' 
    165161 
    166162  write(num_rep_42,*) 
Note: See TracChangeset for help on using the changeset viewer.