Ignore:
Timestamp:
Jul 27, 2016 1:28:04 PM (5 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/average_3d_data.f90

    r1973 r1976  
    101101
    102102    USE radiation_model_mod,                                                   &
    103         ONLY:  rad_net, rad_net_av, rad_lw_in, rad_lw_in_av, rad_lw_out,       &
    104                rad_lw_out_av, rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr,        &
    105                rad_lw_hr_av, rad_sw_in, rad_sw_in_av, rad_sw_out,              &
    106                rad_sw_out_av, rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr,        &
    107                rad_sw_hr_av
     103        ONLY:  radiation, radiation_3d_data_averaging
    108104
    109105
     
    295291             ENDDO
    296292
    297          CASE ( 'rad_net*' )
    298              DO  i = nxlg, nxrg
    299                 DO  j = nysg, nyng
    300                    rad_net_av(j,i) = rad_net_av(j,i) / REAL( average_count_3d, KIND=wp )
    301                 ENDDO
    302              ENDDO
    303 
    304           CASE ( 'rad_lw_in' )
    305              DO  i = nxlg, nxrg
    306                 DO  j = nysg, nyng
    307                    DO  k = nzb, nzt+1
    308                       rad_lw_in_av(k,j,i) = rad_lw_in_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    309                    ENDDO
    310                 ENDDO
    311              ENDDO
    312 
    313           CASE ( 'rad_lw_out' )
    314              DO  i = nxlg, nxrg
    315                 DO  j = nysg, nyng
    316                    DO  k = nzb, nzt+1
    317                       rad_lw_out_av(k,j,i) = rad_lw_out_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    318                    ENDDO
    319                 ENDDO
    320              ENDDO
    321 
    322           CASE ( 'rad_lw_cs_hr' )
    323              DO  i = nxlg, nxrg
    324                 DO  j = nysg, nyng
    325                    DO  k = nzb, nzt+1
    326                       rad_lw_cs_hr_av(k,j,i) = rad_lw_cs_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    327                    ENDDO
    328                 ENDDO
    329              ENDDO
    330 
    331           CASE ( 'rad_lw_hr' )
    332              DO  i = nxlg, nxrg
    333                 DO  j = nysg, nyng
    334                    DO  k = nzb, nzt+1
    335                       rad_lw_hr_av(k,j,i) = rad_lw_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    336                    ENDDO
    337                 ENDDO
    338              ENDDO
    339 
    340           CASE ( 'rad_sw_in' )
    341              DO  i = nxlg, nxrg
    342                 DO  j = nysg, nyng
    343                    DO  k = nzb, nzt+1
    344                       rad_sw_in_av(k,j,i) = rad_sw_in_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    345                    ENDDO
    346                 ENDDO
    347              ENDDO
    348 
    349           CASE ( 'rad_sw_out' )
    350              DO  i = nxlg, nxrg
    351                 DO  j = nysg, nyng
    352                    DO  k = nzb, nzt+1
    353                       rad_sw_out_av(k,j,i) = rad_sw_out_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    354                    ENDDO
    355                 ENDDO
    356              ENDDO
    357 
    358           CASE ( 'rad_sw_cs_hr' )
    359              DO  i = nxlg, nxrg
    360                 DO  j = nysg, nyng
    361                    DO  k = nzb, nzt+1
    362                       rad_sw_cs_hr_av(k,j,i) = rad_sw_cs_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    363                    ENDDO
    364                 ENDDO
    365              ENDDO
    366 
    367           CASE ( 'rad_sw_hr' )
    368              DO  i = nxlg, nxrg
    369                 DO  j = nysg, nyng
    370                    DO  k = nzb, nzt+1
    371                       rad_sw_hr_av(k,j,i) = rad_sw_hr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    372                    ENDDO
    373                 ENDDO
    374              ENDDO
    375 
    376293          CASE ( 'rho' )
    377294             DO  i = nxlg, nxrg
     
    487404
    488405!
     406!--          Radiation quantity
     407             IF ( radiation )  THEN
     408                CALL radiation_3d_data_averaging( 'average', doav(ii) )
     409             ENDIF
     410
     411!
    489412!--          User-defined quantity
    490413             CALL user_3d_data_averaging( 'average', doav(ii) )
Note: See TracChangeset for help on using the changeset viewer.