Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm
- Timestamp:
- Aug 3, 2024, 2:56:58 PM (5 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm
- Files:
-
- 33 edited
-
gppref.F90 (modified) (1 diff)
-
gpxyb.F90 (modified) (3 diffs)
-
lw.F90 (modified) (3 diffs)
-
lwb.F90 (modified) (1 diff)
-
lwbv.F90 (modified) (1 diff)
-
lwc.F90 (modified) (8 diffs)
-
lwu.F90 (modified) (2 diffs)
-
posnam.F90 (modified) (1 diff)
-
radlsw.F90 (modified) (5 diffs)
-
read_rsun_rrtm.F90 (modified) (1 diff)
-
readaerosolstrato1_rrtm.F90 (modified) (7 diffs)
-
readaerosolstrato2_rrtm.F90 (modified) (6 diffs)
-
rrtm_rtrn1a_140gp.F90 (modified) (1 diff)
-
srtm_cldprop.F90 (modified) (4 diffs)
-
srtm_setcoef.F90 (modified) (1 diff)
-
srtm_spcvrt.F90 (modified) (26 diffs)
-
srtm_spcvrt_mcica.F90 (modified) (26 diffs)
-
srtm_srtm_224gp.F90 (modified) (9 diffs)
-
srtm_srtm_224gp_mcica.F90 (modified) (9 diffs)
-
srtm_vrtqdr.F90 (modified) (5 diffs)
-
su_aerw.F90 (modified) (2 diffs)
-
sucst.F90 (modified) (1 diff)
-
suecrad.F90 (modified) (8 diffs)
-
suinit.F90 (modified) (2 diffs)
-
suovlp.F90 (modified) (1 diff)
-
suphec.F90 (modified) (1 diff)
-
suphmf.F90 (modified) (1 diff)
-
suphy.F90 (modified) (1 diff)
-
sw1s.F90 (modified) (3 diffs)
-
swclr.F90 (modified) (1 diff)
-
swtt1.F90 (modified) (2 diffs)
-
swu.F90 (modified) (1 diff)
-
swuvo3.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gppref.F90
r5159 r5160 134 134 IF (LVERTFE) THEN 135 135 DO JLEV=1,KFLEV 136 ! print*,'GPPREF: LVERTFE KFLEV KSTART KPROF JLEV',LVERTFE,KFLEV,KSTART,KPROF,JLEV136 ! PRINT *,'GPPREF: LVERTFE KFLEV KSTART KPROF JLEV',LVERTFE,KFLEV,KSTART,KPROF,JLEV 137 137 PRESF(KSTART:KPROF,JLEV)=VAF(JLEV)+VBF(JLEV)*PRESH(KSTART:KPROF,KFLEV) 138 138 ENDDO -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gpxyb.F90
r1990 r5160 121 121 ! The first block if is for economy (no do loop start up) and the second 122 122 ! for safety. 123 ! print*,'GPXYB: NDLNPR RHYDR0=',NDLNPR,RHYDR0123 !PRINT *,'GPXYB: NDLNPR RHYDR0=',NDLNPR,RHYDR0 124 124 TOPPRES=0.1 !!!!! A REVOIR (MPL) 29042010 passe de 0 a 0.1 comme ARPEGE 125 125 IF(PRES(KSTART,0) <= TOPPRES)THEN … … 139 139 ! -------------------- 140 140 141 ! print*,'NDLNPR LVERTFE',NDLNPR,LVERTFE141 !PRINT *,'NDLNPR LVERTFE',NDLNPR,LVERTFE 142 142 IF(NDLNPR == 0) THEN 143 143 … … 177 177 & *(PVDELB(JLEV)+PVC(JLEV)*PLNPR(JLON,JLEV)*PRDELP(JLON,& 178 178 & JLEV)) 179 ! print*,'GPXYB JLEV JLON JJ PRES ZPRES PDELP ', JLEV,JLON,JJ,PRES(JLON,JLEV),ZRPRES(JLON,JJ),PDELP(JLON,JLEV)180 ! print*,'GPXYB JLEV JLON JM PRDELP PLNPR ', JLEV,JLON,JM,PRDELP(JLON,JLEV),PLNPR (JLON,JLEV)181 ! print*,'GPXYB JLEV JLON JJ PRPRES PALPH ', JLEV,JLON,JJ,PRPRES(JLON,JLEV),PALPH (JLON,JLEV)182 ! print*,'GPXYB JLEV JLON PRPP PRTGR PVDELB PVC ', JLEV,JLON,PRPP (JLON,JLEV),PRTGR (JLON,JLEV),PVDELB(JLEV),PVC(JLEV)179 ! PRINT *,'GPXYB JLEV JLON JJ PRES ZPRES PDELP ', JLEV,JLON,JJ,PRES(JLON,JLEV),ZRPRES(JLON,JJ),PDELP(JLON,JLEV) 180 ! PRINT *,'GPXYB JLEV JLON JM PRDELP PLNPR ', JLEV,JLON,JM,PRDELP(JLON,JLEV),PLNPR (JLON,JLEV) 181 ! PRINT *,'GPXYB JLEV JLON JJ PRPRES PALPH ', JLEV,JLON,JJ,PRPRES(JLON,JLEV),PALPH (JLON,JLEV) 182 ! PRINT *,'GPXYB JLEV JLON PRPP PRTGR PVDELB PVC ', JLEV,JLON,PRPP (JLON,JLEV),PRTGR (JLON,JLEV),PVDELB(JLEV),PVC(JLEV) 183 183 ENDDO 184 184 JTEMP=JM -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lw.F90
r5158 r5160 142 142 143 143 IF (LHOOK) CALL DR_HOOK('LW',0,ZHOOK_HANDLE) 144 print*,' LW: Avant LWU'144 PRINT *,' LW: Avant LWU' 145 145 CALL LWU & 146 146 & ( KIDIA, KFDIA, KLON, KLEV,& … … 156 156 ! --------------------------------- 157 157 158 print*,' LW: Avant LWBV'158 PRINT *,' LW: Avant LWBV' 159 159 CALL LWBV & 160 160 & ( KIDIA, KFDIA, KLON , KLEV , KMODE,& … … 169 169 ! -------------------------------- 170 170 171 print*,' LW: Avant LWC'171 PRINT *,' LW: Avant LWC' 172 172 CALL LWC & 173 173 & ( KIDIA , KFDIA, KLON , KLEV,& -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwb.F90
r1990 r5160 115 115 ! ------------------------------ 116 116 117 print*,'dans LWB'117 PRINT *,'dans LWB' 118 118 IF (LHOOK) CALL DR_HOOK('LWB',0,ZHOOK_HANDLE) 119 119 ILEV2=2*KLEV -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwbv.F90
r1990 r5160 123 123 124 124 IF (LHOOK) CALL DR_HOOK('LWBV',0,ZHOOK_HANDLE) 125 print*,'LWBV: avant LWB'125 PRINT *,'LWBV: avant LWB' 126 126 CALL LWB & 127 127 & ( KIDIA, KFDIA, KLON , KLEV , KMODE,& -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwc.F90
r1990 r5160 114 114 !100 CONTINUE 115 115 116 ! print*,' Enter LWC '116 ! PRINT *,' Enter LWC ' 117 117 IF (LHOOK) CALL DR_HOOK('LWC',0,ZHOOK_HANDLE) 118 118 DO JL = KIDIA,KFDIA … … 152 152 ENDDO 153 153 ENDDO 154 ! print*,' LWC after Initialisation to clear-sky fluxes'154 ! PRINT *,' LWC after Initialisation to clear-sky fluxes' 155 155 156 156 !* 2.1 FLUXES FOR ONE OVERCAST UNITY EMISSIVITY CLOUD … … 211 211 212 212 ENDDO 213 ! print*,' LWC after 213: Fluxes for unity emissivity'213 ! PRINT *,' LWC after 213: Fluxes for unity emissivity' 214 214 215 215 !* 2.2 CLOUD COVER MATRIX … … 228 228 ENDDO 229 229 ENDDO 230 ! print*,' LWC after Initialisation CC matrix'230 ! PRINT *,' LWC after Initialisation CC matrix' 231 231 232 232 !* 2.4 CLOUD COVER BELOW THE LEVEL OF CALCULATION … … 279 279 280 280 ENDDO 281 ! print*,' LWC after 244: CC below level of calculation'281 ! PRINT *,' LWC after 244: CC below level of calculation' 282 282 283 283 !* 2.5 CLOUD COVER ABOVE THE LEVEL OF CALCULATION … … 329 329 ENDDO 330 330 ENDDO 331 ! print*,' LWC after 254: CC above level of calculation'331 ! PRINT *,' LWC after 254: CC above level of calculation' 332 332 333 333 !* 3. FLUXES FOR PARTIAL/MULTIPLE LAYERED CLOUDINESS … … 371 371 372 372 ENDDO 373 ! print*,' LWC after 317: Downward fluxes'373 ! PRINT *,' LWC after 317: Downward fluxes' 374 374 375 375 !* 3.2 UPWARD FLUX AT THE SURFACE … … 413 413 414 414 ENDDO 415 ! print*,' LWC after 337: Upward fluxes'415 ! PRINT *,' LWC after 337: Upward fluxes' 416 416 417 417 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwu.F90
r5159 r5160 223 223 IAE2 = 3 * KLEV + 1 - (IJ + 1) 224 224 IAE3 = 3 * KLEV + 1 - IJPN 225 ! print*,'IAE1= ',IAE1226 ! print*,'IAE2= ',IAE2227 ! print*,'IAE3= ',IAE3228 ! print*,'KIDIA= ',KIDIA229 ! print*,'KFDIA= ',KFDIA230 ! print*,'KLEV= ',KLEV225 ! PRINT *,'IAE1= ',IAE1 226 ! PRINT *,'IAE2= ',IAE2 227 ! PRINT *,'IAE3= ',IAE3 228 ! PRINT *,'KIDIA= ',KIDIA 229 ! PRINT *,'KFDIA= ',KFDIA 230 ! PRINT *,'KLEV= ',KLEV 231 231 DO JAE = 1, 6 232 232 DO JL = KIDIA, KFDIA 233 ! print*,'JL= ',JL,'-JAE= ',JAE,'-JK= ',JK,'-NSIL= ',NSIL233 ! PRINT *,'JL= ',JL,'-JAE= ',JAE,'-JK= ',JK,'-NSIL= ',NSIL 234 234 ZUAER(JL, JAE) = & 235 235 & (RAER(JAE, 1) * PAER(JL, 1, JK) + RAER(JAE, 2) * PAER(JL, 2, JK)& … … 363 363 364 364 ENDDO 365 ! print*,'END OF LWU'365 ! PRINT *,'END OF LWU' 366 366 367 367 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/posnam.F90
r1990 r5160 71 71 ISTATUS=0 72 72 ISCAN=0 73 print*,'On cherche a lire:',CDNAML73 PRINT *,'On cherche a lire:',CDNAML 74 74 DO WHILE (ISTATUS==0 .AND. ISCAN==0) 75 75 READ(KULNAM,'(A)',IOSTAT=ISTATUS) CLINE 76 ! print*,'CLINE,ISTATUS= ',CLINE,ISTATUS76 ! PRINT *,'CLINE,ISTATUS= ',CLINE,ISTATUS 77 77 SELECT CASE (ISTATUS) 78 78 CASE (:-1) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/radlsw.F90
r5159 r5160 1064 1064 ! ------------------------------------ 1065 1065 1066 ! print*,'RADLSW: LPHYLIN, LRRTM',LPHYLIN, LRRTM1066 !PRINT *,'RADLSW: LPHYLIN, LRRTM',LPHYLIN, LRRTM 1067 1067 IF (.NOT.LPHYLIN) THEN 1068 1068 IF ( .NOT. LRRTM) THEN … … 1076 1076 & ZEMIT , PFLUX , PFLUC & 1077 1077 & ) 1078 ! print*,'RADLSW: apres CALL LW'1078 ! PRINT *,'RADLSW: apres CALL LW' 1079 1079 IF(LLDEBUG) THEN 1080 1080 call writefield_phy('radlsw_flux1',PFLUX(:,1,:),klev+1) … … 1103 1103 ENDDO 1104 1104 1105 ! print*,'RADLSW: avant CALL RRTM_RRTM_140GP,PAP=',PAP(1,:)1105 ! PRINT *,'RADLSW: avant CALL RRTM_RRTM_140GP,PAP=',PAP(1,:) 1106 1106 CALL RRTM_RRTM_140GP & 1107 1107 & ( KIDIA , KFDIA , KLON , KLEV,& … … 1113 1113 & PTAU_LW,& 1114 1114 & ZEMIT , PFLUX , PFLUC , ZTCLEAR ) 1115 ! print*,'RADLSW: apres CALL RRTM_RRTM_140GP'1115 ! PRINT *,'RADLSW: apres CALL RRTM_RRTM_140GP' 1116 1116 1117 1117 ENDIF … … 1120 1120 PFLUX(:,:,:)= 0.0_JPRB 1121 1121 PFLUC(:,:,:)= 0.0_JPRB 1122 ! print*,'RADLSW: ZEMIT,PFLUX et PFLUC = 0'1122 ! PRINT *,'RADLSW: ZEMIT,PFLUX et PFLUC = 0' 1123 1123 ENDIF 1124 1124 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/read_rsun_rrtm.F90
r5159 r5160 102 102 solaire=TSI(days_elapsed+1) 103 103 104 print*,'READ_RSUN_RRTM day=', days_elapsed+1,' solaire=', solaire, ' RSUN=', RSUN(1:NSW)104 PRINT *,'READ_RSUN_RRTM day=', days_elapsed+1,' solaire=', solaire, ' RSUN=', RSUN(1:NSW) 105 105 106 106 ENDIF !--fin allocation -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90
r5159 r5160 85 85 86 86 IF (nbands_sw_rrtm.NE.6) THEN 87 print*,'nbands_sw_rrtm doit etre egal a 6 dans readaerosolstrat_rrtm'87 PRINT *,'nbands_sw_rrtm doit etre egal a 6 dans readaerosolstrat_rrtm' 88 88 STOP 89 89 ENDIF … … 95 95 n_lev = size(lev) 96 96 IF (n_lev.NE.klev) THEN 97 print*,'Le nombre de niveaux n est pas egal a klev'97 PRINT *,'Le nombre de niveaux n est pas egal a klev' 98 98 STOP 99 99 ENDIF … … 102 102 CALL nf95_gw_var(ncid_in, varid, latitude) 103 103 n_lat = size(latitude) 104 print*, 'LAT aerosol strato=', n_lat, latitude104 PRINT *, 'LAT aerosol strato=', n_lat, latitude 105 105 106 106 IF (grid_type/=unstructured) THEN 107 107 IF (n_lat.NE.nbp_lat) THEN 108 print*,'Le nombre de lat n est pas egal a nbp_lat'108 PRINT *,'Le nombre de lat n est pas egal a nbp_lat' 109 109 STOP 110 110 ENDIF … … 114 114 CALL nf95_gw_var(ncid_in, varid, longitude) 115 115 n_lon = size(longitude) 116 print*, 'LON aerosol strato=', n_lon, longitude116 PRINT *, 'LON aerosol strato=', n_lon, longitude 117 117 118 118 IF (grid_type/=unstructured) THEN 119 119 IF (n_lon.NE.nbp_lon) THEN 120 print*,'Le nombre de lon n est pas egal a nbp_lon'120 PRINT *,'Le nombre de lon n est pas egal a nbp_lon' 121 121 STOP 122 122 ENDIF … … 127 127 CALL nf95_gw_var(ncid_in, varid, time) 128 128 n_month = size(time) 129 print*, 'TIME aerosol strato=', n_month, time129 PRINT *, 'TIME aerosol strato=', n_month, time 130 130 IF (n_month.NE.12) THEN 131 print*,'Le nombre de month n est pas egal a 12'131 PRINT *,'Le nombre de month n est pas egal a 12' 132 132 STOP 133 133 ENDIF … … 140 140 CALL nf95_inq_varid(ncid_in, "TAUSTRAT", varid) 141 141 ncerr = nf90_get_var(ncid_in, varid, tauaerstrat) 142 print*,'code erreur readaerosolstrato=', ncerr, varid142 PRINT *,'code erreur readaerosolstrato=', ncerr, varid 143 143 144 144 CALL nf95_close(ncid_in) … … 146 146 !---select the correct month 147 147 IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN 148 print*,'probleme avec le mois dans readaerosolstrat =', mth_cur148 PRINT *,'probleme avec le mois dans readaerosolstrat =', mth_cur 149 149 ENDIF 150 150 tauaerstrat_mois(:,:,:) = tauaerstrat(:,:,:,mth_cur) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90
r5159 r5160 92 92 !--check mth_cur 93 93 IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN 94 print*,'probleme avec le mois dans readaerosolstrat =', mth_cur94 PRINT *,'probleme avec le mois dans readaerosolstrat =', mth_cur 95 95 ENDIF 96 96 … … 116 116 IF (grid_type/=unstructured) THEN 117 117 IF (n_lat.NE.nbp_lat) THEN 118 print*, 'latitude=', n_lat, nbp_lat118 PRINT *, 'latitude=', n_lat, nbp_lat 119 119 abort_message='Le nombre de lat n est pas egal a nbp_lat' 120 120 CALL abort_physic(modname,abort_message,1) … … 133 133 CALL nf95_gw_var(ncid_in, varid, wav) 134 134 n_wav = size(wav) 135 print*, 'WAV aerosol strato=', n_wav, wav135 PRINT *, 'WAV aerosol strato=', n_wav, wav 136 136 IF (n_wav.NE.NSW) THEN 137 137 abort_message='Le nombre de wav n est pas egal a NSW' … … 146 146 CALL nf95_inq_varid(ncid_in, "TAU_SUN", varid) 147 147 ncerr = nf90_get_var(ncid_in, varid, tauaerstrat) 148 print*,'code erreur readaerosolstrato=', ncerr, varid148 PRINT *,'code erreur readaerosolstrato=', ncerr, varid 149 149 150 150 !--reading stratospheric aerosol omega per layer 151 151 CALL nf95_inq_varid(ncid_in, "OME_SUN", varid) 152 152 ncerr = nf90_get_var(ncid_in, varid, pizaerstrat) 153 print*,'code erreur readaerosolstrato=', ncerr, varid153 PRINT *,'code erreur readaerosolstrato=', ncerr, varid 154 154 155 155 !--reading stratospheric aerosol g per layer 156 156 CALL nf95_inq_varid(ncid_in, "GGG_SUN", varid) 157 157 ncerr = nf90_get_var(ncid_in, varid, cgaerstrat) 158 print*,'code erreur readaerosolstrato sw=', ncerr, varid158 PRINT *,'code erreur readaerosolstrato sw=', ncerr, varid 159 159 160 160 CALL nf95_close(ncid_in) … … 224 224 CALL nf95_gw_var(ncid_in, varid, wav) 225 225 n_wav = size(wav) 226 print*, 'WAV aerosol strato=', n_wav, wav226 PRINT *, 'WAV aerosol strato=', n_wav, wav 227 227 IF (n_wav.NE.NLW) THEN 228 228 abort_message='Le nombre de wav n est pas egal a NLW' … … 235 235 CALL nf95_inq_varid(ncid_in, "TAU_EAR", varid) 236 236 ncerr = nf90_get_var(ncid_in, varid, taulwaerstrat) 237 print*,'code erreur readaerosolstrato lw=', ncerr, varid237 PRINT *,'code erreur readaerosolstrato lw=', ncerr, varid 238 238 239 239 CALL nf95_close(ncid_in) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/rrtm_rtrn1a_140gp.F90
r5133 r5160 529 529 DO I_LEV = 1, KLEV 530 530 DO IBAND = K_ISTART, K_IEND 531 ! print*,'RTRN1A: I_LEV JPLAY IBAND INDLAY',I_LEV,JPLAY,IBAND,INDLAY(I_LEV)531 ! PRINT *,'RTRN1A: I_LEV JPLAY IBAND INDLAY',I_LEV,JPLAY,IBAND,INDLAY(I_LEV) 532 532 !---- 533 533 !- Calculate the integrated Planck functions for at the -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_cldprop.F90
r1990 r5160 116 116 117 117 IF (I_NDBUG <= 2) THEN 118 print*,'cldprop before loop K_INFLAG, K_ICEFLAG, K_LIQFLAG:',K_INFLAG,K_ICEFLAG,K_LIQFLAG,IB1,IB2118 PRINT *,'cldprop before loop K_INFLAG, K_ICEFLAG, K_LIQFLAG:',K_INFLAG,K_ICEFLAG,K_LIQFLAG,IB1,IB2 119 119 ENDIF 120 120 … … 216 216 ENDDO 217 217 ENDIF 218 print*,'end of ice computations for I_LAY=',I_LAY218 PRINT *,'end of ice computations for I_LAY=',I_LAY 219 219 220 220 ! Calculation of absorption coefficients due to water clouds. … … 270 270 271 271 IF (I_NDBUG <= 1) THEN 272 print*,'end of liquid water computations for I_LAY=',I_LAY272 PRINT *,'end of liquid water computations for I_LAY=',I_LAY 273 273 ENDIF 274 274 … … 332 332 333 333 IF (I_NDBUG <= 1) THEN 334 print*,'about to leave SRTM_CLDPROP'334 PRINT *,'about to leave SRTM_CLDPROP' 335 335 ENDIF 336 336 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_setcoef.F90
r1990 r5160 94 94 95 95 !IF (NDBUG.LE.3) THEN 96 ! print*,'-------- Computed in SETCOEF --------'96 ! PRINT *,'-------- Computed in SETCOEF --------' 97 97 ! print 8990 98 98 8990 format(18x,' T PFAC00, 01, 10, 11 PCO2MULT MOL & -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_spcvrt.F90
r1990 r5160 203 203 IB1=JPB1 204 204 IB2=JPB2 205 ! print*,'IB1, IB2, KSW, KMOL, KLEV: ', IB1,IB2,KSW,KMOL,KLEV205 !PRINT *,'IB1, IB2, KSW, KMOL, KLEV: ', IB1,IB2,KSW,KMOL,KLEV 206 206 207 207 IW=0 … … 215 215 ZINCF14(IBM)=0.0_JPRB 216 216 217 ! print*,'=== spectral band === JB= ',JB,' ====== i.e. IBM= ',IBM,' with IGT= ',IGT217 ! PRINT *,'=== spectral band === JB= ',JB,' ====== i.e. IBM= ',IBM,' with IGT= ',IGT 218 218 219 219 !-- for each band, computes the gaseous and Rayleigh optical thickness … … 229 229 & ZSFLXZEN, ZTAUG , ZTAUR & 230 230 & ) 231 ! print*,'After SRTM_TAUMOL16'231 ! PRINT *,'After SRTM_TAUMOL16' 232 232 233 233 ELSEIF (JB == 17) THEN … … 240 240 & ZSFLXZEN, ZTAUG , ZTAUR & 241 241 & ) 242 ! print*,'After SRTM_TAUMOL17'242 ! PRINT *,'After SRTM_TAUMOL17' 243 243 244 244 ELSEIF (JB == 18) THEN … … 251 251 & ZSFLXZEN, ZTAUG , ZTAUR & 252 252 & ) 253 ! print*,'After SRTM_TAUMOL18'253 ! PRINT *,'After SRTM_TAUMOL18' 254 254 255 255 ELSEIF (JB == 19) THEN … … 262 262 & ZSFLXZEN, ZTAUG , ZTAUR & 263 263 & ) 264 ! print*,'After SRTM_TAUMOL19'264 ! PRINT *,'After SRTM_TAUMOL19' 265 265 266 266 ELSEIF (JB == 20) THEN … … 273 273 & ZSFLXZEN, ZTAUG , ZTAUR & 274 274 & ) 275 ! print*,'After SRTM_TAUMOL20'275 ! PRINT *,'After SRTM_TAUMOL20' 276 276 277 277 ELSEIF (JB == 21) THEN … … 284 284 & ZSFLXZEN, ZTAUG , ZTAUR & 285 285 & ) 286 ! print*,'After SRTM_TAUMOL21'286 ! PRINT *,'After SRTM_TAUMOL21' 287 287 288 288 ELSEIF (JB == 22) THEN … … 295 295 & ZSFLXZEN, ZTAUG , ZTAUR & 296 296 & ) 297 ! print*,'After SRTM_TAUMOL22'297 ! PRINT *,'After SRTM_TAUMOL22' 298 298 299 299 ELSEIF (JB == 23) THEN … … 306 306 & ZSFLXZEN, ZTAUG , ZTAUR & 307 307 & ) 308 ! print*,'After SRTM_TAUMOL23'308 ! PRINT *,'After SRTM_TAUMOL23' 309 309 310 310 ELSEIF (JB == 24) THEN … … 317 317 & ZSFLXZEN, ZTAUG , ZTAUR & 318 318 & ) 319 ! print*,'After SRTM_TAUMOL24'319 ! PRINT *,'After SRTM_TAUMOL24' 320 320 321 321 ELSEIF (JB == 25) THEN … … 329 329 & ZSFLXZEN, ZTAUG , ZTAUR & 330 330 & ) 331 ! print*,'After SRTM_TAUMOL25'331 ! PRINT *,'After SRTM_TAUMOL25' 332 332 333 333 ELSEIF (JB == 26) THEN … … 341 341 & ZSFLXZEN, ZTAUG , ZTAUR & 342 342 & ) 343 ! print*,'After SRTM_TAUMOL26'343 ! PRINT *,'After SRTM_TAUMOL26' 344 344 345 345 ELSEIF (JB == 27) THEN … … 353 353 & ZSFLXZEN, ZTAUG , ZTAUR & 354 354 & ) 355 ! print*,'After SRTM_TAUMOL27'355 ! PRINT *,'After SRTM_TAUMOL27' 356 356 357 357 ELSEIF (JB == 28) THEN … … 365 365 & ZSFLXZEN, ZTAUG , ZTAUR & 366 366 & ) 367 ! print*,'After SRTM_TAUMOL28'367 ! PRINT *,'After SRTM_TAUMOL28' 368 368 369 369 ELSEIF (JB == 29) THEN … … 376 376 & ZSFLXZEN , ZTAUG , ZTAUR & 377 377 & ) 378 ! print*,'After SRTM_TAUMOL29'378 ! PRINT *,'After SRTM_TAUMOL29' 379 379 380 380 ENDIF 381 381 382 382 ! IF (NDBUG.LE.3) THEN 383 ! print*,'Incident Solar Flux'383 ! PRINT *,'Incident Solar Flux' 384 384 ! PRINT 9010,(ZSFLXZEN(JG),JG=1,16) 385 385 9010 format(1x,'SolFlx ',16F8.4) 386 ! print *,'Optical thickness for molecular absorption for JB= ',JB386 ! PRINT *,'Optical thickness for molecular absorption for JB= ',JB 387 387 ! DO JK=1,KLEV 388 388 ! PRINT 9011,JK,(ZTAUG(JK,JG),JG=1,16) 389 389 9011 format(1x,'TauGas ',I3,16E9.2) 390 390 ! ENDDO 391 ! print *,'Optical thickness for Rayleigh scattering for JB= ',JB391 ! PRINT *,'Optical thickness for Rayleigh scattering for JB= ',JB 392 392 ! DO JK=1,KLEV 393 393 ! PRINT 9012,JK,(ZTAUR(JK,JG),JG=1,16) 394 394 9012 format(1x,'TauRay ',I3,16E9.2) 395 395 ! ENDDO 396 ! print*,'Cloud optical properties for JB= ',JB396 ! PRINT *,'Cloud optical properties for JB= ',JB 397 397 ! DO JK=1,KLEV 398 398 ! PRINT 9013,JK,PFRCL(JK),PTAUC(JK,IBM),POMGC(JK,IBM),PASYC(JK,IBM) … … 405 405 406 406 ! IF (NDBUG.LE.1) THEN 407 ! print*,' === JG= ',JG,' === for JB= ',JB,' with IW, IBM, JPLAY, KLEV=',IW,IBM,JPLAY,KLEV407 ! PRINT *,' === JG= ',JG,' === for JB= ',JB,' with IW, IBM, JPLAY, KLEV=',IW,IBM,JPLAY,KLEV 408 408 ! ENDIF 409 409 … … 460 460 ZRUP(KLEV+1) =PALBP(IBM) 461 461 ZRUPD(KLEV+1)=PALBD(IBM) 462 ! if (NDBUG < 2) print*,'SWSPCTRL after 1 with JB,JG,IBM and IW= ',JB,JG,IBM,IW462 ! if (NDBUG < 2) PRINT *,'SWSPCTRL after 1 with JB,JG,IBM and IW= ',JB,JG,IBM,IW 463 463 464 464 DO JK=1,KLEV … … 511 511 ! end if 512 512 ENDDO 513 ! if (NDBUG < 2) print*,'SWSPCTRL after 2'513 ! if (NDBUG < 2) PRINT *,'SWSPCTRL after 2' 514 514 515 515 !-- Delta scaling for clear-sky / aerosol optical quantities … … 525 525 & LLRTCHK, ZGCC , PRMU0, ZTAUC , ZOMCC ,& 526 526 & ZREFC , ZREFDC, ZTRAC, ZTRADC ) 527 ! if (NDBUG < 2) print*,'SWSPCTR after SWREFTRA for clear-sky'527 ! if (NDBUG < 2) PRINT *,'SWSPCTR after SWREFTRA for clear-sky' 528 528 529 529 !-- Delta scaling for cloudy quantities … … 545 545 546 546 ENDDO 547 ! if (NDBUG < 2) print*,'SWSPCTR after Delta scaling'547 ! if (NDBUG < 2) PRINT *,'SWSPCTR after Delta scaling' 548 548 549 549 CALL SRTM_REFTRA ( KLEV, I_KMODTS ,& 550 550 & LLRTCHK, ZGCO , PRMU0, ZTAUO , ZOMCO ,& 551 551 & ZREFO , ZREFDO, ZTRAO, ZTRADO ) 552 ! if (NDBUG < 2) print*,'SWSPCTR after SWREFTRA for cloudy'552 ! if (NDBUG < 2) PRINT *,'SWSPCTR after SWREFTRA for cloudy' 553 553 554 554 DO JK=1,KLEV … … 582 582 583 583 ENDDO 584 ! if (NDBUG < 2) print*,'SRTM_SPCVRT after combining clear and cloudy'584 ! if (NDBUG < 2) PRINT *,'SRTM_SPCVRT after combining clear and cloudy' 585 585 586 586 !-- vertical quadrature producing clear-sky fluxes 587 587 588 ! print*,'SRTM_SPCVRT after 3 before SRTM_VRTQDR clear'588 ! PRINT *,'SRTM_SPCVRT after 3 before SRTM_VRTQDR clear' 589 589 590 590 CALL SRTM_VRTQDR ( KLEV, IW ,& … … 594 594 595 595 ! IF (NDBUG < 2) THEN 596 ! print *,'SRTM_SPCVRT out of SRTM_VRTQDR for clear IW=',IW596 ! PRINT *,'SRTM_SPCVRT out of SRTM_VRTQDR for clear IW=',IW 597 597 ! DO JK=1,KLEV+1 598 598 ! print 9201,JK,ZCD(JK,IW),ZCU(JK,IW) … … 603 603 !-- vertical quadrature producing cloudy fluxes 604 604 605 ! print*,'SRTM_SPCVRT after 4 before SRTM_VRTQDR cloudy'605 ! PRINT *,'SRTM_SPCVRT after 4 before SRTM_VRTQDR cloudy' 606 606 607 607 CALL SRTM_VRTQDR ( KLEV, IW ,& … … 611 611 612 612 ! IF (NDBUG < 2) THEN 613 ! print*,'SRTM_SPCVRT out of SRTM_VRTQDR for cloudy IW=',IW613 ! PRINT *,'SRTM_SPCVRT out of SRTM_VRTQDR for cloudy IW=',IW 614 614 ! DO JK=1,KLEV+1 615 615 ! print 9202,JK,ZFD(JK,IW),ZFU(JK,IW) … … 655 655 ENDDO 656 656 657 ! if (NDBUG < 2) print*,'SRTM_SPCVRT end of JG=',JG,' for JB=',JB,' i.e. IW=',IW657 ! if (NDBUG < 2) PRINT *,'SRTM_SPCVRT end of JG=',JG,' for JB=',JB,' i.e. IW=',IW 658 658 ENDDO 659 659 !-- end loop on JG 660 660 661 ! print*,' --- JB= ',JB,' with IB1, IB2= ',IB1,IB2661 ! PRINT *,' --- JB= ',JB,' with IB1, IB2= ',IB1,IB2 662 662 ENDDO 663 663 !-- end loop on JB 664 !if (NDBUG < 2) print*,'SRTM_SPCVRT about to come out'664 !if (NDBUG < 2) PRINT *,'SRTM_SPCVRT about to come out' 665 665 666 666 !DO IBM=1,14 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_spcvrt_mcica.F90
r1990 r5160 210 210 IB1=JPB1 211 211 IB2=JPB2 212 ! print*,'IB1, IB2, KSW, KMOL, KLEV: ', IB1,IB2,KSW,KMOL,KLEV212 !PRINT *,'IB1, IB2, KSW, KMOL, KLEV: ', IB1,IB2,KSW,KMOL,KLEV 213 213 214 214 IW=0 … … 222 222 ZINCF14(IBM)=0.0_JPRB 223 223 224 ! print*,'=== spectral band === JB= ',JB,' ====== i.e. IBM= ',IBM,' with IGT= ',IGT224 ! PRINT *,'=== spectral band === JB= ',JB,' ====== i.e. IBM= ',IBM,' with IGT= ',IGT 225 225 226 226 !-- for each band, computes the gaseous and Rayleigh optical thickness … … 236 236 & ZSFLXZEN, ZTAUG , ZTAUR & 237 237 & ) 238 ! print*,'After SRTM_TAUMOL16'238 ! PRINT *,'After SRTM_TAUMOL16' 239 239 240 240 ELSEIF (JB == 17) THEN … … 247 247 & ZSFLXZEN, ZTAUG , ZTAUR & 248 248 & ) 249 ! print*,'After SRTM_TAUMOL17'249 ! PRINT *,'After SRTM_TAUMOL17' 250 250 251 251 ELSEIF (JB == 18) THEN … … 258 258 & ZSFLXZEN, ZTAUG , ZTAUR & 259 259 & ) 260 ! print*,'After SRTM_TAUMOL18'260 ! PRINT *,'After SRTM_TAUMOL18' 261 261 262 262 ELSEIF (JB == 19) THEN … … 269 269 & ZSFLXZEN, ZTAUG , ZTAUR & 270 270 & ) 271 ! print*,'After SRTM_TAUMOL19'271 ! PRINT *,'After SRTM_TAUMOL19' 272 272 273 273 ELSEIF (JB == 20) THEN … … 280 280 & ZSFLXZEN, ZTAUG , ZTAUR & 281 281 & ) 282 ! print*,'After SRTM_TAUMOL20'282 ! PRINT *,'After SRTM_TAUMOL20' 283 283 284 284 ELSEIF (JB == 21) THEN … … 291 291 & ZSFLXZEN, ZTAUG , ZTAUR & 292 292 & ) 293 ! print*,'After SRTM_TAUMOL21'293 ! PRINT *,'After SRTM_TAUMOL21' 294 294 295 295 ELSEIF (JB == 22) THEN … … 302 302 & ZSFLXZEN, ZTAUG , ZTAUR & 303 303 & ) 304 ! print*,'After SRTM_TAUMOL22'304 ! PRINT *,'After SRTM_TAUMOL22' 305 305 306 306 ELSEIF (JB == 23) THEN … … 313 313 & ZSFLXZEN, ZTAUG , ZTAUR & 314 314 & ) 315 ! print*,'After SRTM_TAUMOL23'315 ! PRINT *,'After SRTM_TAUMOL23' 316 316 317 317 ELSEIF (JB == 24) THEN … … 324 324 & ZSFLXZEN, ZTAUG , ZTAUR & 325 325 & ) 326 ! print*,'After SRTM_TAUMOL24'326 ! PRINT *,'After SRTM_TAUMOL24' 327 327 328 328 ELSEIF (JB == 25) THEN … … 336 336 & ZSFLXZEN, ZTAUG , ZTAUR & 337 337 & ) 338 ! print*,'After SRTM_TAUMOL25'338 ! PRINT *,'After SRTM_TAUMOL25' 339 339 340 340 ELSEIF (JB == 26) THEN … … 348 348 & ZSFLXZEN, ZTAUG , ZTAUR & 349 349 & ) 350 ! print*,'After SRTM_TAUMOL26'350 ! PRINT *,'After SRTM_TAUMOL26' 351 351 352 352 ELSEIF (JB == 27) THEN … … 360 360 & ZSFLXZEN, ZTAUG , ZTAUR & 361 361 & ) 362 ! print*,'After SRTM_TAUMOL27'362 ! PRINT *,'After SRTM_TAUMOL27' 363 363 364 364 ELSEIF (JB == 28) THEN … … 372 372 & ZSFLXZEN, ZTAUG , ZTAUR & 373 373 & ) 374 ! print*,'After SRTM_TAUMOL28'374 ! PRINT *,'After SRTM_TAUMOL28' 375 375 376 376 ELSEIF (JB == 29) THEN … … 383 383 & ZSFLXZEN , ZTAUG , ZTAUR & 384 384 & ) 385 ! print*,'After SRTM_TAUMOL29'385 ! PRINT *,'After SRTM_TAUMOL29' 386 386 387 387 ENDIF 388 388 389 389 ! IF (NDBUG.LE.3) THEN 390 ! print*,'Incident Solar Flux'390 ! PRINT *,'Incident Solar Flux' 391 391 ! PRINT 9010,(ZSFLXZEN(JG),JG=1,16) 392 392 9010 format(1x,'SolFlx ',16F8.4) 393 ! print *,'Optical thickness for molecular absorption for JB= ',JB393 ! PRINT *,'Optical thickness for molecular absorption for JB= ',JB 394 394 ! DO JK=1,KLEV 395 395 ! PRINT 9011,JK,(ZTAUG(JK,JG),JG=1,16) 396 396 9011 format(1x,'TauGas ',I3,16E9.2) 397 397 ! ENDDO 398 ! print *,'Optical thickness for Rayleigh scattering for JB= ',JB398 ! PRINT *,'Optical thickness for Rayleigh scattering for JB= ',JB 399 399 ! DO JK=1,KLEV 400 400 ! PRINT 9012,JK,(ZTAUR(JK,JG),JG=1,16) … … 407 407 408 408 ! IF (NDBUG.LE.1) THEN 409 ! print*,' === JG= ',JG,' === for JB= ',JB,' with IW, IBM, JPLAY, KLEV=',IW,IBM,JPLAY,KLEV409 ! PRINT *,' === JG= ',JG,' === for JB= ',JB,' with IW, IBM, JPLAY, KLEV=',IW,IBM,JPLAY,KLEV 410 410 ! ENDIF 411 411 ! IF (NDBUG.LE.3) THEN 412 ! print*,'Cloud optical properties for JB= ',JB412 ! PRINT *,'Cloud optical properties for JB= ',JB 413 413 ! DO JK=1,KLEV 414 414 ! PRINT 9013,JK,PFRCL(IW,JK),PTAUC(JK,IW),POMGC(JK,IW),PASYC(JK,IW) … … 469 469 ZRUP(KLEV+1) =PALBP(IBM) 470 470 ZRUPD(KLEV+1)=PALBD(IBM) 471 ! if (NDBUG < 2) print*,'SWSPCTRL after 1 with JB,JG,IBM and IW= ',JB,JG,IBM,IW471 ! if (NDBUG < 2) PRINT *,'SWSPCTRL after 1 with JB,JG,IBM and IW= ',JB,JG,IBM,IW 472 472 473 473 … … 523 523 ! end if 524 524 ENDDO 525 ! if (NDBUG < 2) print*,'SWSPCTRL after 2'525 ! if (NDBUG < 2) PRINT *,'SWSPCTRL after 2' 526 526 527 527 !-- Delta scaling for clear-sky / aerosol optical quantities … … 537 537 & LLRTCHK, ZGCC , PRMU0, ZTAUC , ZOMCC ,& 538 538 & ZREFC , ZREFDC, ZTRAC, ZTRADC ) 539 ! if (NDBUG < 2) print*,'SWSPCTR after SWREFTRA for clear-sky'539 ! if (NDBUG < 2) PRINT *,'SWSPCTR after SWREFTRA for clear-sky' 540 540 541 541 !-- Delta scaling for cloudy quantities … … 557 557 558 558 ENDDO 559 ! if (NDBUG < 2) print*,'SWSPCTR after Delta scaling'559 ! if (NDBUG < 2) PRINT *,'SWSPCTR after Delta scaling' 560 560 561 561 CALL SRTM_REFTRA ( KLEV, I_KMODTS ,& 562 562 & LLRTCHK, ZGCO , PRMU0, ZTAUO , ZOMCO ,& 563 563 & ZREFO , ZREFDO, ZTRAO, ZTRADO ) 564 ! if (NDBUG < 2) print*,'SWSPCTR after SWREFTRA for cloudy'564 ! if (NDBUG < 2) PRINT *,'SWSPCTR after SWREFTRA for cloudy' 565 565 566 566 DO JK=1,KLEV … … 606 606 607 607 ENDDO 608 ! if (NDBUG < 2) print*,'SRTM_SPCVRT after combining clear and cloudy'608 ! if (NDBUG < 2) PRINT *,'SRTM_SPCVRT after combining clear and cloudy' 609 609 610 610 !-- vertical quadrature producing clear-sky fluxes 611 611 612 ! print*,'SRTM_SPCVRT after 3 before SRTM_VRTQDR clear'612 ! PRINT *,'SRTM_SPCVRT after 3 before SRTM_VRTQDR clear' 613 613 614 614 CALL SRTM_VRTQDR ( KLEV, IW ,& … … 618 618 619 619 ! IF (NDBUG < 2) THEN 620 ! print *,'SRTM_SPCVRT out of SRTM_VRTQDR for clear IW=',IW620 ! PRINT *,'SRTM_SPCVRT out of SRTM_VRTQDR for clear IW=',IW 621 621 ! DO JK=1,KLEV+1 622 622 ! print 9201,JK,ZCD(JK,IW),ZCU(JK,IW) … … 627 627 !-- vertical quadrature producing cloudy fluxes 628 628 629 ! print*,'SRTM_SPCVRT after 4 before SRTM_VRTQDR cloudy'629 ! PRINT *,'SRTM_SPCVRT after 4 before SRTM_VRTQDR cloudy' 630 630 631 631 CALL SRTM_VRTQDR ( KLEV, IW ,& … … 635 635 636 636 ! IF (NDBUG < 2) THEN 637 ! print*,'SRTM_SPCVRT out of SRTM_VRTQDR for cloudy IW=',IW637 ! PRINT *,'SRTM_SPCVRT out of SRTM_VRTQDR for cloudy IW=',IW 638 638 ! DO JK=1,KLEV+1 639 639 ! print 9202,JK,ZFD(JK,IW),ZFU(JK,IW) … … 679 679 ENDDO 680 680 681 ! if (NDBUG < 2) print*,'SRTM_SPCVRT end of JG=',JG,' for JB=',JB,' i.e. IW=',IW681 ! if (NDBUG < 2) PRINT *,'SRTM_SPCVRT end of JG=',JG,' for JB=',JB,' i.e. IW=',IW 682 682 ENDDO 683 683 !-- end loop on JG 684 684 685 ! print*,' --- JB= ',JB,' with IB1, IB2= ',IB1,IB2685 ! PRINT *,' --- JB= ',JB,' with IB1, IB2= ',IB1,IB2 686 686 ENDDO 687 687 !-- end loop on JB 688 !if (NDBUG < 2) print*,'SRTM_SPCVRT about to come out'689 ! print*,'SRTM_SPCVRT about to come out'688 !if (NDBUG < 2) PRINT *,'SRTM_SPCVRT about to come out' 689 !PRINT *,'SRTM_SPCVRT about to come out' 690 690 691 691 !DO IBM=1,14 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_srtm_224gp.F90
r5159 r5160 133 133 IOVLP = 3 134 134 135 ! print*,'Entering srtm_srtm_224gp'135 !PRINT *,'Entering srtm_srtm_224gp' 136 136 137 137 ICLDATM = 1 … … 148 148 !- coefficients related to the cloud optical properties (original RRTM_SW) 149 149 150 ! print*,'just before SRTM_CLDPROP'150 ! PRINT *,'just before SRTM_CLDPROP' 151 151 152 152 ! DO JK=1,KLEV … … 218 218 ENDDO 219 219 220 ! print*,'ZTOTCC ZCLEAR : ',ZTOTCC,' ',ZCLEAR220 ! PRINT *,'ZTOTCC ZCLEAR : ',ZTOTCC,' ',ZCLEAR 221 221 222 222 DO IMOL = 1, I_NMOL … … 236 236 ! ENDIF 237 237 238 ! print*,'just before SRTM_SETCOEF'238 ! PRINT *,'just before SRTM_SETCOEF' 239 239 240 240 ZFRCL(1:KLEV) = PFRCL(JL, 1:KLEV) … … 253 253 &) 254 254 255 ! print*,'just after SRTM_SETCOEF'255 ! PRINT *,'just after SRTM_SETCOEF' 256 256 257 257 !- call the radiation transfer routine … … 271 271 !- mixing of aerosols 272 272 273 ! print*,'Aerosol optical properties computations'273 ! PRINT *,'Aerosol optical properties computations' 274 274 ! DO JSW=1,KSW 275 275 ! print 9012,JSW,(JAE,RSRTAUA(JSW,JAE),RSRPIZA(JSW,JAE),RSRASYA(JSW,JAE),JAE=1,6) … … 335 335 ENDDO 336 336 337 ! print*,'just before calling STRM_SPCVRT for JL=',JL,' and ZRMU0=',ZRMU0337 ! PRINT *,'just before calling STRM_SPCVRT for JL=',JL,' and ZRMU0=',ZRMU0 338 338 339 339 CALL SRTM_SPCVRT & … … 351 351 &) 352 352 353 ! print*,'SRTM_SRTM_224GP before potential scaling'353 ! PRINT *,'SRTM_SRTM_224GP before potential scaling' 354 354 ! IF (IOVLP == 3) THEN 355 355 ! DO JK=1,KLEV+1 … … 362 362 ! ENDDO 363 363 ! ELSE 364 ! print*,'SRTM_SRTM_224GP after potential scaling'364 ! PRINT *,'SRTM_SRTM_224GP after potential scaling' 365 365 DO JK = 1, KLEV + 1 366 366 PFSUC(JL, 1, JK) = ZADJI0 * ZBBCU(JK) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_srtm_224gp_mcica.F90
r5154 r5160 144 144 IOVLP = 3 145 145 146 ! print*,'Entering srtm_srtm_224gp_mcica'146 !PRINT *,'Entering srtm_srtm_224gp_mcica' 147 147 148 148 ICLDATM = 1 … … 159 159 !- coefficients related to the cloud optical properties (original RRTM_SW) 160 160 161 ! print*,'just before SRTM_CLDPROP'161 ! PRINT *,'just before SRTM_CLDPROP' 162 162 163 163 ! DO JK=1,KLEV … … 214 214 ENDDO 215 215 216 ! print*,'ZTOTCC ZCLEAR : ',ZTOTCC,' ',ZCLEAR216 ! PRINT *,'ZTOTCC ZCLEAR : ',ZTOTCC,' ',ZCLEAR 217 217 218 218 DO IMOL = 1, ITMOL … … 222 222 ENDDO 223 223 224 ! print*,'just before SRTM_SETCOEF'224 ! PRINT *,'just before SRTM_SETCOEF' 225 225 226 226 CALL SRTM_SETCOEF & … … 235 235 &) 236 236 237 ! print*,'just after SRTM_SETCOEF'237 ! PRINT *,'just after SRTM_SETCOEF' 238 238 239 239 !- call the radiation transfer routine … … 265 265 !- mixing of aerosols 266 266 267 ! print*,'Aerosol optical properties computations'267 ! PRINT *,'Aerosol optical properties computations' 268 268 ! DO JSW=1,KSW 269 269 ! print 9012,JSW,(JAE,RSRTAUA(JSW,JAE),RSRPIZA(JSW,JAE),RSRASYA(JSW,JAE),JAE=1,6) … … 329 329 ENDDO 330 330 331 ! print*,'just before calling STRM_SPCVRT for JL=',JL,' and ZRMU0=',ZRMU0331 ! PRINT *,'just before calling STRM_SPCVRT for JL=',JL,' and ZRMU0=',ZRMU0 332 332 333 333 CALL SRTM_SPCVRT_MCICA & … … 347 347 ! & ) 348 348 349 ! print*,'SRTM_SRTM_224GP before potential scaling'349 ! PRINT *,'SRTM_SRTM_224GP before potential scaling' 350 350 ! IF (IOVLP == 3) THEN 351 351 ! DO JK=1,KLEV+1 … … 358 358 ! ENDDO 359 359 ! ELSE 360 ! print*,'SRTM_SRTM_224GP after potential scaling'360 ! PRINT *,'SRTM_SRTM_224GP after potential scaling' 361 361 DO JK = 1, KLEV + 1 362 362 PFSUC(JL, 1, JK) = ZADJI0 * ZBBCU(JK) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_vrtqdr.F90
r1990 r5160 114 114 ! print 9201,PRUP(KLEV),PRUPD(KLEV) 115 115 9201 format(1x,'link surf:',6E13.6) 116 ! print*,'SRTM_VRTQDR after linking with surface layer'116 ! PRINT *,'SRTM_VRTQDR after linking with surface layer' 117 117 !END IF 118 118 … … 134 134 9203 format(1x,'bot2top:',6E13.6) 135 135 ENDDO 136 ! print*,'SRTM_VRTQDR after passing from bottom to top'136 !PRINT *,'SRTM_VRTQDR after passing from bottom to top' 137 137 138 138 !-- upper boundary conditions … … 146 146 ! print 9204,ZTDN(1),PRDND(1),ZTDN(2),PRDND(2) 147 147 9204 format(1x,'link upper bound:',6E13.6) 148 ! print*,'SRTM_VRTQDR after upper boundary conditions'148 ! PRINT *,'SRTM_VRTQDR after upper boundary conditions' 149 149 !END IF 150 150 … … 166 166 167 167 ENDDO 168 ! print*,'SRTM_VRTQDR after passing from top to bottom'168 !PRINT *,'SRTM_VRTQDR after passing from top to bottom' 169 169 170 170 !-- up and down-welling fluxes at levels … … 190 190 191 191 ENDDO 192 ! print*,'SRTM_VRTQDR after up and down flux'192 !PRINT *,'SRTM_VRTQDR after up and down flux' 193 193 194 ! print*,'SRTM_VRTQDR about to come out'194 !PRINT *,'SRTM_VRTQDR about to come out' 195 195 ! ------------------------------------------------------------------ 196 196 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerw.F90
r5158 r5160 90 90 ! ---------------------------- 91 91 92 print*,'DANS SU_AERW'92 PRINT *,'DANS SU_AERW' 93 93 NBINAER(:) = (/ 3, 3, 2, 2, 1, 1, 1, 1, 1 /) 94 94 … … 166 166 167 167 CALL SU_AERP 168 print*,'SU_AERW: apres SU_AERP'168 PRINT *,'SU_AERW: apres SU_AERP' 169 169 CALL SU_AEROP 170 print*,'SU_AERW: apres SU_AEROP'170 PRINT *,'SU_AERW: apres SU_AEROP' 171 171 172 172 IF (LEPAERO) THEN -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sucst.F90
r5158 r5160 202 202 RBETW=RLVTT/RV+RGAMW*RTT 203 203 RALPW=LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT) 204 print*,'SUCST: RESTT,RBETW,RTT,RGAMW',RESTT,RBETW,RTT,RGAMW205 print*,'SUCST: RALPW',RALPW204 PRINT *,'SUCST: RESTT,RBETW,RTT,RGAMW',RESTT,RBETW,RTT,RGAMW 205 PRINT *,'SUCST: RALPW',RALPW 206 206 RGAMS=(RCS-RCPV)/RV 207 207 RBETS=RLSTT/RV+RGAMS*RTT 208 208 RALPS=LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT) 209 print*,'SUCST: RESTT,RBETS,RTT,RGAMS',RESTT,RBETS,RTT,RGAMS210 print*,'SUCST: RALPS',RALPS209 PRINT *,'SUCST: RESTT,RBETS,RTT,RGAMS',RESTT,RBETS,RTT,RGAMS 210 PRINT *,'SUCST: RALPS',RALPS 211 211 RGAMS=(RCS-RCPV)/RV 212 212 RGAMD=RGAMS-RGAMW -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suecrad.F90
r5159 r5160 295 295 NRADIP = 3 ! before 3?R1 default=2 3 296 296 NRADLP = 2 ! before 3?R1 default=2 2 297 print*, 'SUECRAD: NRADLP, NRADIP=', NRADLP, NRADIP297 PRINT *, 'SUECRAD: NRADLP, NRADIP=', NRADLP, NRADIP 298 298 RRe2De = 0.64952_JPRB ! before 3?R1 default=0.5_JPRB 299 299 … … 345 345 !- interaction radiation / prognostic O3 off by default 346 346 LEPO3RA = .FALSE. 347 print*, 'SUECRAD-0'347 PRINT *, 'SUECRAD-0' 348 348 IF (.NOT.YO3%LGP) THEN 349 349 LEPO3RA = .FALSE. … … 407 407 NOVLP = 1 408 408 END SELECT 409 print*, 'SUECRAD: NOVLP=', NOVLP409 PRINT *, 'SUECRAD: NOVLP=', NOVLP 410 410 NLW = 16 411 411 NTSW = 14 … … 462 462 RCCFC11 = CFC11_ppt * 1.0e-12 463 463 RCCFC12 = CFC12_ppt * 1.0e-12 464 ! print*,'LMDZSUECRAD-1 RCCO2=',RCCO2465 ! print*,'LMDZSUECRAD-1 RCCH4=',RCCH4466 ! print*,'LMDZSUECRAD-1 RCN2O=',RCN2O467 ! print*,'LMDZSUECRAD-1 RCCFC11=',RCCFC11468 ! print*,'LMDZSUECRAD-1 RCCFC12=',RCCFC12464 !PRINT *,'LMDZSUECRAD-1 RCCO2=',RCCO2 465 !PRINT *,'LMDZSUECRAD-1 RCCH4=',RCCH4 466 !PRINT *,'LMDZSUECRAD-1 RCN2O=',RCN2O 467 !PRINT *,'LMDZSUECRAD-1 RCCFC11=',RCCFC11 468 !PRINT *,'LMDZSUECRAD-1 RCCFC12=',RCCFC12 469 469 ! ------------------------------------------------------------------ 470 470 … … 474 474 !CALL POSNAM(NULNAM,'NAERAD') 475 475 !READ (NULNAM,NAERAD) 476 print*, 'SUECRAD-2'476 PRINT *, 'SUECRAD-2' 477 477 478 478 !CALL POSNAM(NULNAM,'NAEAER') … … 510 510 RSWINHF = 0.7_JPRB 511 511 ENDIF 512 print*, 'SUECRAD-3'512 PRINT *, 'SUECRAD-3' 513 513 514 514 !- for McICA computations, make sure these parameters are as follows ... … … 520 520 CALL SU_McICA 521 521 ENDIF 522 print*, 'SUECRAD-4'522 PRINT *, 'SUECRAD-4' 523 523 524 524 IF(LLDEBUG)THEN … … 613 613 ENDIF 614 614 ENDIF 615 print*, 'SUECRAD-5'615 PRINT *, 'SUECRAD-5' 616 616 617 617 ! test if radiation grid resolution has been set -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suinit.F90
r3435 r5160 94 94 ! NDLNPR=1: formulation of delta used in non hydrostatic model, 95 95 NDLNPR=0 96 print*,'SUINIT: RHYDR0 NDLNPR',RHYDR0,NDLNPR96 PRINT *,'SUINIT: RHYDR0 NDLNPR',RHYDR0,NDLNPR 97 97 98 98 !---------------------------------------------------------------- … … 114 114 NSSSSS=0 ! LMDZ demarre tjrs a 00h -- MPL 15.04.09 115 115 CALL SUCST(6,NINDAT,NSSSSS,1) 116 print*,'SUINIT: NINDAT, NSSSSS',NINDAT, NSSSSS116 PRINT *,'SUINIT: NINDAT, NSSSSS',NINDAT, NSSSSS 117 117 118 118 IF (LLDEBUG) THEN -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suovlp.F90
r1990 r5160 43 43 DO JK=1,KLEV 44 44 RA1OVLP(JK)=RAOVLP*STZ(JK)+RBOVLP 45 print*,'SU_OVLP: JK RAOVLP STZ RBOVLP:',JK,RAOVLP,STZ(JK),RBOVLP45 PRINT *,'SU_OVLP: JK RAOVLP STZ RBOVLP:',JK,RAOVLP,STZ(JK),RBOVLP 46 46 ENDDO 47 47 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphec.F90
r5159 r5160 190 190 ! VAH(JLEV)=ap(JLEV+1)ap(JLEV+1) 191 191 ! VBH(JLEV)=bp(JLEV+1) 192 ! print*,'SUPHEC: jlev ap bp',JLEV,ap(JLEV+1),bp(JLEV+1)192 ! PRINT *,'SUPHEC: jlev ap bp',JLEV,ap(JLEV+1),bp(JLEV+1) 193 193 VAH(JLEV)=ap(NFLEVG+1-JLEV) 194 194 VBH(JLEV)=bp(NFLEVG+1-JLEV) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphmf.F90
r1990 r5160 93 93 94 94 IF (LHOOK) CALL DR_HOOK('SUPHMF',0,ZHOOK_HANDLE) 95 print*,'SUPHMF: avant SU0PHY'95 PRINT *,'SUPHMF: avant SU0PHY' 96 96 CALL SU0PHY(KULOUT) 97 print*,'SUPHMF: avant SUPHY0'97 PRINT *,'SUPHMF: avant SUPHY0' 98 98 CALL SUPHY0(KULOUT) 99 print*,'SUPHMF: avant SUPHY1'99 PRINT *,'SUPHMF: avant SUPHY1' 100 100 CALL SUPHY1(KULOUT) 101 print*,'SUPHMF: avant SUPHY2'101 PRINT *,'SUPHMF: avant SUPHY2' 102 102 CALL SUPHY2(KULOUT) 103 print*,'SUPHMF: avant SUPHY3'103 PRINT *,'SUPHMF: avant SUPHY3' 104 104 CALL SUPHY3(KULOUT) 105 print*,'SUPHMF: avant SUTOPH'105 PRINT *,'SUPHMF: avant SUTOPH' 106 106 CALL SUTOPH(KULOUT) 107 print*,'SUPHMF: avant VAL923'107 PRINT *,'SUPHMF: avant VAL923' 108 108 109 109 CALL VAL923(LSOLV) 110 110 111 print*,'SUPHMF: avant SUCAPE'111 PRINT *,'SUPHMF: avant SUCAPE' 112 112 CALL SUCAPE(KULOUT) 113 113 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy.F90
r5159 r5160 97 97 ! -------------------- 98 98 99 print*,'---- SUPHY: avant SUPHMF'99 PRINT *,'---- SUPHY: avant SUPHMF' 100 100 CALL SUPHMF(KULOUT) 101 101 102 print*,'---- SUPHY: avant SUGFL'102 PRINT *,'---- SUPHY: avant SUGFL' 103 103 !SUGFL: Set up unified_treatment grid-point fields 104 104 CALL SUGFL -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sw1s.F90
r5154 r5160 191 191 ! IKL = KLEV+1-JK 192 192 ! DO JL = KIDIA,KFDIA 193 ! print*,'Apres SWCLR,SWR RMU0 RMUE ',ZRMU0(JL,IKL),ZRMUE(JL,IKL)193 ! PRINT *,'Apres SWCLR,SWR RMU0 RMUE ',ZRMU0(JL,IKL),ZRMUE(JL,IKL) 194 194 ! ENDDO 195 195 ! ENDDO … … 303 303 304 304 ELSEIF (NSW == 6) THEN 305 ! print*,'... dans SW1S: NSW=',NSW305 !PRINT *,'... dans SW1S: NSW=',NSW 306 306 307 307 !* 3.2 SIX SPECTRAL INTERVALS … … 425 425 PCU(JL,JK) = ZDIRF(JL) * RSUN(KNU) 426 426 !WRITE(*,'("---> Dans SW1S:")') 427 ! print*,'===JL= ',jl427 !PRINT *,'===JL= ',jl 428 428 !WRITE(*,'("ZR1",10E12.5)') (ZR(JL,1)) 429 429 !WRITE(*,'("ZR2",10E12.5)') (ZR(JL,2)) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/swclr.F90
r5154 r5160 229 229 ! MPLFH : ZTRAY N'EST PAS INITIALISE !!!!! A REVOIR (MPL) 230 230 ZTRAY= PRAYL(JL) * PDSIG(JL,JK) 231 ! print*,'>>>>>>> swclr: ZTRAY ',ZTRAY231 ! PRINT *,'>>>>>>> swclr: ZTRAY ',ZTRAY 232 232 ZDENB = ZTRAY + PTAUAZ(JL,JK)*(1.0_JPRB-PPIZAZ(JL,JK)*ZFF) 233 233 ZRATIO=ZTRAY/ZDENB -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/swtt1.F90
r1990 r5160 91 91 DO JA = 1,KABS 92 92 IA=KIND(JA) 93 ! print*,'SWTT1: KNU', KNU93 ! PRINT *,'SWTT1: KNU', KNU 94 94 DO JL = KIDIA,KFDIA 95 95 ZU(JL) = PU(JL,JA) … … 98 98 & * ( APAD(KNU,IA,5) + ZU(JL) * (APAD(KNU,IA,6) + ZU(JL)& 99 99 & * ( APAD(KNU,IA,7) )))))) 100 ! print*,'SWTT1 ZU APAD',IA,ZU(JL),APAD(KNU,IA,1),APAD(KNU,IA,2),&100 ! PRINT *,'SWTT1 ZU APAD',IA,ZU(JL),APAD(KNU,IA,1),APAD(KNU,IA,2),& 101 101 ! &APAD(KNU,IA,3),APAD(KNU,IA,4),APAD(KNU,IA,5),APAD(KNU,IA,6),APAD(KNU,IA,7) 102 102 -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/swu.F90
r5154 r5160 180 180 ZO175(JL) = ZN175(JL) 181 181 ZO190(JL) = ZN190(JL) 182 ! print*,'SWU: RTH2O RTDH2O RTUMG RTDUMG',RTH2O,RTDH2O,RTUMG,RTDUMG183 ! print*,'SWU: RPNH ZDSH2O ZWH2O ZRTH',RPNH,ZDSH2O,ZWH2O,ZRTH184 ! print*,'SWU: RPNU ZDSCO2 PCARDI ZRTU',RPNU,ZDSCO2,PCARDI,ZRTU182 !PRINT *,'SWU: RTH2O RTDH2O RTUMG RTDUMG',RTH2O,RTDH2O,RTUMG,RTDUMG 183 !PRINT *,'SWU: RPNH ZDSH2O ZWH2O ZRTH',RPNH,ZDSH2O,ZWH2O,ZRTH 184 !PRINT *,'SWU: RPNU ZDSCO2 PCARDI ZRTU',RPNU,ZDSCO2,PCARDI,ZRTU 185 185 186 186 !++MODIFCODE -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/swuvo3.F90
r5133 r5160 90 90 LLDEBUG=.FALSE. 91 91 92 ! print*,'Dans SWUVO3, N_VMASS= ',N_VMASS92 !PRINT *,'Dans SWUVO3, N_VMASS= ',N_VMASS 93 93 IF(N_VMASS > 0) THEN 94 94 JLEN=KFDIA-KIDIA+N_VMASS-MOD(KFDIA-KIDIA,N_VMASS)
Note: See TracChangeset
for help on using the changeset viewer.
