Ignore:
Timestamp:
Jul 27, 2016 1:28:04 PM (8 years ago)
Author:
maronga
Message:

further modularization of land surface model (2D/3D output and restart data). Bugfix for restart runs without land surface model

File:
1 edited

Legend:

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

    r1973 r1976  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! Output of radiation quantities is now done directly in the respective module
    2222!
    2323! Former revisions:
     
    165165
    166166    USE radiation_model_mod,                                                   &
    167         ONLY:  rad_lw_in, rad_lw_in_av, rad_lw_out, rad_lw_out_av,             &
    168                rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr, rad_lw_hr_av,         &
    169                rad_sw_in, rad_sw_in_av, rad_sw_out, rad_sw_out_av,             &
    170                rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr, rad_sw_hr_av
     167        ONLY:  radiation, radiation_data_output_3d
    171168
    172169
     
    494491             ENDIF
    495492
    496           CASE ( 'rad_sw_in' )
    497              IF ( av == 0 )  THEN
    498                 to_be_resorted => rad_sw_in
    499              ELSE
    500                 to_be_resorted => rad_sw_in_av
    501              ENDIF
    502 
    503           CASE ( 'rad_sw_out' )
    504              IF ( av == 0 )  THEN
    505                 to_be_resorted => rad_sw_out
    506              ELSE
    507                 to_be_resorted => rad_sw_out_av
    508              ENDIF
    509 
    510           CASE ( 'rad_sw_cs_hr' )
    511              IF ( av == 0 )  THEN
    512                 to_be_resorted => rad_sw_cs_hr
    513              ELSE
    514                 to_be_resorted => rad_sw_cs_hr_av
    515              ENDIF
    516 
    517           CASE ( 'rad_sw_hr' )
    518              IF ( av == 0 )  THEN
    519                 to_be_resorted => rad_sw_hr
    520              ELSE
    521                 to_be_resorted => rad_sw_hr_av
    522              ENDIF
    523 
    524           CASE ( 'rad_lw_in' )
    525              IF ( av == 0 )  THEN
    526                 to_be_resorted => rad_lw_in
    527              ELSE
    528                 to_be_resorted => rad_lw_in_av
    529              ENDIF
    530 
    531           CASE ( 'rad_lw_out' )
    532              IF ( av == 0 )  THEN
    533                 to_be_resorted => rad_lw_out
    534              ELSE
    535                 to_be_resorted => rad_lw_out_av
    536              ENDIF
    537 
    538           CASE ( 'rad_lw_cs_hr' )
    539              IF ( av == 0 )  THEN
    540                 to_be_resorted => rad_lw_cs_hr
    541              ELSE
    542                 to_be_resorted => rad_lw_cs_hr_av
    543              ENDIF
    544 
    545           CASE ( 'rad_lw_hr' )
    546              IF ( av == 0 )  THEN
    547                 to_be_resorted => rad_lw_hr
    548              ELSE
    549                 to_be_resorted => rad_lw_hr_av
    550              ENDIF
    551 
    552493          CASE ( 'rho' )
    553494             IF ( av == 0 )  THEN
     
    613554
    614555                CALL lsm_data_output_3d( av, do3d(av,if), found, local_pf )
     556                resorted = .TRUE.
     557
     558!
     559!--             If no soil model variable was found, re-allocate local_pf
     560                IF ( .NOT. found )  THEN
     561                   nzb_do = nzb
     562                   nzt_do = nz_do3d
     563
     564                   DEALLOCATE ( local_pf )
     565                   ALLOCATE( local_pf(nxlg:nxrg,nysg:nyng,nzb_do:nzt_do) )                 
     566                ENDIF
     567
     568             ENDIF
     569
     570!
     571!--          Radiation quantity
     572             IF ( .NOT. found  .AND.  radiation )  THEN
     573                CALL radiation_data_output_3d( av, do3d(av,if), found,         &
     574                                               local_pf )
    615575                resorted = .TRUE.
    616576             ENDIF
Note: See TracChangeset for help on using the changeset viewer.