Changeset 1664 for LMDZ5/branches/testing/libf/phylmd/conf_phys.F90
- Timestamp:
- Oct 9, 2012, 3:29:15 PM (12 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk (added) merged: 1555-1557,1559-1566,1568-1575
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/conf_phys.F90
r1539 r1664 37 37 38 38 include "thermcell.h" 39 include "iniprint.h" 39 40 40 41 !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12 … … 101 102 102 103 ! Local 103 integer :: numout = 6104 104 real :: zzz 105 105 … … 199 199 call getin('OCEAN', ocean_omp) 200 200 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' 204 204 CALL abort_gcm('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1) 205 205 END IF … … 1351 1351 !Config Help = 1352 1352 ! 1353 ecrit_tra_omp = 30.1353 ecrit_tra_omp = 0. 1354 1354 call getin('ecrit_tra',ecrit_tra_omp) 1355 1355 ! … … 1696 1696 ! Test of coherence between type_ocean and version_ocean 1697 1697 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' 1699 1699 CALL abort_gcm('conf_phys','version_ocean not valid',1) 1700 1700 END IF … … 1703 1703 version_ocean='sicOBS' 1704 1704 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' 1706 1706 CALL abort_gcm('conf_phys','version_ocean not valid',1) 1707 1707 END IF … … 1717 1717 !$OMP MASTER 1718 1718 1719 write( numout,*)' ##############################################'1720 write( numout,*)' Configuration des parametres de la physique: '1721 write( numout,*)' Type ocean = ', type_ocean1722 write( numout,*)' Version ocean = ', version_ocean1723 write( numout,*)' Config veget = ', ok_veget1724 write( numout,*)' Sortie journaliere = ', ok_journe1725 write( numout,*)' Sortie haute frequence = ', ok_hf1726 write( numout,*)' Sortie mensuelle = ', ok_mensuel1727 write( numout,*)' Sortie instantanee = ', ok_instan1728 write( numout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP1729 write( numout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP1730 write( numout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP1731 write( numout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy1732 write( numout,*)' Excentricite = ',R_ecc1733 write( numout,*)' Equinoxe = ',R_peri1734 write( numout,*)' Inclinaison =',R_incl1735 write( numout,*)' Constante solaire =',solaire1736 write( numout,*)' co2_ppm =',co2_ppm1737 write( numout,*)' RCO2_act = ',RCO2_act1738 write( numout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act1739 write( numout,*)' N2O_ppb =',N2O_ppb,' RN2O_act= ',RN2O_act1740 write( numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act= ',RCFC11_act1741 write( numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act= ',RCFC12_act1742 write( numout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per1743 write( numout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per1744 write( numout,*)' RCFC12_per = ',RCFC12_per1745 write( numout,*)' cvl_corr=', cvl_corr1746 write( numout,*)'ok_lic_melt=', ok_lic_melt1747 write( numout,*)'cycle_diurne=',cycle_diurne1748 write( numout,*)'soil_model=',soil_model1749 write( numout,*)'new_oliq=',new_oliq1750 write( numout,*)'ok_orodr=',ok_orodr1751 write( numout,*)'ok_orolf=',ok_orolf1752 write( numout,*)'ok_limitvrai=',ok_limitvrai1753 write( numout,*)'nbapp_rad=',nbapp_rad1754 write( numout,*)'iflag_con=',iflag_con1755 write( numout,*)' epmax = ', epmax1756 write( numout,*)' ok_adj_ema = ', ok_adj_ema1757 write( numout,*)' iflag_clw = ', iflag_clw1758 write( numout,*)' cld_lc_lsc = ', cld_lc_lsc1759 write( numout,*)' cld_lc_con = ', cld_lc_con1760 write( numout,*)' cld_tau_lsc = ', cld_tau_lsc1761 write( numout,*)' cld_tau_con = ', cld_tau_con1762 write( numout,*)' ffallv_lsc = ', ffallv_lsc1763 write( numout,*)' ffallv_con = ', ffallv_con1764 write( numout,*)' coef_eva = ', coef_eva1765 write( numout,*)' reevap_ice = ', reevap_ice1766 write( numout,*)' iflag_pdf = ', iflag_pdf1767 write( numout,*)' iflag_cldcon = ', iflag_cldcon1768 write( numout,*)' iflag_radia = ', iflag_radia1769 write( numout,*)' iflag_rrtm = ', iflag_rrtm1770 write( numout,*)' iflag_ratqs = ', iflag_ratqs1771 write( numout,*)' seuil_inversion = ', seuil_inversion1772 write( numout,*)' fact_cldcon = ', fact_cldcon1773 write( numout,*)' facttemps = ', facttemps1774 write( numout,*)' ok_newmicro = ',ok_newmicro1775 write( numout,*)' ratqsbas = ',ratqsbas1776 write( numout,*)' ratqshaut = ',ratqshaut1777 write( numout,*)' tau_ratqs = ',tau_ratqs1778 write( numout,*)' top_height = ',top_height1779 write( numout,*)' rad_froid = ',rad_froid1780 write( numout,*)' rad_chau1 = ',rad_chau11781 write( numout,*)' rad_chau2 = ',rad_chau21782 write( numout,*)' t_glace_min = ',t_glace_min1783 write( numout,*)' t_glace_max = ',t_glace_max1784 write( numout,*)' rei_min = ',rei_min1785 write( numout,*)' rei_max = ',rei_max1786 write( numout,*)' overlap = ',overlap1787 write( numout,*)' cdmmax = ',cdmmax1788 write( numout,*)' cdhmax = ',cdhmax1789 write( numout,*)' ksta = ',ksta1790 write( numout,*)' ksta_ter = ',ksta_ter1791 write( numout,*)' ok_kzmin = ',ok_kzmin1792 write( numout,*)' fmagic = ',fmagic1793 write( numout,*)' pmagic = ',pmagic1794 write( numout,*)' ok_ade = ',ok_ade1795 write( numout,*)' ok_aie = ',ok_aie1796 write( numout,*)' aerosol_couple = ', aerosol_couple1797 write( numout,*)' flag_aerosol = ', flag_aerosol1798 write( numout,*)' new_aod = ', new_aod1799 write( numout,*)' aer_type = ',aer_type1800 write( numout,*)' bl95_b0 = ',bl95_b01801 write( numout,*)' bl95_b1 = ',bl95_b11802 write( numout,*)' lev_histhf = ',lev_histhf1803 write( numout,*)' lev_histday = ',lev_histday1804 write( numout,*)' lev_histmth = ',lev_histmth1805 write( numout,*)' lev_histins = ',lev_histins1806 write( numout,*)' lev_histLES = ',lev_histLES1807 write( numout,*)' lev_histdayNMC = ',lev_histdayNMC1808 write( numout,*)' ok_histNMC = ',ok_histNMC1809 write( numout,*)' freq_outNMC = ',freq_outNMC1810 write( numout,*)' freq_calNMC = ',freq_calNMC1811 write( numout,*)' iflag_pbl = ', iflag_pbl1812 write( numout,*)' iflag_thermals = ', iflag_thermals1813 write( numout,*)' iflag_thermals_ed = ', iflag_thermals_ed1814 write( numout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux1815 write( numout,*)' iflag_clos = ', iflag_clos1816 write( numout,*)' type_run = ',type_run1817 write( numout,*)' ok_isccp = ',ok_isccp1818 write( numout,*)' ok_cosp = ',ok_cosp1819 write( numout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP1820 write( numout,*)' ok_journeCOSP = ',ok_journeCOSP1821 write( numout,*)' ok_hfCOSP =',ok_hfCOSP1822 write( numout,*)' solarlong0 = ', solarlong01823 write( numout,*)' qsol0 = ', qsol01824 write( numout,*)' inertie_sol = ', inertie_sol1825 write( numout,*)' inertie_ice = ', inertie_ice1826 write( numout,*)' inertie_sno = ', inertie_sno1827 write( numout,*)' f_cdrag_ter = ',f_cdrag_ter1828 write( numout,*)' f_cdrag_oce = ',f_cdrag_oce1829 write( numout,*)' f_rugoro = ',f_rugoro1830 write( numout,*)' supcrit1 = ', supcrit11831 write( numout,*)' supcrit2 = ', supcrit21832 write( numout,*)' iflag_mix = ', iflag_mix1833 write( numout,*)' scut = ', scut1834 write( numout,*)' qqa1 = ', qqa11835 write( numout,*)' qqa2 = ', qqa21836 write( numout,*)' gammas = ', gammas1837 write( numout,*)' Fmax = ', Fmax1838 write( numout,*)' alphas = ', alphas1839 write( numout,*)' iflag_wake = ', iflag_wake1840 write( numout,*)' alp_offset = ', alp_offset1841 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 =',& 1843 1843 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',& 1845 1845 ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES 1846 1846 1847 write( numout,*) 'ok_strato = ', ok_strato1848 write( numout,*) 'ok_hines = ', ok_hines1849 write( numout,*) 'read_climoz = ', read_climoz1850 write( numout,*) 'carbon_cycle_tr = ', carbon_cycle_tr1851 write( numout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl1847 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 1852 1852 1853 1853 !$OMP END MASTER … … 1866 1866 use IOIPSL 1867 1867 implicit none 1868 1868 include "iniprint.h" 1869 1869 ! Configuration de l'interace atm/surf 1870 1870 ! … … 1874 1874 REAL,SAVE :: tau_calv_omp 1875 1875 1876 ! Local1877 integer :: numout = 61878 1876 ! 1879 1877 !Config Key = tau_calv … … 1891 1889 1892 1890 !$OMP MASTER 1893 write( numout,*)' ##############################################'1894 WRITE( numout,*)' Configuration de l''interface atm/surfaces : '1895 WRITE( numout,*)' tau_calv = ',tau_calv1891 write(lunout,*)' ##############################################' 1892 WRITE(lunout,*)' Configuration de l''interface atm/surfaces : ' 1893 WRITE(lunout,*)' tau_calv = ',tau_calv 1896 1894 !$OMP END MASTER 1897 1895
Note: See TracChangeset
for help on using the changeset viewer.