Ignore:
Timestamp:
May 27, 2020, 2:59:10 PM (4 years ago)
Author:
fhourdin
Message:

Retropedalage sur le 1D.
Annulation de toutes les modifs de la veille.
Frederic

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

    r3682 r3686  
    627627      RETURN
    628628      END
    629 !     
    630       SUBROUTINE conf_unicol_std
    631 !
    632 #ifdef CPP_IOIPSL
    633       use IOIPSL
    634 #else
    635 ! if not using IOIPSL, we still need to use (a local version of) getin
    636       use ioipsl_getincom
    637 #endif
    638       USE print_control_mod, ONLY: lunout
    639       IMPLICIT NONE
    640 !-----------------------------------------------------------------------
    641 !     Auteurs :   A. Lahellec  - adaptation au format standard.
    642 !
    643 !   Declarations :
    644 !   --------------
    645 
    646 #include "compar1d_std.h"
    647 #include "flux_arp.h"
    648 #include "tsoilnudge.h"
    649 #include "fcg_gcssold.h"
    650 #include "fcg_racmo.h"
    651 !
    652 !
    653 !   local:
    654 !   ------
    655 
    656 !      CHARACTER ch1*72,ch2*72,ch3*72,ch4*12
    657      
    658 !
    659 !  -------------------------------------------------------------------
    660 !
    661 !      .........    Initilisation parametres du lmdz1D      ..........
    662 !
    663 !---------------------------------------------------------------------
    664 !   initialisations:
    665 !   ----------------
    666 
    667 !Config  Key  = lunout
    668 !Config  Desc = unite de fichier pour les impressions
    669 !Config  Def  = 6
    670 !Config  Help = unite de fichier pour les impressions
    671 !Config         (defaut sortie standard = 6)
    672       lunout=6
    673 !      CALL getin('lunout', lunout)
    674       IF (lunout /= 5 .and. lunout /= 6) THEN
    675         OPEN(lunout,FILE='lmdz.out')
    676       ENDIF
    677 
    678 !Config  Key  = prt_level
    679 !Config  Desc = niveau d'impressions de debogage
    680 !Config  Def  = 0
    681 !Config  Help = Niveau d'impression pour le debogage
    682 !Config         (0 = minimum d'impression)
    683 !      prt_level = 0
    684 !      CALL getin('prt_level',prt_level)
    685 
    686 !-----------------------------------------------------------------------
    687 !  Parametres de controle du run:
    688 !-----------------------------------------------------------------------
    689 
    690 !Config  Key  = restart
    691 !Config  Desc = on repart des startphy et start1dyn
    692 !Config  Def  = false
    693 !Config  Help = les fichiers restart doivent etre renomme en start
    694        restart =.false.
    695        CALL getin('restart',restart)
    696 
    697 !----------------------------------------------------------
    698 ! Parametres de forcage pour les forcages communs:
    699 ! Voir ici: https://github.com/romainroehrig/DEPHY-SCM/blob/master/DEPHY_Format_Version_0.pdf
    700 ! Pour les forcages communs: ces entiers valent 0 ou 1
    701 ! adv_temp= advection tempe, theta ou thetal, qv,qt,rv ou rt
    702 ! rad_temp= 0 (rayonnement actif) ou 1 (prescrit par tend_rad) ou adv (prescir et contenu dans les adv_temp)
    703 ! idem rad_theta et rad_thetal 
    704 ! forcages en omega, w, vent geostrophique ou ustar
    705 ! Parametres de nudging en u,v,temp, theta, thetal,qv, qt, rv, rt valent 0 ou 1 ou le temps de nudging
    706 ! p_nudging_xxx pression (Pa) a partir de laquelle appliquer le nudging de xxx
    707 ! ou z_nudging_xxx hauteur(m) a partir de laquelle appliquer le nudging de xxx
    708 !----------------------------------------------------------
    709 !
    710 !Parametres de forcage
    711 !Config  Key  = adv_temp
    712 !Config  Desc = forcage ou non par advection de T
    713 !Config  Def  = false
    714 !Config  Help = forcage ou non par advection de T
    715        adv_temp =0
    716        CALL getin('adv_temp',adv_temp)
    717 
    718 !
    719 !Parametres de forcage
    720 !Config  Key  = adv_theta
    721 !Config  Desc = forcage ou non par advection de Theta
    722 !Config  Def  = false
    723 !Config  Help = forcage ou non par advection de Theta
    724        adv_theta =0
    725        CALL getin('adv_theta',adv_theta)
    726 
    727 !
    728 !Parametres de forcage
    729 !Config  Key  = adv_thetal
    730 !Config  Desc = forcage ou non par advection de Thetal
    731 !Config  Def  = false
    732 !Config  Help = forcage ou non par advection de Thetal
    733        adv_thetal =0
    734        CALL getin('adv_thetal',adv_thetal)
    735 
    736 !
    737 !Parametres de forcage
    738 !Config  Key  = rad_temp
    739 !Config  Desc = forcage par tendance radiative en tempe
    740 !Config  Def  = false
    741 !Config  Help = forcage par tendance radiative en tempe
    742        rad_temp ="0"
    743        CALL getin('rad_temp',rad_temp)
    744 
    745 !
    746 !Parametres de forcage
    747 !Config  Key  = rad_theta
    748 !Config  Desc = forcage par tendance radiative en theta
    749 !Config  Def  = false
    750 !Config  Help = forcage par tendance radiative en theta
    751        rad_theta ="0"
    752        CALL getin('rad_theta',rad_theta)
    753 
    754 !
    755 !Parametres de forcage
    756 !Config  Key  = rad_thetal
    757 !Config  Desc = forcage par tendance radiative en thetal
    758 !Config  Def  = false
    759 !Config  Help = forcage par tendance radiative en thetal
    760        rad_thetal ="0"
    761        CALL getin('rad_thetal',rad_thetal)
    762 
    763 !
    764 !Parametres de forcage
    765 !Config  Key  = adv_qv
    766 !Config  Desc = forcage ou non par advection de qv
    767 !Config  Def  = false
    768 !Config  Help = forcage ou non par advection de qv
    769        adv_qv =0
    770        CALL getin('adv_qv',adv_qv)
    771 
    772 !
    773 !Parametres de forcage
    774 !Config  Key  = adv_qt
    775 !Config  Desc = forcage ou non par advection de qt
    776 !Config  Def  = false
    777 !Config  Help = forcage ou non par advection de qt
    778        adv_qt =0
    779        CALL getin('adv_qt',adv_qt)
    780 
    781 !
    782 !Parametres de forcage
    783 !Config  Key  = adv_rv
    784 !Config  Desc = forcage ou non par advection de rv
    785 !Config  Def  = false
    786 !Config  Help = forcage ou non par advection de rv
    787        adv_rv =0
    788        CALL getin('adv_rv',adv_rv)
    789 
    790 !
    791 !Parametres de forcage
    792 !Config  Key  = adv_rt
    793 !Config  Desc = forcage ou non par advection de rt
    794 !Config  Def  = false
    795 !Config  Help = forcage ou non par advection de rt
    796        adv_rt =0
    797        CALL getin('adv_rt',adv_rt)
    798 
    799 !
    800 !Parametres de forcage
    801 !Config  Key  = nudging_temp
    802 !Config  Desc = forcage ou non par advection de tempe
    803 !Config  Def  = false
    804 !Config  Help = forcage ou non par advection de tempe
    805        nudging_temp =0
    806        CALL getin('nudging_temp',nudging_temp)
    807 
    808 !
    809 !Parametres de forcage
    810 !Config  Key  = nudging_theta
    811 !Config  Desc = forcage ou non par advection de theta
    812 !Config  Def  = false
    813 !Config  Help = forcage ou non par advection de theta
    814        nudging_theta =0
    815        CALL getin('nudging_theta',nudging_theta)
    816 
    817 !
    818 !Parametres de forcage
    819 !Config  Key  = nudging_thetal
    820 !Config  Desc = forcage ou non par advection de thetal
    821 !Config  Def  = false
    822 !Config  Help = forcage ou non par advection de thetal
    823        nudging_thetal =0
    824        CALL getin('nudging_thetal',nudging_thetal)
    825 
    826 !
    827 !Parametres de forcage
    828 !Config  Key  = nudging_qv
    829 !Config  Desc = forcage ou non par advection de qv
    830 !Config  Def  = false
    831 !Config  Help = forcage ou non par advection de qv
    832        nudging_qv =0
    833        CALL getin('nudging_qv',nudging_qv)
    834 
    835 !
    836 !Parametres de forcage
    837 !Config  Key  = nudging_qt
    838 !Config  Desc = forcage ou non par advection de qt
    839 !Config  Def  = false
    840 !Config  Help = forcage ou non par advection de qt
    841        nudging_qt =0
    842        CALL getin('nudging_qt',nudging_qt)
    843 
    844 !
    845 !Parametres de forcage
    846 !Config  Key  = nudging_rv
    847 !Config  Desc = forcage ou non par advection de rv
    848 !Config  Def  = false
    849 !Config  Help = forcage ou non par advection de rv
    850        nudging_rv =0
    851        CALL getin('nudging_rv',nudging_rv)
    852 
    853 !
    854 !Parametres de forcage
    855 !Config  Key  = nudging_rt
    856 !Config  Desc = forcage ou non par advection de rt
    857 !Config  Def  = false
    858 !Config  Help = forcage ou non par advection de rt
    859        nudging_rt =0
    860        CALL getin('nudging_rt',nudging_rt)
    861 
    862 !
    863 !Parametres de forcage
    864 !Config  Key  = nudging_u
    865 !Config  Desc = forcage ou non par advection de u wind
    866 !Config  Def  = false
    867 !Config  Help = forcage ou non par advection de u wind
    868        nudging_u =0
    869        CALL getin('nudging_u',nudging_u)
    870 
    871 !
    872 !Parametres de forcage
    873 !Config  Key  = nudging_v
    874 !Config  Desc = forcage ou non par advection de v wind
    875 !Config  Def  = false
    876 !Config  Help = forcage ou non par advection de v wind
    877        nudging_v =0
    878        CALL getin('nudging_v',nudging_v)
    879 
    880 !
    881 !Parametres de forcage
    882 !Config  Key  = p_nudging_temp
    883 !Config  Desc = Pressure (Pa) above which tempe should be nudged
    884 !Config  Def  = false
    885 !Config  Help = Pressure (Pa) above which tempe should be nudged
    886        p_nudging_temp =11000.
    887        CALL getin('p_nudging_tempe',p_nudging_temp)
    888 
    889 !
    890 !Parametres de forcage
    891 !Config  Key  = p_nudging_theta
    892 !Config  Desc = Pressure (Pa) above which theta should be nudged
    893 !Config  Def  = false
    894 !Config  Help = Pressure (Pa) above which theta should be nudged
    895        p_nudging_theta =11000.
    896        CALL getin('p_nudging_theta',p_nudging_theta)
    897 
    898 !
    899 !Parametres de forcage
    900 !Config  Key  = p_nudging_thetal
    901 !Config  Desc = Pressure (Pa) above which thetal should be nudged
    902 !Config  Def  = false
    903 !Config  Help = Pressure (Pa) above which thetal should be nudged
    904        p_nudging_thetal =11000.
    905        CALL getin('p_nudging_thetal',p_nudging_thetal)
    906 
    907 !
    908 !Parametres de forcage
    909 !Config  Key  = p_nudging_qv
    910 !Config  Desc = Pressure (Pa) above which qv should be nudged
    911 !Config  Def  = false
    912 !Config  Help = Pressure (Pa) above which qv should be nudged
    913        p_nudging_qv =11000.
    914        CALL getin('p_nudging_qv',p_nudging_qv)
    915 
    916 !
    917 !Parametres de forcage
    918 !Config  Key  = p_nudging_qt
    919 !Config  Desc = Pressure (Pa) above which qt should be nudged
    920 !Config  Def  = false
    921 !Config  Help = Pressure (Pa) above which qt should be nudged
    922        p_nudging_qt =11000.
    923        CALL getin('p_nudging_qt',p_nudging_qt)
    924 
    925 !
    926 !Parametres de forcage
    927 !Config  Key  = p_nudging_rv
    928 !Config  Desc = Pressure (Pa) above which rv should be nudged
    929 !Config  Def  = false
    930 !Config  Help = Pressure (Pa) above which rv should be nudged
    931        p_nudging_rv =11000.
    932        CALL getin('p_nudging_rv',p_nudging_rv)
    933 
    934 !
    935 !Parametres de forcage
    936 !Config  Key  = p_nudging_rt
    937 !Config  Desc = Pressure (Pa) above which rt should be nudged
    938 !Config  Def  = false
    939 !Config  Help = Pressure (Pa) above which rt should be nudged
    940        p_nudging_rt =11000.
    941        CALL getin('p_nudging_rt',p_nudging_rt)
    942 
    943 !
    944 !Parametres de forcage
    945 !Config  Key  = p_nudging_u
    946 !Config  Desc = Pressure (Pa) above which u should be nudged
    947 !Config  Def  = false
    948 !Config  Help = Pressure (Pa) above which u should be nudged
    949        p_nudging_u =11000.
    950        CALL getin('p_nudging_u',p_nudging_u)
    951 
    952 !
    953 !Parametres de forcage
    954 !Config  Key  = p_nudging_v
    955 !Config  Desc = Pressure (Pa) above which v should be nudged
    956 !Config  Def  = false
    957 !Config  Help = Pressure (Pa) above which v should be nudged
    958        p_nudging_v =11000.
    959        CALL getin('p_nudging_v',p_nudging_v)
    960 
    961 !
    962 !Parametres de forcage
    963 !Config  Key  = z_nudging_temp
    964 !Config  Desc = Height (m) above which tempe should be nudged
    965 !Config  Def  = false
    966 !Config  Help = Height (m) above which tempe should be nudged
    967        z_nudging_temp=0.
    968        CALL getin('z_nudging_tempe',z_nudging_temp)
    969 
    970 !
    971 !Parametres de forcage
    972 !Config  Key  = z_nudging_theta
    973 !Config  Desc = Height (m) above which theta should be nudged
    974 !Config  Def  = false
    975 !Config  Help = Height (m) above which theta should be nudged
    976        z_nudging_theta=0.
    977        CALL getin('z_nudging_theta',z_nudging_theta)
    978 
    979 !
    980 !Parametres de forcage
    981 !Config  Key  = z_nudging_thetal
    982 !Config  Desc = Height (m) above which thetal should be nudged
    983 !Config  Def  = false
    984 !Config  Help = Height (m) above which thetal should be nudged
    985        z_nudging_thetal=0.
    986        CALL getin('z_nudging_thetal',z_nudging_thetal)
    987 
    988 !
    989 !Parametres de forcage
    990 !Config  Key  = z_nudging_qv
    991 !Config  Desc = Height (m) above which qv should be nudged
    992 !Config  Def  = false
    993 !Config  Help = Height (m) above which qv should be nudged
    994        z_nudging_qv=0.
    995        CALL getin('z_nudging_qv',z_nudging_qv)
    996 
    997 !
    998 !Parametres de forcage
    999 !Config  Key  = z_nudging_qt
    1000 !Config  Desc = Height (m) above which qt should be nudged
    1001 !Config  Def  = false
    1002 !Config  Help = Height (m) above which qt should be nudged
    1003        z_nudging_qt=0.
    1004        CALL getin('z_nudging_qt',z_nudging_qt)
    1005 
    1006 !
    1007 !Parametres de forcage
    1008 !Config  Key  = z_nudging_rv
    1009 !Config  Desc = Height (m) above which rv should be nudged
    1010 !Config  Def  = false
    1011 !Config  Help = Height (m) above which rv should be nudged
    1012        z_nudging_rv=0.
    1013        CALL getin('z_nudging_rv',z_nudging_rv)
    1014 
    1015 !
    1016 !Parametres de forcage
    1017 !Config  Key  = z_nudging_rt
    1018 !Config  Desc = Height (m) above which rt should be nudged
    1019 !Config  Def  = false
    1020 !Config  Help = Height (m) above which rt should be nudged
    1021        z_nudging_rt=0.
    1022        CALL getin('z_nudging_rt',z_nudging_rt)
    1023 
    1024 !
    1025 !Parametres de forcage
    1026 !Config  Key  = z_nudging_u
    1027 !Config  Desc = Height (m) above which u should be nudged
    1028 !Config  Def  = false
    1029 !Config  Help = Height (m) above which u should be nudged
    1030        z_nudging_u=0.
    1031        CALL getin('z_nudging_u',z_nudging_u)
    1032 
    1033 !
    1034 !Parametres de forcage
    1035 !Config  Key  = z_nudging_v
    1036 !config  desc = height (m) above which v should be nudged
    1037 !config  def  = false
    1038 !config  help = height (m) above which v should be nudged
    1039        z_nudging_v=0.
    1040        call getin('z_nudging_v',z_nudging_v)
    1041 
    1042 !config  key  = ok_flux_surf
    1043 !config  desc = forcage ou non par les flux de surface
    1044 !config  def  = false
    1045 !config  help = forcage ou non par les flux de surface
    1046        ok_flux_surf =.false.
    1047        call getin('ok_flux_surf',ok_flux_surf)
    1048 
    1049 !config  key  = ok_prescr_ust
    1050 !config  desc = ustar impose ou non
    1051 !config  def  = false
    1052 !config  help = ustar impose ou non
    1053        ok_prescr_ust = .false.
    1054        call getin('ok_prescr_ust',ok_prescr_ust)
    1055 
    1056 !config  key  = ok_old_disvert
    1057 !config  desc = utilisation de l ancien programme disvert0 (dans 1dutils.h)
    1058 !config  def  = false
    1059 !config  help = utilisation de l ancien programme disvert0 (dans 1dutils.h)
    1060        ok_old_disvert = .false.
    1061        call getin('ok_old_disvert',ok_old_disvert)
    1062 
    1063 !config  key  = time_ini
    1064 !config  desc = meaningless in this  case
    1065 !config  def  = 0.
    1066 !config  help =
    1067        tsurf = 0.
    1068        call getin('time_ini',time_ini)
    1069 
    1070 !config  key  = rlat et rlon
    1071 !config  desc = latitude et longitude
    1072 !config  def  = 0.0  0.0
    1073 !config  help = fixe la position de la colonne
    1074        xlat = 0.
    1075        xlon = 0.
    1076        call getin('rlat',xlat)
    1077        call getin('rlon',xlon)
    1078 
    1079 !config  key  = airephy
    1080 !config  desc = grid cell area
    1081 !config  def  = 1.e11
    1082 !config  help =
    1083        airefi = 1.e11
    1084        call getin('airephy',airefi)
    1085 
    1086 !config  key  = nat_surf
    1087 !config  desc = surface type
    1088 !config  def  = 0 (ocean)
    1089 !config  help = 0=ocean,1=land,2=glacier,3=banquise
    1090        nat_surf = 0.
    1091        call getin('nat_surf',nat_surf)
    1092 
    1093 !config  key  = tsurf
    1094 !config  desc = surface temperature
    1095 !config  def  = 290.
    1096 !config  help = not used if type_ts_forcing=1 in lmdz1d.f
    1097        tsurf = 290.
    1098        call getin('tsurf',tsurf)
    1099 
    1100 !config  key  = psurf
    1101 !config  desc = surface pressure
    1102 !config  def  = 102400.
    1103 !config  help =
    1104        psurf = 102400.
    1105        call getin('psurf',psurf)
    1106 
    1107 !config  key  = zsurf
    1108 !config  desc = surface altitude
    1109 !config  def  = 0.
    1110 !config  help =
    1111        zsurf = 0.
    1112        call getin('zsurf',zsurf)
    1113 
    1114 !config  key  = rugos
    1115 !config  desc = coefficient de frottement
    1116 !config  def  = 0.0001
    1117 !config  help = calcul du cdrag
    1118        rugos = 0.0001
    1119        call getin('rugos',rugos)
    1120 ! fh/2020/04/08/confinement: pour le nouveau format standard, la rugosite s'appelle z0
    1121        call getin('z0',rugos)
    1122 
    1123 !config  key  = rugosh
    1124 !config  desc = coefficient de frottement
    1125 !config  def  = rugos
    1126 !config  help = calcul du cdrag
    1127        rugosh = rugos
    1128        call getin('rugosh',rugosh)
    1129 
    1130 !config  key  = snowmass
    1131 !config  desc = mass de neige de la surface en kg/m2
    1132 !config  def  = 0.0000
    1133 !config  help = snowmass
    1134        snowmass = 0.0000
    1135        call getin('snowmass',snowmass)
    1136 
    1137 !config  key  = wtsurf et wqsurf
    1138 !config  desc = ???
    1139 !config  def  = 0.0 0.0
    1140 !config  help =
    1141        wtsurf = 0.0
    1142        wqsurf = 0.0
    1143        call getin('wtsurf',wtsurf)
    1144        call getin('wqsurf',wqsurf)
    1145 
    1146 !config  key  = albedo
    1147 !config  desc = albedo
    1148 !config  def  = 0.09
    1149 !config  help =
    1150        albedo = 0.09
    1151        call getin('albedo',albedo)
    1152 
    1153 !config  key  = agesno
    1154 !config  desc = age de la neige
    1155 !config  def  = 30.0
    1156 !config  help =
    1157        xagesno = 30.0
    1158        call getin('agesno',xagesno)
    1159 
    1160 !config  key  = restart_runoff
    1161 !config  desc = age de la neige
    1162 !config  def  = 30.0
    1163 !config  help =
    1164        restart_runoff = 0.0
    1165        call getin('restart_runoff',restart_runoff)
    1166 
    1167 !config  key  = qsolinp
    1168 !config  desc = initial bucket water content (kg/m2) when land (5std)
    1169 !config  def  = 30.0
    1170 !config  help =
    1171        qsolinp = 1.
    1172        call getin('qsolinp',qsolinp)
    1173 
    1174 !config  key  = zpicinp
    1175 !config  desc = denivellation orographie
    1176 !config  def  = 0.
    1177 !config  help =  input brise
    1178        zpicinp = 0.
    1179        call getin('zpicinp',zpicinp)
    1180 
    1181 !config key = nudge_tsoil
    1182 !config  desc = activation of soil temperature nudging
    1183 !config  def  = .false.
    1184 !config  help = ...
    1185        nudge_tsoil=.false.
    1186        call getin('nudge_tsoil',nudge_tsoil)
    1187 
    1188 !config key = isoil_nudge
    1189 !config  desc = level number where soil temperature is nudged
    1190 !config  def  = 3
    1191 !config  help = ...
    1192        isoil_nudge=3
    1193        call getin('isoil_nudge',isoil_nudge)
    1194 
    1195 !config key = tsoil_nudge
    1196 !config  desc = target temperature for tsoil(isoil_nudge)
    1197 !config  def  = 300.
    1198 !config  help = ...
    1199        tsoil_nudge=300.
    1200        call getin('tsoil_nudge',tsoil_nudge)
    1201 
    1202 !config key = tau_soil_nudge
    1203 !config  desc = nudging relaxation time for tsoil
    1204 !config  def  = 3600.
    1205 !config  help = ...
    1206        tau_soil_nudge=3600.
    1207        call getin('tau_soil_nudge',tau_soil_nudge)
    1208 
    1209 !config  key  = forc_omega
    1210 !config  desc = forcage ou non par omega
    1211 !config  def  = false
    1212 !config  help = forcage ou non par omega
    1213        forc_omega =0
    1214        call getin('forc_omega',forc_omega)
    1215 
    1216 !config  key  = forc_w
    1217 !config  desc = forcage ou non par w
    1218 !config  def  = false
    1219 !config  help = forcage ou non par w
    1220        forc_w =0
    1221        call getin('forc_w',forc_w)
    1222 
    1223 !config  key  = forc_geo
    1224 !config  desc = forcage ou non par geo
    1225 !config  def  = false
    1226 !config  help = forcage ou non par geo
    1227        forc_geo =0
    1228        call getin('forc_geo',forc_geo)
    1229 
    1230 ! meme chose que ok_precr_ust
    1231 !config  key  = forc_ustar
    1232 !config  desc = forcage ou non par ustar
    1233 !config  def  = false
    1234 !config  help = forcage ou non par ustar
    1235        forc_ustar =0
    1236        call getin('forc_ustar',forc_ustar)
    1237        if (forc_ustar .eq. 1) ok_prescr_ust=.true.
    1238 
    1239 
    1240       write(lunout,*)' +++++++++++++++++++++++++++++++++++++++'
    1241       write(lunout,*)' configuration des parametres du gcm1d: '
    1242       write(lunout,*)' +++++++++++++++++++++++++++++++++++++++'
    1243       write(lunout,*)' restart = ', restart
    1244       write(lunout,*)' forcing_type = ', forcing_type
    1245       write(lunout,*)' time_ini = ', time_ini
    1246       write(lunout,*)' rlat = ', xlat
    1247       write(lunout,*)' rlon = ', xlon
    1248       write(lunout,*)' airephy = ', airefi
    1249       write(lunout,*)' nat_surf = ', nat_surf
    1250       write(lunout,*)' tsurf = ', tsurf
    1251       write(lunout,*)' psurf = ', psurf
    1252       write(lunout,*)' zsurf = ', zsurf
    1253       write(lunout,*)' rugos = ', rugos
    1254       write(lunout,*)' snowmass=', snowmass
    1255       write(lunout,*)' wtsurf = ', wtsurf
    1256       write(lunout,*)' wqsurf = ', wqsurf
    1257       write(lunout,*)' albedo = ', albedo
    1258       write(lunout,*)' xagesno = ', xagesno
    1259       write(lunout,*)' restart_runoff = ', restart_runoff
    1260       write(lunout,*)' qsolinp = ', qsolinp
    1261       write(lunout,*)' zpicinp = ', zpicinp
    1262       write(lunout,*)' nudge_tsoil = ', nudge_tsoil
    1263       write(lunout,*)' isoil_nudge = ', isoil_nudge
    1264       write(lunout,*)' tsoil_nudge = ', tsoil_nudge
    1265       write(lunout,*)' tau_soil_nudge = ', tau_soil_nudge
    1266       write(lunout,*)' adv_temp =      ', adv_temp
    1267       write(lunout,*)' adv_theta =      ', adv_theta
    1268       write(lunout,*)' adv_thetal =      ', adv_thetal
    1269       write(lunout,*)' rad_temp =      ', rad_temp
    1270       write(lunout,*)' rad_theta =      ', rad_theta
    1271       write(lunout,*)' rad_thetal =      ', rad_thetal
    1272       write(lunout,*)' adv_qv =      ', adv_qv
    1273       write(lunout,*)' adv_qt =      ', adv_qt
    1274       write(lunout,*)' adv_rv =      ', adv_rv
    1275       write(lunout,*)' adv_rt =      ', adv_rt
    1276       write(lunout,*)' nudging_temp =      ', nudging_temp
    1277       write(lunout,*)' nudging_theta =      ', nudging_theta
    1278       write(lunout,*)' nudging_thetal =      ', nudging_thetal
    1279       write(lunout,*)' nudging_qv =      ', nudging_qv
    1280       write(lunout,*)' nudging_qt =      ', nudging_qt
    1281       write(lunout,*)' nudging_rv =      ', nudging_rv
    1282       write(lunout,*)' nudging_rt =      ', nudging_rt
    1283       write(lunout,*)' p_nudging_temp =      ', p_nudging_temp
    1284       write(lunout,*)' p_nudging_theta =      ', p_nudging_theta
    1285       write(lunout,*)' p_nudging_thetal =      ', p_nudging_thetal
    1286       write(lunout,*)' p_nudging_qv =      ', p_nudging_qv
    1287       write(lunout,*)' p_nudging_qt =      ', p_nudging_qt
    1288       write(lunout,*)' p_nudging_rv =      ', p_nudging_rv
    1289       write(lunout,*)' p_nudging_rt =      ', p_nudging_rt
    1290       write(lunout,*)' z_nudging_temp =      ', z_nudging_temp
    1291       write(lunout,*)' z_nudging_theta =      ', z_nudging_theta
    1292       write(lunout,*)' z_nudging_thetal =      ', z_nudging_thetal
    1293       write(lunout,*)' z_nudging_qv =      ', z_nudging_qv
    1294       write(lunout,*)' z_nudging_qt =      ', z_nudging_qt
    1295       write(lunout,*)' z_nudging_rv =      ', z_nudging_rv
    1296       write(lunout,*)' z_nudging_rt =      ', z_nudging_rt
    1297       write(lunout,*)' forc_omega = ', forc_omega
    1298       write(lunout,*)' forc_w     = ', forc_w
    1299       write(lunout,*)' forc_geo   = ', forc_geo
    1300       write(lunout,*)' forc_ustar = ', forc_ustar
    1301       IF (forcing_type .eq.40) THEN
    1302         write(lunout,*) '--- Forcing type GCSS Old --- with:'
    1303         write(lunout,*)'imp_fcg',imp_fcg_gcssold
    1304         write(lunout,*)'ts_fcg',ts_fcg_gcssold
    1305         write(lunout,*)'tp_fcg',Tp_fcg_gcssold
    1306         write(lunout,*)'tp_ini',Tp_ini_gcssold
    1307         write(lunout,*)'xturb_fcg',xTurb_fcg_gcssold
    1308       ENDIF
    1309 
    1310       write(lunout,*)' +++++++++++++++++++++++++++++++++++++++'
    1311       write(lunout,*)
    1312 !
    1313       RETURN
    1314       END
    1315629!
    1316630! $Id: dyn1deta0.F 1279 2010/07/30 A Lahellec$
Note: See TracChangeset for help on using the changeset viewer.