Ignore:
Timestamp:
Apr 30, 2015 7:05:52 AM (6 years ago)
Author:
maronga
Message:

Added support for RRTMG radiation code

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/header.f90

    r1576 r1585  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Further output for radiation model(s).
    2323!
    2424! Former revisions:
     
    245245
    246246    USE radiation_model_mod,                                                   &
    247         ONLY:  albedo, day_init, dt_radiation, lambda, net_radiation,          &
    248                radiation, radiation_scheme, time_utc_init
     247        ONLY:  albedo, albedo_type, albedo_type_name, constant_albedo,         &
     248               day_init, dt_radiation, lambda, lw_radiation, net_radiation,    &
     249               radiation, radiation_scheme, sw_radiation, time_utc_init
    249250   
    250251    USE spectrum,                                                              &
     
    927928    IF ( radiation )  THEN
    928929!
    929 !--    Write land surface model header
     930!--    Write radiation model header
    930931       WRITE( io, 444 )
    931932
     
    934935       ELSEIF ( radiation_scheme == "clear-sky" )  THEN
    935936          WRITE( io, 446 )
    936        ELSE
    937           WRITE( io, 447 ) radiation_scheme
    938        ENDIF
    939 
    940        WRITE( io, 448 ) albedo
     937       ELSEIF ( radiation_scheme == "rrtmg" )  THEN
     938          WRITE( io, 447 )
     939          IF ( .NOT. lw_radiation )  WRITE( io, 458 )
     940          IF ( .NOT. sw_radiation )  WRITE( io, 459 )
     941       ENDIF
     942
     943       IF ( albedo_type == 0 )  THEN
     944          WRITE( io, 448 ) albedo
     945       ELSE
     946          WRITE( io, 456 ) albedo_type_name(albedo_type)
     947       ENDIF
     948       IF ( constant_albedo )  THEN
     949          WRITE( io, 457 )
     950       ENDIF
    941951       WRITE( io, 449 ) dt_radiation
    942 
    943952    ENDIF
    944953
     
    22622271446 FORMAT (' --> Simple radiation scheme for clear sky is used (no clouds,',  &
    22632272                   ' default)')
    2264 447 FORMAT (' --> Radiation scheme:', A)
    2265 448 FORMAT (/'    Surface albedo: albedo = ', F5.3)
    2266 449 FORMAT  ('    Timestep: dt_radiation = ', F5.2, '  s')
     2273447 FORMAT (' --> RRTMG scheme is used')
     2274448 FORMAT (/'     User-specific surface albedo: albedo = ', F5.3)
     2275449 FORMAT  ('     Timestep: dt_radiation = ', F5.2, '  s')
    22672276
    22682277450 FORMAT (//' LES / Turbulence quantities:'/ &
     
    22732282454 FORMAT ('    TKE is not allowed to fall below ',E9.2,' (m/s)**2')
    22742283455 FORMAT ('    initial TKE is prescribed as ',E9.2,' (m/s)**2')
     2284456 FORMAT (/'    Albedo is set for land surface type: ', A)
     2285457 FORMAT (/'    --> Albedo is fixed during the run')
     2286458 FORMAT (/'    --> Longwave radiation is disabled')
     2287459 FORMAT (/'    --> Shortwave radiation is disabled.')
    22752288470 FORMAT (//' Actions during the simulation:'/ &
    22762289              ' -----------------------------'/)
Note: See TracChangeset for help on using the changeset viewer.