Changeset 1585 for palm/trunk/SOURCE/header.f90

Timestamp:

Apr 30, 2015 7:05:52 AM (9 years ago)

Author:

maronga

Message:

Added support for RRTMG radiation code

File:

• palm/trunk/SOURCE/header.f90 (modified) (6 diffs)

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Further output for radiation model(s).
 !
 ! Former revisions:
@@ -245 +245 @@
 
     USE radiation_model_mod,                                                   &
-        ONLY:  albedo, day_init, dt_radiation, lambda, net_radiation,          &
-               radiation, radiation_scheme, time_utc_init
+        ONLY:  albedo, albedo_type, albedo_type_name, constant_albedo,         &
+               day_init, dt_radiation, lambda, lw_radiation, net_radiation,    &
+               radiation, radiation_scheme, sw_radiation, time_utc_init
     
     USE spectrum,                                                              &
@@ -927 +928 @@
     IF ( radiation )  THEN
 !
-!--    Write land surface model header
+!--    Write radiation model header
        WRITE( io, 444 )
 
@@ -934 +935 @@
        ELSEIF ( radiation_scheme == "clear-sky" )  THEN
           WRITE( io, 446 )
-       ELSE
-          WRITE( io, 447 ) radiation_scheme
-       ENDIF
-
-       WRITE( io, 448 ) albedo
+       ELSEIF ( radiation_scheme == "rrtmg" )  THEN
+          WRITE( io, 447 )
+          IF ( .NOT. lw_radiation )  WRITE( io, 458 )
+          IF ( .NOT. sw_radiation )  WRITE( io, 459 )
+       ENDIF
+
+       IF ( albedo_type == 0 )  THEN
+          WRITE( io, 448 ) albedo
+       ELSE
+          WRITE( io, 456 ) albedo_type_name(albedo_type)
+       ENDIF
+       IF ( constant_albedo )  THEN
+          WRITE( io, 457 )
+       ENDIF
        WRITE( io, 449 ) dt_radiation
-
     ENDIF
 
@@ -2262 +2271 @@
 446 FORMAT (' --> Simple radiation scheme for clear sky is used (no clouds,',  &
                    ' default)')
-447 FORMAT (' --> Radiation scheme:', A)
-448 FORMAT (/'    Surface albedo: albedo = ', F5.3)
-449 FORMAT  ('    Timestep: dt_radiation = ', F5.2, '  s')
+447 FORMAT (' --> RRTMG scheme is used')
+448 FORMAT (/'     User-specific surface albedo: albedo = ', F5.3)
+449 FORMAT  ('     Timestep: dt_radiation = ', F5.2, '  s')
 
 450 FORMAT (//' LES / Turbulence quantities:'/ &
@@ -2273 +2282 @@
 454 FORMAT ('    TKE is not allowed to fall below ',E9.2,' (m/s)**2')
 455 FORMAT ('    initial TKE is prescribed as ',E9.2,' (m/s)**2')
+456 FORMAT (/'    Albedo is set for land surface type: ', A)
+457 FORMAT (/'    --> Albedo is fixed during the run')
+458 FORMAT (/'    --> Longwave radiation is disabled')
+459 FORMAT (/'    --> Shortwave radiation is disabled.')
 470 FORMAT (//' Actions during the simulation:'/ &
               ' -----------------------------'/)