Ignore:
Timestamp:
Sep 21, 2011, 3:57:48 PM (13 years ago)
Author:
jghattas
Message:
  • Added suffix _mpi_rank to lmdz.out text file. Each processus now write into seperate file but only if lunout/=6 (lunout is set in run.def).
  • Change some print* into write(lunout,*)
  • phytrac.F90 : always include ini_histrac and write_histrac. The file histrac.nc is written if ecrit_tra> 0 (set in physiq.def). Change default value of ecrit_tra into 0.
Location:
LMDZ5/trunk/libf/phylmd
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/concvl.F

    r1574 r1575  
    238238cc$$$         close (56)
    239239c
    240          print*, 'supcrit1, supcrit2' ,supcrit1, supcrit2
     240         IF (prt_level .ge. 10)
     241     &       WRITE(lunout,*) 'supcrit1, supcrit2' ,supcrit1, supcrit2
    241242C
    242243C===========================================================================
    … …  
    370371      endif 
    371372C------------------------------------------------------------------
    372       print *,' cva_driver -> cbmf,plcl,plfc,wbeff ',
    373      .          cbmf(1),plcl(1),plfc(1),wbeff(1)
     373      IF (prt_level .ge. 10)
     374     .   WRITE(lunout,*) ' cva_driver -> cbmf,plcl,plfc,wbeff ',
     375     .                     cbmf(1),plcl(1),plfc(1),wbeff(1)
    374376
    375377      DO i = 1,klon

  • LMDZ5/trunk/libf/phylmd/conf_phys.F90

    r1539 r1575  
    3737
    3838 include "thermcell.h"
     39 include "iniprint.h"
    3940
    4041!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
    … …  
    101102
    102103! Local
    103   integer              :: numout = 6
    104104  real                 :: zzz
    105105
    … …  
    199199  call getin('OCEAN', ocean_omp)
    200200  IF (ocean_omp /= 'yyyyyy') THEN
    201      WRITE(numout,*)'ERROR!! Old variable name OCEAN used in parmeter file.'
    202      WRITE(numout,*)'Variable OCEAN has been replaced by the variable type_ocean.'
    203      WRITE(numout,*)'You have to update your parameter file physiq.def to succed running'
     201     WRITE(lunout,*)'ERROR!! Old variable name OCEAN used in parmeter file.'
     202     WRITE(lunout,*)'Variable OCEAN has been replaced by the variable type_ocean.'
     203     WRITE(lunout,*)'You have to update your parameter file physiq.def to succed running'
    204204     CALL abort_gcm('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1)
    205205  END IF
    … …  
    13511351!Config Help =
    13521352!
    1353   ecrit_tra_omp = 30.
     1353  ecrit_tra_omp = 0.
    13541354  call getin('ecrit_tra',ecrit_tra_omp)
    13551355!
    … …  
    16961696! Test of coherence between type_ocean and version_ocean
    16971697    IF (type_ocean=='couple' .AND. (version_ocean/='opa8' .AND. version_ocean/='nemo') ) THEN
    1698        WRITE(numout,*)' ERROR version_ocean=',version_ocean,' not valid in coupled configuration'
     1698       WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid in coupled configuration'
    16991699       CALL abort_gcm('conf_phys','version_ocean not valid',1)
    17001700    END IF
    … …  
    17031703       version_ocean='sicOBS'
    17041704    ELSE IF (type_ocean=='slab' .AND. version_ocean/='sicOBS') THEN
    1705        WRITE(numout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean'
     1705       WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean'
    17061706       CALL abort_gcm('conf_phys','version_ocean not valid',1)
    17071707    END IF
    … …  
    17171717!$OMP MASTER
    17181718
    1719   write(numout,*)' ##############################################'
    1720   write(numout,*)' Configuration des parametres de la physique: '
    1721   write(numout,*)' Type ocean = ', type_ocean
    1722   write(numout,*)' Version ocean = ', version_ocean
    1723   write(numout,*)' Config veget = ', ok_veget
    1724   write(numout,*)' Sortie journaliere = ', ok_journe
    1725   write(numout,*)' Sortie haute frequence = ', ok_hf
    1726   write(numout,*)' Sortie mensuelle = ', ok_mensuel
    1727   write(numout,*)' Sortie instantanee = ', ok_instan
    1728   write(numout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
    1729   write(numout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
    1730   write(numout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
    1731   write(numout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
    1732   write(numout,*)' Excentricite = ',R_ecc
    1733   write(numout,*)' Equinoxe = ',R_peri
    1734   write(numout,*)' Inclinaison =',R_incl
    1735   write(numout,*)' Constante solaire =',solaire
    1736   write(numout,*)' co2_ppm =',co2_ppm
    1737   write(numout,*)' RCO2_act = ',RCO2_act
    1738   write(numout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
    1739   write(numout,*)' N2O_ppb =',N2O_ppb,' RN2O_act=  ',RN2O_act
    1740   write(numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act=  ',RCFC11_act
    1741   write(numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act=  ',RCFC12_act
    1742   write(numout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
    1743   write(numout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
    1744   write(numout,*)' RCFC12_per = ',RCFC12_per
    1745   write(numout,*)' cvl_corr=', cvl_corr
    1746   write(numout,*)'ok_lic_melt=', ok_lic_melt
    1747   write(numout,*)'cycle_diurne=',cycle_diurne
    1748   write(numout,*)'soil_model=',soil_model
    1749   write(numout,*)'new_oliq=',new_oliq
    1750   write(numout,*)'ok_orodr=',ok_orodr
    1751   write(numout,*)'ok_orolf=',ok_orolf
    1752   write(numout,*)'ok_limitvrai=',ok_limitvrai
    1753   write(numout,*)'nbapp_rad=',nbapp_rad
    1754   write(numout,*)'iflag_con=',iflag_con
    1755   write(numout,*)' epmax = ', epmax
    1756   write(numout,*)' ok_adj_ema = ', ok_adj_ema
    1757   write(numout,*)' iflag_clw = ', iflag_clw
    1758   write(numout,*)' cld_lc_lsc = ', cld_lc_lsc
    1759   write(numout,*)' cld_lc_con = ', cld_lc_con
    1760   write(numout,*)' cld_tau_lsc = ', cld_tau_lsc
    1761   write(numout,*)' cld_tau_con = ', cld_tau_con
    1762   write(numout,*)' ffallv_lsc = ', ffallv_lsc
    1763   write(numout,*)' ffallv_con = ', ffallv_con
    1764   write(numout,*)' coef_eva = ', coef_eva
    1765   write(numout,*)' reevap_ice = ', reevap_ice
    1766   write(numout,*)' iflag_pdf = ', iflag_pdf
    1767   write(numout,*)' iflag_cldcon = ', iflag_cldcon
    1768   write(numout,*)' iflag_radia = ', iflag_radia
    1769   write(numout,*)' iflag_rrtm = ', iflag_rrtm
    1770   write(numout,*)' iflag_ratqs = ', iflag_ratqs
    1771   write(numout,*)' seuil_inversion = ', seuil_inversion
    1772   write(numout,*)' fact_cldcon = ', fact_cldcon
    1773   write(numout,*)' facttemps = ', facttemps
    1774   write(numout,*)' ok_newmicro = ',ok_newmicro
    1775   write(numout,*)' ratqsbas = ',ratqsbas
    1776   write(numout,*)' ratqshaut = ',ratqshaut
    1777   write(numout,*)' tau_ratqs = ',tau_ratqs
    1778   write(numout,*)' top_height = ',top_height
    1779   write(numout,*)' rad_froid = ',rad_froid
    1780   write(numout,*)' rad_chau1 = ',rad_chau1
    1781   write(numout,*)' rad_chau2 = ',rad_chau2
    1782   write(numout,*)' t_glace_min = ',t_glace_min
    1783   write(numout,*)' t_glace_max = ',t_glace_max
    1784   write(numout,*)' rei_min = ',rei_min
    1785   write(numout,*)' rei_max = ',rei_max
    1786   write(numout,*)' overlap = ',overlap
    1787   write(numout,*)' cdmmax = ',cdmmax
    1788   write(numout,*)' cdhmax = ',cdhmax
    1789   write(numout,*)' ksta = ',ksta
    1790   write(numout,*)' ksta_ter = ',ksta_ter
    1791   write(numout,*)' ok_kzmin = ',ok_kzmin
    1792   write(numout,*)' fmagic = ',fmagic
    1793   write(numout,*)' pmagic = ',pmagic
    1794   write(numout,*)' ok_ade = ',ok_ade
    1795   write(numout,*)' ok_aie = ',ok_aie
    1796   write(numout,*)' aerosol_couple = ', aerosol_couple
    1797   write(numout,*)' flag_aerosol = ', flag_aerosol
    1798   write(numout,*)' new_aod = ', new_aod
    1799   write(numout,*)' aer_type = ',aer_type
    1800   write(numout,*)' bl95_b0 = ',bl95_b0
    1801   write(numout,*)' bl95_b1 = ',bl95_b1
    1802   write(numout,*)' lev_histhf = ',lev_histhf
    1803   write(numout,*)' lev_histday = ',lev_histday
    1804   write(numout,*)' lev_histmth = ',lev_histmth
    1805   write(numout,*)' lev_histins = ',lev_histins
    1806   write(numout,*)' lev_histLES = ',lev_histLES
    1807   write(numout,*)' lev_histdayNMC = ',lev_histdayNMC
    1808   write(numout,*)' ok_histNMC = ',ok_histNMC
    1809   write(numout,*)' freq_outNMC = ',freq_outNMC
    1810   write(numout,*)' freq_calNMC = ',freq_calNMC
    1811   write(numout,*)' iflag_pbl = ', iflag_pbl
    1812   write(numout,*)' iflag_thermals = ', iflag_thermals
    1813   write(numout,*)' iflag_thermals_ed = ', iflag_thermals_ed
    1814   write(numout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux
    1815   write(numout,*)' iflag_clos = ', iflag_clos
    1816   write(numout,*)' type_run = ',type_run
    1817   write(numout,*)' ok_isccp = ',ok_isccp
    1818   write(numout,*)' ok_cosp = ',ok_cosp
    1819   write(numout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP
    1820   write(numout,*)' ok_journeCOSP = ',ok_journeCOSP
    1821   write(numout,*)' ok_hfCOSP =',ok_hfCOSP
    1822   write(numout,*)' solarlong0 = ', solarlong0
    1823   write(numout,*)' qsol0 = ', qsol0
    1824   write(numout,*)' inertie_sol = ', inertie_sol
    1825   write(numout,*)' inertie_ice = ', inertie_ice
    1826   write(numout,*)' inertie_sno = ', inertie_sno
    1827   write(numout,*)' f_cdrag_ter = ',f_cdrag_ter
    1828   write(numout,*)' f_cdrag_oce = ',f_cdrag_oce
    1829   write(numout,*)' f_rugoro = ',f_rugoro
    1830   write(numout,*)' supcrit1 = ', supcrit1
    1831   write(numout,*)' supcrit2 = ', supcrit2
    1832   write(numout,*)' iflag_mix = ', iflag_mix
    1833   write(numout,*)' scut = ', scut
    1834   write(numout,*)' qqa1 = ', qqa1
    1835   write(numout,*)' qqa2 = ', qqa2
    1836   write(numout,*)' gammas = ', gammas
    1837   write(numout,*)' Fmax = ', Fmax
    1838   write(numout,*)' alphas = ', alphas
    1839   write(numout,*)' iflag_wake = ', iflag_wake
    1840   write(numout,*)' alp_offset = ', alp_offset
    1841 
    1842   write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
     1719  write(lunout,*)' ##############################################'
     1720  write(lunout,*)' Configuration des parametres de la physique: '
     1721  write(lunout,*)' Type ocean = ', type_ocean
     1722  write(lunout,*)' Version ocean = ', version_ocean
     1723  write(lunout,*)' Config veget = ', ok_veget
     1724  write(lunout,*)' Sortie journaliere = ', ok_journe
     1725  write(lunout,*)' Sortie haute frequence = ', ok_hf
     1726  write(lunout,*)' Sortie mensuelle = ', ok_mensuel
     1727  write(lunout,*)' Sortie instantanee = ', ok_instan
     1728  write(lunout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
     1729  write(lunout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
     1730  write(lunout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
     1731  write(lunout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
     1732  write(lunout,*)' Excentricite = ',R_ecc
     1733  write(lunout,*)' Equinoxe = ',R_peri
     1734  write(lunout,*)' Inclinaison =',R_incl
     1735  write(lunout,*)' Constante solaire =',solaire
     1736  write(lunout,*)' co2_ppm =',co2_ppm
     1737  write(lunout,*)' RCO2_act = ',RCO2_act
     1738  write(lunout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
     1739  write(lunout,*)' N2O_ppb =',N2O_ppb,' RN2O_act=  ',RN2O_act
     1740  write(lunout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act=  ',RCFC11_act
     1741  write(lunout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act=  ',RCFC12_act
     1742  write(lunout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
     1743  write(lunout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
     1744  write(lunout,*)' RCFC12_per = ',RCFC12_per
     1745  write(lunout,*)' cvl_corr=', cvl_corr
     1746  write(lunout,*)'ok_lic_melt=', ok_lic_melt
     1747  write(lunout,*)'cycle_diurne=',cycle_diurne
     1748  write(lunout,*)'soil_model=',soil_model
     1749  write(lunout,*)'new_oliq=',new_oliq
     1750  write(lunout,*)'ok_orodr=',ok_orodr
     1751  write(lunout,*)'ok_orolf=',ok_orolf
     1752  write(lunout,*)'ok_limitvrai=',ok_limitvrai
     1753  write(lunout,*)'nbapp_rad=',nbapp_rad
     1754  write(lunout,*)'iflag_con=',iflag_con
     1755  write(lunout,*)' epmax = ', epmax
     1756  write(lunout,*)' ok_adj_ema = ', ok_adj_ema
     1757  write(lunout,*)' iflag_clw = ', iflag_clw
     1758  write(lunout,*)' cld_lc_lsc = ', cld_lc_lsc
     1759  write(lunout,*)' cld_lc_con = ', cld_lc_con
     1760  write(lunout,*)' cld_tau_lsc = ', cld_tau_lsc
     1761  write(lunout,*)' cld_tau_con = ', cld_tau_con
     1762  write(lunout,*)' ffallv_lsc = ', ffallv_lsc
     1763  write(lunout,*)' ffallv_con = ', ffallv_con
     1764  write(lunout,*)' coef_eva = ', coef_eva
     1765  write(lunout,*)' reevap_ice = ', reevap_ice
     1766  write(lunout,*)' iflag_pdf = ', iflag_pdf
     1767  write(lunout,*)' iflag_cldcon = ', iflag_cldcon
     1768  write(lunout,*)' iflag_radia = ', iflag_radia
     1769  write(lunout,*)' iflag_rrtm = ', iflag_rrtm
     1770  write(lunout,*)' iflag_ratqs = ', iflag_ratqs
     1771  write(lunout,*)' seuil_inversion = ', seuil_inversion
     1772  write(lunout,*)' fact_cldcon = ', fact_cldcon
     1773  write(lunout,*)' facttemps = ', facttemps
     1774  write(lunout,*)' ok_newmicro = ',ok_newmicro
     1775  write(lunout,*)' ratqsbas = ',ratqsbas
     1776  write(lunout,*)' ratqshaut = ',ratqshaut
     1777  write(lunout,*)' tau_ratqs = ',tau_ratqs
     1778  write(lunout,*)' top_height = ',top_height
     1779  write(lunout,*)' rad_froid = ',rad_froid
     1780  write(lunout,*)' rad_chau1 = ',rad_chau1
     1781  write(lunout,*)' rad_chau2 = ',rad_chau2
     1782  write(lunout,*)' t_glace_min = ',t_glace_min
     1783  write(lunout,*)' t_glace_max = ',t_glace_max
     1784  write(lunout,*)' rei_min = ',rei_min
     1785  write(lunout,*)' rei_max = ',rei_max
     1786  write(lunout,*)' overlap = ',overlap
     1787  write(lunout,*)' cdmmax = ',cdmmax
     1788  write(lunout,*)' cdhmax = ',cdhmax
     1789  write(lunout,*)' ksta = ',ksta
     1790  write(lunout,*)' ksta_ter = ',ksta_ter
     1791  write(lunout,*)' ok_kzmin = ',ok_kzmin
     1792  write(lunout,*)' fmagic = ',fmagic
     1793  write(lunout,*)' pmagic = ',pmagic
     1794  write(lunout,*)' ok_ade = ',ok_ade
     1795  write(lunout,*)' ok_aie = ',ok_aie
     1796  write(lunout,*)' aerosol_couple = ', aerosol_couple
     1797  write(lunout,*)' flag_aerosol = ', flag_aerosol
     1798  write(lunout,*)' new_aod = ', new_aod
     1799  write(lunout,*)' aer_type = ',aer_type
     1800  write(lunout,*)' bl95_b0 = ',bl95_b0
     1801  write(lunout,*)' bl95_b1 = ',bl95_b1
     1802  write(lunout,*)' lev_histhf = ',lev_histhf
     1803  write(lunout,*)' lev_histday = ',lev_histday
     1804  write(lunout,*)' lev_histmth = ',lev_histmth
     1805  write(lunout,*)' lev_histins = ',lev_histins
     1806  write(lunout,*)' lev_histLES = ',lev_histLES
     1807  write(lunout,*)' lev_histdayNMC = ',lev_histdayNMC
     1808  write(lunout,*)' ok_histNMC = ',ok_histNMC
     1809  write(lunout,*)' freq_outNMC = ',freq_outNMC
     1810  write(lunout,*)' freq_calNMC = ',freq_calNMC
     1811  write(lunout,*)' iflag_pbl = ', iflag_pbl
     1812  write(lunout,*)' iflag_thermals = ', iflag_thermals
     1813  write(lunout,*)' iflag_thermals_ed = ', iflag_thermals_ed
     1814  write(lunout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux
     1815  write(lunout,*)' iflag_clos = ', iflag_clos
     1816  write(lunout,*)' type_run = ',type_run
     1817  write(lunout,*)' ok_isccp = ',ok_isccp
     1818  write(lunout,*)' ok_cosp = ',ok_cosp
     1819  write(lunout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP
     1820  write(lunout,*)' ok_journeCOSP = ',ok_journeCOSP
     1821  write(lunout,*)' ok_hfCOSP =',ok_hfCOSP
     1822  write(lunout,*)' solarlong0 = ', solarlong0
     1823  write(lunout,*)' qsol0 = ', qsol0
     1824  write(lunout,*)' inertie_sol = ', inertie_sol
     1825  write(lunout,*)' inertie_ice = ', inertie_ice
     1826  write(lunout,*)' inertie_sno = ', inertie_sno
     1827  write(lunout,*)' f_cdrag_ter = ',f_cdrag_ter
     1828  write(lunout,*)' f_cdrag_oce = ',f_cdrag_oce
     1829  write(lunout,*)' f_rugoro = ',f_rugoro
     1830  write(lunout,*)' supcrit1 = ', supcrit1
     1831  write(lunout,*)' supcrit2 = ', supcrit2
     1832  write(lunout,*)' iflag_mix = ', iflag_mix
     1833  write(lunout,*)' scut = ', scut
     1834  write(lunout,*)' qqa1 = ', qqa1
     1835  write(lunout,*)' qqa2 = ', qqa2
     1836  write(lunout,*)' gammas = ', gammas
     1837  write(lunout,*)' Fmax = ', Fmax
     1838  write(lunout,*)' alphas = ', alphas
     1839  write(lunout,*)' iflag_wake = ', iflag_wake
     1840  write(lunout,*)' alp_offset = ', alp_offset
     1841
     1842  write(lunout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
    18431843   lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
    1844   write(numout,*)' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
     1844  write(lunout,*)' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
    18451845   ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
    18461846
    1847   write(numout,*) 'ok_strato = ', ok_strato
    1848   write(numout,*) 'ok_hines = ',  ok_hines
    1849   write(numout,*) 'read_climoz = ', read_climoz
    1850   write(numout,*) 'carbon_cycle_tr = ', carbon_cycle_tr
    1851   write(numout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl
     1847  write(lunout,*) 'ok_strato = ', ok_strato
     1848  write(lunout,*) 'ok_hines = ',  ok_hines
     1849  write(lunout,*) 'read_climoz = ', read_climoz
     1850  write(lunout,*) 'carbon_cycle_tr = ', carbon_cycle_tr
     1851  write(lunout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl
    18521852 
    18531853!$OMP END MASTER
    … …  
    18661866   use IOIPSL
    18671867   implicit none
    1868 
     1868   include "iniprint.h"
    18691869! Configuration de l'interace atm/surf
    18701870!
    … …  
    18741874  REAL,SAVE     :: tau_calv_omp
    18751875
    1876 ! Local
    1877   integer              :: numout = 6
    18781876!
    18791877!Config Key  = tau_calv
    … …  
    18911889 
    18921890!$OMP MASTER
    1893   write(numout,*)' ##############################################'
    1894   WRITE(numout,*)' Configuration de l''interface atm/surfaces  : '
    1895   WRITE(numout,*)' tau_calv = ',tau_calv
     1891  write(lunout,*)' ##############################################'
     1892  WRITE(lunout,*)' Configuration de l''interface atm/surfaces  : '
     1893  WRITE(lunout,*)' tau_calv = ',tau_calv
    18961894!$OMP END MASTER
    18971895

  • LMDZ5/trunk/libf/phylmd/fisrtilp.F90

    r1507 r1575  
    147147  IF (appel1er) THEN
    148148     !
    149      PRINT*, 'fisrtilp, ninter:', ninter
    150      PRINT*, 'fisrtilp, evap_prec:', evap_prec
    151      PRINT*, 'fisrtilp, cpartiel:', cpartiel
     149     WRITE(lunout,*) 'fisrtilp, ninter:', ninter
     150     WRITE(lunout,*) 'fisrtilp, evap_prec:', evap_prec
     151     WRITE(lunout,*) 'fisrtilp, cpartiel:', cpartiel
    152152     IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN
    153         PRINT*, 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
    154         PRINT*, 'Je prefere un sous-intervalle de 6 minutes'
     153        WRITE(lunout,*) 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
     154        WRITE(lunout,*) 'Je prefere un sous-intervalle de 6 minutes'
    155155        !         CALL abort
    156156     ENDIF
    … …  
    616616
    617617  if (ncoreczq>0) then
    618      print*,'WARNING : ZQ dans fisrtilp ',ncoreczq,' val < 1.e-15.'
     618     WRITE(lunout,*)'WARNING : ZQ dans fisrtilp ',ncoreczq,' val < 1.e-15.'
    619619  endif
    620620

  • LMDZ5/trunk/libf/phylmd/fisrtilp_tr.F

    r1403 r1575  
    2424#include "YOMCST.h"
    2525#include "tracstoke.h"
     26#include "iniprint.h"
    2627c
    2728c Arguments:
    … …  
    136137      IF (appel1er) THEN
    137138c
    138          PRINT*, 'fisrtilp, calcrat:', calcrat
    139          PRINT*, 'fisrtilp, ninter:', ninter
    140          PRINT*, 'fisrtilp, evap_prec:', evap_prec
    141          PRINT*, 'fisrtilp, cpartiel:', cpartiel
     139         WRITE(lunout,*) 'fisrtilp, calcrat:', calcrat
     140         WRITE(lunout,*) 'fisrtilp, ninter:', ninter
     141         WRITE(lunout,*) 'fisrtilp, evap_prec:', evap_prec
     142         WRITE(lunout,*) 'fisrtilp, cpartiel:', cpartiel
    142143         IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN
    143           PRINT*, 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
    144           PRINT*, 'Je prefere un sous-intervalle de 6 minutes'
     144          WRITE(lunout,*)
     145     $    'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
     146          WRITE(lunout,*) 'Je prefere un sous-intervalle de 6 minutes'
    145147          CALL abort
    146148         ENDIF

  • LMDZ5/trunk/libf/phylmd/hgardfou.F

    r1550 r1575  
    1111#include "YOMCST.h"
    1212#include "indicesol.h"
     13#include "iniprint.h"
    1314      REAL t(klon,klev), tsol(klon,nbsrf)
    1415      CHARACTER*(*) text
    … …  
    2728
    2829      IF (firstcall) THEN
    29          PRINT*, 'hgardfou garantit la temperature dans [100,370] K'
     30         WRITE(lunout,*)
     31     $  'hgardfou garantit la temperature dans [100,370] K'
    3032         firstcall = .FALSE.
    3133c        DO i = 1, klon
    32 c         print*,'i=',i,'rlon=',rlon(i),'rlat=',rlat(i)
     34c         WRITE(lunout,*)'i=',i,'rlon=',rlon(i),'rlat=',rlat(i)
    3335c        ENDDO
    3436c
    … …  
    5456           ok = .FALSE.
    5557           DO i = 1, jbad
    56              PRINT *,'i,k,temperature,lon,lat,pourc ter,oce,lic,sic =',
     58             WRITE(lunout,*)
     59     $       'i,k,temperature,lon,lat,pourc ter,oce,lic,sic =',
    5760     $       jadrs(i),k,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i)),
    5861     $       (pctsrf(jadrs(i),nsrf),nsrf=1,nbsrf)
    … …  
    7477           ok = .FALSE.
    7578           DO i = 1, jbad
    76              PRINT *,'i,k,temperature,lon,lat,pourc ter,oce,lic,sic =',
     79             WRITE(lunout,*)
     80     $       'i,k,temperature,lon,lat,pourc ter,oce,lic,sic =',
    7781     $       jadrs(i),k,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i)),
    7882     $       (pctsrf(jadrs(i),nsrf),nsrf=1,nbsrf)
    … …  
    99103           ok = .FALSE.
    100104           DO i = 1, jbad
    101             PRINT *,'i,nsrf,temperature,lon,lat,pourc ter,oce,lic,sic ='
     105            WRITE(lunout,*)
     106     $      'i,nsrf,temperature,lon,lat,pourc ter,oce,lic,sic ='
    102107     $      ,jadrs(i),nsrf,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i))
    103108     $      ,pctsrf(jadrs(i),nsrf)
    … …  
    119124           ok = .FALSE.
    120125           DO i = 1, jbad
    121             PRINT *,'i,nsrf,temperature,lon,lat,pourc ter,oce,lic,sic ='
     126            WRITE(lunout,*)
     127     $      'i,nsrf,temperature,lon,lat,pourc ter,oce,lic,sic ='
    122128     $      ,jadrs(i),nsrf,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i))
    123129     $      ,pctsrf(jadrs(i),nsrf)

  • LMDZ5/trunk/libf/phylmd/mod_phys_lmdz_mpi_data.F90

    r1001 r1575  
    167167  SUBROUTINE print_module_data
    168168  IMPLICIT NONE
    169  
    170  
    171     PRINT *, 'ii_begin =', ii_begin
    172     PRINT *, 'ii_end =', ii_end
    173     PRINT *, 'jj_begin =',jj_begin
    174     PRINT *, 'jj_end =', jj_end
    175     PRINT *, 'jj_nb =', jj_nb
    176     PRINT *, 'ij_begin =', ij_begin
    177     PRINT *, 'ij_end =', ij_end
    178     PRINT *, 'ij_nb =', ij_nb
    179     PRINT *, 'klon_mpi_begin =', klon_mpi_begin
    180     PRINT *, 'klon_mpi_end =', klon_mpi_end
    181     PRINT *, 'klon_mpi =', klon_mpi
    182     PRINT *, 'jj_para_nb =', jj_para_nb
    183     PRINT *, 'jj_para_begin =', jj_para_begin
    184     PRINT *, 'jj_para_end =', jj_para_end
    185     PRINT *, 'ii_para_begin =', ii_para_begin
    186     PRINT *, 'ii_para_end =', ii_para_end
    187     PRINT *, 'ij_para_nb =', ij_para_nb
    188     PRINT *, 'ij_para_begin =', ij_para_begin
    189     PRINT *, 'ij_para_end =', ij_para_end
    190     PRINT *, 'klon_mpi_para_nb =', klon_mpi_para_nb
    191     PRINT *, 'klon_mpi_para_begin =', klon_mpi_para_begin
    192     PRINT *, 'klon_mpi_para_end  =', klon_mpi_para_end
    193     PRINT *, 'mpi_rank =', mpi_rank
    194     PRINT *, 'mpi_size =', mpi_size
    195     PRINT *, 'mpi_root =', mpi_root
    196     PRINT *, 'is_mpi_root =', is_mpi_root
    197     PRINT *, 'is_north_pole =', is_north_pole
    198     PRINT *, 'is_south_pole =', is_south_pole
    199     PRINT *, 'COMM_LMDZ_PHY =', COMM_LMDZ_PHY
     169  INCLUDE "iniprint.h"
     170 
     171    WRITE(lunout,*) 'ii_begin =', ii_begin
     172    WRITE(lunout,*) 'ii_end =', ii_end
     173    WRITE(lunout,*) 'jj_begin =',jj_begin
     174    WRITE(lunout,*) 'jj_end =', jj_end
     175    WRITE(lunout,*) 'jj_nb =', jj_nb
     176    WRITE(lunout,*) 'ij_begin =', ij_begin
     177    WRITE(lunout,*) 'ij_end =', ij_end
     178    WRITE(lunout,*) 'ij_nb =', ij_nb
     179    WRITE(lunout,*) 'klon_mpi_begin =', klon_mpi_begin
     180    WRITE(lunout,*) 'klon_mpi_end =', klon_mpi_end
     181    WRITE(lunout,*) 'klon_mpi =', klon_mpi
     182    WRITE(lunout,*) 'jj_para_nb =', jj_para_nb
     183    WRITE(lunout,*) 'jj_para_begin =', jj_para_begin
     184    WRITE(lunout,*) 'jj_para_end =', jj_para_end
     185    WRITE(lunout,*) 'ii_para_begin =', ii_para_begin
     186    WRITE(lunout,*) 'ii_para_end =', ii_para_end
     187    WRITE(lunout,*) 'ij_para_nb =', ij_para_nb
     188    WRITE(lunout,*) 'ij_para_begin =', ij_para_begin
     189    WRITE(lunout,*) 'ij_para_end =', ij_para_end
     190    WRITE(lunout,*) 'klon_mpi_para_nb =', klon_mpi_para_nb
     191    WRITE(lunout,*) 'klon_mpi_para_begin =', klon_mpi_para_begin
     192    WRITE(lunout,*) 'klon_mpi_para_end  =', klon_mpi_para_end
     193    WRITE(lunout,*) 'mpi_rank =', mpi_rank
     194    WRITE(lunout,*) 'mpi_size =', mpi_size
     195    WRITE(lunout,*) 'mpi_root =', mpi_root
     196    WRITE(lunout,*) 'is_mpi_root =', is_mpi_root
     197    WRITE(lunout,*) 'is_north_pole =', is_north_pole
     198    WRITE(lunout,*) 'is_south_pole =', is_south_pole
     199    WRITE(lunout,*) 'COMM_LMDZ_PHY =', COMM_LMDZ_PHY
    200200 
    201201  END SUBROUTINE print_module_data

  • LMDZ5/trunk/libf/phylmd/mod_phys_lmdz_omp_data.F90

    r1403 r1575  
    9191  SUBROUTINE Print_module_data
    9292  IMPLICIT NONE
     93  INCLUDE "iniprint.h"
    9394
    9495!$OMP CRITICAL 
    95   PRINT *,'--------> TASK ',omp_rank
    96   PRINT *,'omp_size =',omp_size
    97   PRINT *,'omp_rank =',omp_rank
    98   PRINT *,'is_omp_root =',is_omp_root
    99   PRINT *,'klon_omp_para_nb =',klon_omp_para_nb
    100   PRINT *,'klon_omp_para_begin =',klon_omp_para_begin
    101   PRINT *,'klon_omp_para_end =',klon_omp_para_end   
    102   PRINT *,'klon_omp =',klon_omp
    103   PRINT *,'klon_omp_begin =',klon_omp_begin
    104   PRINT *,'klon_omp_end =',klon_omp_end   
     96  WRITE(lunout,*)'--------> TASK ',omp_rank
     97  WRITE(lunout,*)'omp_size =',omp_size
     98  WRITE(lunout,*)'omp_rank =',omp_rank
     99  WRITE(lunout,*)'is_omp_root =',is_omp_root
     100  WRITE(lunout,*)'klon_omp_para_nb =',klon_omp_para_nb
     101  WRITE(lunout,*)'klon_omp_para_begin =',klon_omp_para_begin
     102  WRITE(lunout,*)'klon_omp_para_end =',klon_omp_para_end   
     103  WRITE(lunout,*)'klon_omp =',klon_omp
     104  WRITE(lunout,*)'klon_omp_begin =',klon_omp_begin
     105  WRITE(lunout,*)'klon_omp_end =',klon_omp_end   
    105106!$OMP END CRITICAL
    106107

  • LMDZ5/trunk/libf/phylmd/mod_phys_lmdz_para.F90

    r1327 r1575  
    4747  USE mod_grid_phy_lmdz
    4848  IMPLICIT NONE
    49  
     49    INCLUDE "iniprint.h"
     50 
    5051    REAL :: Test_Field1d_glo(klon_glo,nbp_lev)
    5152    REAL :: tmp1d_glo(klon_glo,nbp_lev)
    … …  
    8081!$OMP MASTER 
    8182      Checksum=sum(Test_Field1d_glo-tmp1d_glo)
    82       PRINT *, "------> Checksum =",Checksum," MUST BE 0"
     83      WRITE(lunout,*) "------> Checksum =",Checksum," MUST BE 0"
    8384!$OMP END MASTER
    8485    ENDIF
    … …  
    9293!$OMP MASTER 
    9394      Checksum=sum(Test_Field1d_glo-tmp1d_glo)
    94       PRINT *, "------> Checksum =",Checksum," MUST BE 0"
     95      WRITE(lunout,*) "------> Checksum =",Checksum," MUST BE 0"
    9596!$OMP END MASTER
    9697    ENDIF
    … …  
    102103!$OMP MASTER 
    103104      Checksum=sum(Test_Field1d_glo*omp_size*mpi_size-tmp1d_glo)
    104       PRINT *, "------> Checksum =",Checksum," MUST BE 0"
     105      WRITE(lunout,*) "------> Checksum =",Checksum," MUST BE 0"
    105106!$OMP END MASTER
    106107    ENDIF

  • LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

    r1574 r1575  
    631631    include "thermcell.h"
    632632    include "comvert.h"
     633    include "iniprint.h"
    633634
    634635    real,dimension(klon),intent(in) :: rlon
    … …  
    687688    real, dimension(nfiles), save     :: phys_out_latmax        = (/ 90., 90., 90., 90., 90., 90. /)
    688689
    689     print*,'Debut phys_output_mod.F90'
     690    write(lunout,*) 'Debut phys_output_mod.F90'
    690691    ! Initialisations (Valeurs par defaut
    691692
    … …  
    761762    type_ecri_files(:)=type_ecri(:)
    762763
    763     print*,'phys_out_lonmin=',phys_out_lonmin
    764     print*,'phys_out_lonmax=',phys_out_lonmax
    765     print*,'phys_out_latmin=',phys_out_latmin
    766     print*,'phys_out_latmax=',phys_out_latmax
    767     print*,'phys_out_filenames=',phys_out_filenames
    768     print*,'phys_out_filetypes=',type_ecri
    769     print*,'phys_out_filekeys=',clef_files
    770     print*,'phys_out_filestations=',clef_stations
    771     print*,'phys_out_filelevels=',lev_files
     764    write(lunout,*)'phys_out_lonmin=',phys_out_lonmin
     765    write(lunout,*)'phys_out_lonmax=',phys_out_lonmax
     766    write(lunout,*)'phys_out_latmin=',phys_out_latmin
     767    write(lunout,*)'phys_out_latmax=',phys_out_latmax
     768    write(lunout,*)'phys_out_filenames=',phys_out_filenames
     769    write(lunout,*)'phys_out_filetypes=',type_ecri
     770    write(lunout,*)'phys_out_filekeys=',clef_files
     771    write(lunout,*)'phys_out_filestations=',clef_stations
     772    write(lunout,*)'phys_out_filelevels=',lev_files
    772773
    773774!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    785786    enddo
    786787    !          if(prt_level.ge.1) then
    787     print*,'Ap Hybrid = ',Ahyb(1:klev)
    788     print*,'Bp Hybrid = ',Bhyb(1:klev)
    789     print*,'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
     788    write(lunout,*)'Ap Hybrid = ',Ahyb(1:klev)
     789    write(lunout,*)'Bp Hybrid = ',Bhyb(1:klev)
     790    write(lunout,*)'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
    790791    !          endif
    791792    DO iff=1,nfiles
    … …  
    799800             call convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))
    800801          endif
    801           print*,'ecrit_files(',iff,')= ',ecrit_files(iff)
     802          write(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
    802803
    803804          zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
    … …  
    817818             ! correction abderr       
    818819             do i=1,iim
    819                 print*,'io_lon(i)=',io_lon(i)
     820                write(lunout,*)'io_lon(i)=',io_lon(i)
    820821                if (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
    821822                if (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1
    … …  
    823824
    824825             do j=1,jjmp1
    825                 print*,'io_lat(j)=',io_lat(j)
     826                write(lunout,*)'io_lat(j)=',io_lat(j)
    826827                if (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1
    827828                if (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j
    828829             enddo
    829830
    830              print*,'On stoke le fichier histoire numero ',iff,' sur ', &
     831             write(lunout,*)'On stoke le fichier histoire numero ',iff,' sur ', &
    831832                  imin_ins,imax_ins,jmin_ins,jmax_ins
    832              print*,'longitudes : ', &
     833             write(lunout,*)'longitudes : ', &
    833834                  io_lon(imin_ins),io_lon(imax_ins), &
    834835                  'latitudes : ', &
    … …  
    843844          else if (clef_stations(iff)) THEN
    844845
    845              print*,'phys_output_mod phys_out_filenames=',phys_out_filenames(iff)
     846             write(lunout,*)'phys_output_mod phys_out_filenames=',phys_out_filenames(iff)
    846847
    847848             call histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &
    … …  
    17091710
    17101711    ENDDO !  iff
    1711     print*,'Fin phys_output_mod.F90'
     1712    write(lunout,*)'Fin phys_output_mod.F90'
    17121713  end subroutine phys_output_open
    17131714
    … …  
    18411842    include "temps.h"
    18421843    include "comconst.h"
     1844    include "iniprint.h"
    18431845
    18441846    ipos=scan(str,'0123456789.',.true.)
    18451847    ! 
    18461848    il=len_trim(str)
    1847     print*,ipos,il
     1849    write(lunout,*)ipos,il
    18481850    read(str(1:ipos),*) ttt
    1849     print*,ttt
     1851    write(lunout,*)ttt
    18501852    type=str(ipos+1:il)
    18511853
    … …  
    18571859    if ( type == 'day'.or.type == 'days'.or.type == 'jours'.or.type == 'jour' ) timestep = ttt * dayseconde
    18581860    if ( type == 'mounths'.or.type == 'mth'.or.type == 'mois' ) then
    1859        print*,'annee_ref,day_ref mon_len',annee_ref,day_ref,ioget_mon_len(annee_ref,day_ref)
     1861       write(lunout,*)'annee_ref,day_ref mon_len',annee_ref,day_ref,ioget_mon_len(annee_ref,day_ref)
    18601862       timestep = ttt * dayseconde * mth_len
    18611863    endif
    … …  
    18651867    if ( type == 'TS' ) timestep = ttt * dtime
    18661868
    1867     print*,'type =      ',type
    1868     print*,'nb j/h/m =  ',ttt
    1869     print*,'timestep(s)=',timestep
     1869    write(lunout,*)'type =      ',type
     1870    write(lunout,*)'nb j/h/m =  ',ttt
     1871    write(lunout,*)'timestep(s)=',timestep
    18701872
    18711873  END SUBROUTINE convers_timesteps

  • LMDZ5/trunk/libf/phylmd/phytrac.F90

    r1565 r1575  
    4646  INCLUDE "paramet.h"
    4747  INCLUDE "thermcell.h"
     48  INCLUDE "iniprint.h"
    4849!==========================================================================
    4950!                   -- ARGUMENT DESCRIPTION --
    … …  
    202203!######################################################################
    203204  IF (debutphy) THEN
    204      WRITE(*,*) 'FIRST TIME IN PHYTRAC : pdtphys(sec) = ',pdtphys,'ecrit_tra (sec) = ',ecrit_tra
     205     IF (prt_level >9) WRITE(lunout,*) 'FIRST TIME IN PHYTRAC : pdtphys(sec) = ',pdtphys,'ecrit_tra (sec) = ',ecrit_tra
    205206     ALLOCATE( source(klon,nbtr), stat=ierr)
    206207     IF (ierr /= 0) CALL abort_gcm('phytrac', 'pb in allocation 1',1)
    … …  
    224225! ----------------------------
    225226#ifdef CPP_IOIPSL
    226 !     INCLUDE "ini_histrac.h"
     227     INCLUDE "ini_histrac.h"
    227228#endif
    228229  END IF
    … …  
    414415!=============================================================
    415416#ifdef CPP_IOIPSL
    416 !  INCLUDE "write_histrac.h"
     417  INCLUDE "write_histrac.h"
    417418#endif
    418419

  • LMDZ5/trunk/libf/phylmd/radiation_AR4.F

    r1568 r1575  
    1919cym#include "raddim.h"
    2020#include "YOMCST.h"
     21#include "iniprint.h"
    2122C
    2223C     ------------------------------------------------------------------
    … …  
    171172
    172173      IF (appel1er) THEN
    173          PRINT*, 'SW calling frequency : ', swpas
    174          PRINT*, "   In general, it should be 1"
     174         WRITE(lunout,*) 'SW calling frequency : ', swpas
     175         WRITE(lunout,*) "   In general, it should be 1"
    175176         appel1er = .FALSE.
    176177      ENDIF
    … …  
    531532cym#include "dimensions.h"
    532533cym#include "dimphy.h"
    533 cym#include "raddim.h"
     534cym#include "raddim.h"i
     535#include "iniprint.h"
    534536C
    535537C     ------------------------------------------------------------------
    … …  
    623625            RSUN(2) = RSUNTIME(2)
    624626         ENDIF
    625          PRINT*,'RSUN(1): ',RSUN(1)
     627         WRITE(lunout,*)'RSUN(1): ',RSUN(1)
    626628#endif
    627629      END IF
    … …  
    23382340#include "raddimlw.h"
    23392341#include "YOMCST.h"
     2342#include "iniprint.h"
    23402343C
    23412344C-----------------------------------------------------------------------
    … …  
    24432446C     ------------------------------------------------------------------
    24442447      IF (appel1er) THEN
    2445          PRINT*, "LW clear-sky calling frequency: ", lw0pas
    2446          PRINT*, "LW cloudy-sky calling frequency: ", lwpas
    2447          PRINT*, "   In general, they should be 1"
     2448         WRITE(lunout,*) "LW clear-sky calling frequency: ", lw0pas
     2449         WRITE(lunout,*) "LW cloudy-sky calling frequency: ", lwpas
     2450         WRITE(lunout,*) "   In general, they should be 1"
    24482451cym
    24492452         allocate(ZFLUX(KDLON,2,KFLEV+1) )

  • LMDZ5/trunk/libf/phylmd/soil.F90

    r996 r1575  
    5454  INCLUDE "indicesol.h"
    5555  INCLUDE "comsoil.h"
     56  INCLUDE "iniprint.h"
    5657!-----------------------------------------------------------------------
    5758! Arguments
    … …  
    111112           READ(99,*) min_period
    112113           READ(99,*) dalph_soil
    113            PRINT*,'Discretization for the soil model'
    114            PRINT*,'First level e-folding depth',min_period, &
     114           WRITE(lunout,*)'Discretization for the soil model'
     115           WRITE(lunout,*)'First level e-folding depth',min_period, &
    115116                '   dalph',dalph_soil
    116117           CLOSE(99)
    … …  
    135136     ENDDO
    136137     lambda=fz(.5)*dz1(1)
    137      PRINT*,'full layers, intermediate layers (seconds)'
     138     WRITE(lunout,*)'full layers, intermediate layers (seconds)'
    138139     DO jk=1,nsoilmx
    139140        rk=jk
    140141        rk1=jk+.5
    141142        rk2=jk-.5
    142         PRINT *,'fz=', &
     143        WRITE(lunout,*)'fz=', &
    143144             fz(rk1)*fz(rk2)*3.14,fz(rk)*fz(rk)*3.14
    144145     ENDDO
    … …  
    175176     ENDDO
    176177  ELSE
    177      PRINT*, "valeur d indice non prevue", indice
     178     WRITE(lunout,*) "valeur d indice non prevue", indice
    178179     CALL abort
    179180  ENDIF

  • LMDZ5/trunk/libf/phylmd/sw_aeroAR4.F90

    r1307 r1575  
    2626#include "YOMCST.h"
    2727#include "clesphys.h"
     28#include "iniprint.h"
    2829  !
    2930  !     ------------------------------------------------------------------
    … …  
    225226
    226227  IF (appel1er) THEN
    227      PRINT*, 'SW calling frequency : ', swpas
    228      PRINT*, "   In general, it should be 1"
     228     WRITE(lunout,*) 'SW calling frequency : ', swpas
     229     WRITE(lunout,*) "   In general, it should be 1"
    229230     appel1er = .FALSE.
    230231  ENDIF
Note: See TracChangeset for help on using the changeset viewer.