Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_spcvrt.F90

Timestamp:

Aug 3, 2024, 2:56:58 PM (7 weeks ago)

Author:

abarral

Message:

Put .h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_spcvrt.F90 (modified) (26 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_spcvrt.F90

@@ -203 +203 @@
 IB1=JPB1
 IB2=JPB2
-!print *,'IB1, IB2, KSW, KMOL, KLEV: ', IB1,IB2,KSW,KMOL,KLEV
+!PRINT *,'IB1, IB2, KSW, KMOL, KLEV: ', IB1,IB2,KSW,KMOL,KLEV
 
 IW=0
@@ -215 +215 @@
   ZINCF14(IBM)=0.0_JPRB
 
-!  print *,'=== spectral band === JB= ',JB,' ====== i.e. IBM= ',IBM,' with IGT= ',IGT
+!  PRINT *,'=== spectral band === JB= ',JB,' ====== i.e. IBM= ',IBM,' with IGT= ',IGT
         
 !-- for each band, computes the gaseous and Rayleigh optical thickness 
@@ -229 +229 @@
      &   ZSFLXZEN, ZTAUG    , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL16'
+!    PRINT *,'After  SRTM_TAUMOL16'
 
   ELSEIF (JB == 17) THEN
@@ -240 +240 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL17'
+!    PRINT *,'After  SRTM_TAUMOL17'
 
   ELSEIF (JB == 18) THEN
@@ -251 +251 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL18'
+!    PRINT *,'After  SRTM_TAUMOL18'
 
   ELSEIF (JB == 19) THEN
@@ -262 +262 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL19'
+!    PRINT *,'After  SRTM_TAUMOL19'
 
   ELSEIF (JB == 20) THEN
@@ -273 +273 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL20'
+!    PRINT *,'After  SRTM_TAUMOL20'
 
   ELSEIF (JB == 21) THEN
@@ -284 +284 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL21'
+!    PRINT *,'After  SRTM_TAUMOL21'
 
   ELSEIF (JB == 22) THEN
@@ -295 +295 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL22'
+!    PRINT *,'After  SRTM_TAUMOL22'
 
   ELSEIF (JB == 23) THEN
@@ -306 +306 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL23'
+!    PRINT *,'After  SRTM_TAUMOL23'
 
   ELSEIF (JB == 24) THEN
@@ -317 +317 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL24'
+!    PRINT *,'After  SRTM_TAUMOL24'
 
   ELSEIF (JB == 25) THEN
@@ -329 +329 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR   &
      & )  
-!    print *,'After  SRTM_TAUMOL25'
+!    PRINT *,'After  SRTM_TAUMOL25'
 
   ELSEIF (JB == 26) THEN
@@ -341 +341 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL26'
+!    PRINT *,'After  SRTM_TAUMOL26'
 
   ELSEIF (JB == 27) THEN
@@ -353 +353 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL27'
+!    PRINT *,'After  SRTM_TAUMOL27'
 
   ELSEIF (JB == 28) THEN
@@ -365 +365 @@
      &   ZSFLXZEN, ZTAUG   , ZTAUR  &
      & )  
-!    print *,'After  SRTM_TAUMOL28'
+!    PRINT *,'After  SRTM_TAUMOL28'
 
   ELSEIF (JB == 29) THEN
@@ -376 +376 @@
      &   ZSFLXZEN , ZTAUG   , ZTAUR    &
      & )  
-!    print *,'After  SRTM_TAUMOL29'
+!    PRINT *,'After  SRTM_TAUMOL29'
 
   ENDIF
 
 !  IF (NDBUG.LE.3) THEN
-!    print *,'Incident Solar Flux'
+!    PRINT *,'Incident Solar Flux'
 !    PRINT 9010,(ZSFLXZEN(JG),JG=1,16)
   9010 format(1x,'SolFlx ',16F8.4)
-!    print *,'Optical thickness for molecular absorption for JB= ',JB 
+!    PRINT *,'Optical thickness for molecular absorption for JB= ',JB
 !    DO JK=1,KLEV
 !      PRINT 9011,JK,(ZTAUG(JK,JG),JG=1,16)
   9011  format(1x,'TauGas ',I3,16E9.2)
 !    ENDDO
-!    print *,'Optical thickness for Rayleigh scattering for JB= ',JB 
+!    PRINT *,'Optical thickness for Rayleigh scattering for JB= ',JB
 !    DO JK=1,KLEV
 !      PRINT 9012,JK,(ZTAUR(JK,JG),JG=1,16)
   9012  format(1x,'TauRay ',I3,16E9.2)
 !    ENDDO
-!    print *,'Cloud optical properties for JB= ',JB
+!    PRINT *,'Cloud optical properties for JB= ',JB
 !    DO JK=1,KLEV
 !      PRINT 9013,JK,PFRCL(JK),PTAUC(JK,IBM),POMGC(JK,IBM),PASYC(JK,IBM)
@@ -405 +405 @@
 
 !    IF (NDBUG.LE.1) THEN
-!      print *,' === JG= ',JG,' === for JB= ',JB,' with IW, IBM, JPLAY, KLEV=',IW,IBM,JPLAY,KLEV
+!      PRINT *,' === JG= ',JG,' === for JB= ',JB,' with IW, IBM, JPLAY, KLEV=',IW,IBM,JPLAY,KLEV
 !    ENDIF
 
@@ -460 +460 @@
     ZRUP(KLEV+1) =PALBP(IBM)
     ZRUPD(KLEV+1)=PALBD(IBM)
-!    if (NDBUG < 2) print *,'SWSPCTRL after 1 with JB,JG,IBM and IW= ',JB,JG,IBM,IW
+!    if (NDBUG < 2) PRINT *,'SWSPCTRL after 1 with JB,JG,IBM and IW= ',JB,JG,IBM,IW
     
     DO JK=1,KLEV
@@ -511 +511 @@
 !      end if
     ENDDO    
-!    if (NDBUG < 2) print *,'SWSPCTRL after 2'
+!    if (NDBUG < 2) PRINT *,'SWSPCTRL after 2'
 
 !-- Delta scaling for clear-sky / aerosol optical quantities
@@ -525 +525 @@
      &   LLRTCHK, ZGCC  , PRMU0, ZTAUC , ZOMCC ,&
      &   ZREFC  , ZREFDC, ZTRAC, ZTRADC )  
-!    if (NDBUG < 2) print *,'SWSPCTR after SWREFTRA for clear-sky'
+!    if (NDBUG < 2) PRINT *,'SWSPCTR after SWREFTRA for clear-sky'
     
 !-- Delta scaling for cloudy quantities
@@ -545 +545 @@
 
     ENDDO
-!    if (NDBUG < 2) print *,'SWSPCTR after Delta scaling'
+!    if (NDBUG < 2) PRINT *,'SWSPCTR after Delta scaling'
     
     CALL SRTM_REFTRA ( KLEV, I_KMODTS ,&
      &   LLRTCHK, ZGCO  , PRMU0, ZTAUO , ZOMCO ,&
      &   ZREFO , ZREFDO, ZTRAO, ZTRADO )  
-!    if (NDBUG < 2) print *,'SWSPCTR after SWREFTRA for cloudy'
+!    if (NDBUG < 2) PRINT *,'SWSPCTR after SWREFTRA for cloudy'
 
     DO JK=1,KLEV
@@ -582 +582 @@
 
     ENDDO           
-!    if (NDBUG < 2) print *,'SRTM_SPCVRT after combining clear and cloudy'
+!    if (NDBUG < 2) PRINT *,'SRTM_SPCVRT after combining clear and cloudy'
                  
 !-- vertical quadrature producing clear-sky fluxes
 
-!    print *,'SRTM_SPCVRT after 3 before SRTM_VRTQDR clear'
+!    PRINT *,'SRTM_SPCVRT after 3 before SRTM_VRTQDR clear'
     
     CALL SRTM_VRTQDR ( KLEV, IW ,&
@@ -594 +594 @@
       
 !    IF (NDBUG < 2) THEN
-!      print *,'SRTM_SPCVRT out of SRTM_VRTQDR for clear IW=',IW  
+!      PRINT *,'SRTM_SPCVRT out of SRTM_VRTQDR for clear IW=',IW
 !      DO JK=1,KLEV+1
 !        print 9201,JK,ZCD(JK,IW),ZCU(JK,IW)
@@ -603 +603 @@
 !-- vertical quadrature producing cloudy fluxes
 
-!    print *,'SRTM_SPCVRT after 4 before SRTM_VRTQDR cloudy'
+!    PRINT *,'SRTM_SPCVRT after 4 before SRTM_VRTQDR cloudy'
     
     CALL SRTM_VRTQDR ( KLEV, IW ,&
@@ -611 +611 @@
  
 !    IF (NDBUG < 2) THEN     
-!      print *,'SRTM_SPCVRT out of SRTM_VRTQDR for cloudy IW=',IW
+!      PRINT *,'SRTM_SPCVRT out of SRTM_VRTQDR for cloudy IW=',IW
 !      DO JK=1,KLEV+1
 !        print 9202,JK,ZFD(JK,IW),ZFU(JK,IW)
@@ -655 +655 @@
     ENDDO
 
-!    if (NDBUG < 2) print *,'SRTM_SPCVRT end of JG=',JG,' for JB=',JB,' i.e. IW=',IW
+!    if (NDBUG < 2) PRINT *,'SRTM_SPCVRT end of JG=',JG,' for JB=',JB,' i.e. IW=',IW
   ENDDO             
 !-- end loop on JG
 
-!  print *,' --- JB= ',JB,' with IB1, IB2= ',IB1,IB2
+!  PRINT *,' --- JB= ',JB,' with IB1, IB2= ',IB1,IB2
 ENDDO                    
 !-- end loop on JB
-!if (NDBUG < 2) print *,'SRTM_SPCVRT about to come out'
+!if (NDBUG < 2) PRINT *,'SRTM_SPCVRT about to come out'
 
 !DO IBM=1,14