Changeset 1853 for palm/trunk
- Timestamp:
- Apr 11, 2016 9:00:35 AM (9 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/flow_statistics.f90
r1852 r1853 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 21 ! Adjusted for use with radiation_scheme = constant 22 22 ! 23 23 ! Former revisions: … … 798 798 ENDIF 799 799 800 IF ( radiation ) THEN800 IF ( radiation .AND. radiation_scheme /= 'constant' ) THEN 801 801 sums_l(nzb,101,tn) = sums_l(nzb,101,tn) + rad_net(j,i) 802 802 sums_l(nzb,102,tn) = sums_l(nzb,102,tn) + rad_lw_in(nzb,j,i) -
palm/trunk/SOURCE/radiation_model_mod.f90
r1852 r1853 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 21 ! Added routine for radiation_scheme = constant. 22 22 ! 23 23 ! Former revisions: … … 374 374 MODULE PROCEDURE radiation_clearsky 375 375 END INTERFACE radiation_clearsky 376 376 377 INTERFACE radiation_constant 378 MODULE PROCEDURE radiation_constant 379 END INTERFACE radiation_constant 380 377 381 INTERFACE radiation_header 378 382 MODULE PROCEDURE radiation_header … … 403 407 !-- Public functions 404 408 PUBLIC radiation_check_data_output, radiation_check_data_output_pr, & 405 radiation_check_parameters, radiation_clearsky, radiation_header, & 406 radiation_init, radiation_parin, radiation_rrtmg, radiation_tendency 409 radiation_check_parameters, radiation_clearsky, radiation_constant, & 410 radiation_header, radiation_init, radiation_parin, radiation_rrtmg, & 411 radiation_tendency 407 412 408 413 ! … … 757 762 IF ( radiation_scheme == 'constant' ) THEN 758 763 rad_net = net_radiation 759 radiation = .FALSE.764 ! radiation = .FALSE. 760 765 ! 761 766 !-- Calculate orbital constants … … 769 774 770 775 776 IF ( radiation_scheme == 'constant' ) THEN 777 778 IF ( .NOT. ALLOCATED ( rad_lw_out ) ) THEN 779 ALLOCATE ( rad_lw_out(0:0,nysg:nyng,nxlg:nxrg) ) 780 ENDIF 781 782 ENDIF 783 771 784 IF ( radiation_scheme == 'clear-sky' ) THEN 772 785 … … 789 802 IF ( .NOT. ALLOCATED ( rad_lw_in ) ) THEN 790 803 ALLOCATE ( rad_lw_in(0:0,nysg:nyng,nxlg:nxrg) ) 791 ENDIF792 IF ( .NOT. ALLOCATED ( rad_lw_out ) ) THEN793 ALLOCATE ( rad_lw_out(0:0,nysg:nyng,nxlg:nxrg) )794 804 ENDIF 795 805 … … 989 999 !-- Calculate radiative fluxes at model start 990 1000 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 991 IF ( radiation_scheme == 'clear-sky' ) THEN 992 CALL radiation_clearsky 993 ELSEIF ( radiation_scheme == 'rrtmg' ) THEN 994 CALL radiation_rrtmg 995 ENDIF 1001 1002 SELECT CASE ( radiation_scheme ) 1003 CASE ( 'rrtmg' ) 1004 CALL radiation_rrtmg 1005 CASE ( 'clear-sky' ) 1006 CALL radiation_clearsky 1007 CASE ( 'constant' ) 1008 CALL radiation_constant 1009 CASE DEFAULT 1010 END SELECT 1011 996 1012 ENDIF 997 1013 … … 1008 1024 SUBROUTINE radiation_clearsky 1009 1025 1010 USE indices, &1011 ONLY: nbgp1012 1026 1013 1027 IMPLICIT NONE … … 1057 1071 END SUBROUTINE radiation_clearsky 1058 1072 1073 1074 !------------------------------------------------------------------------------! 1075 ! Description: 1076 ! ------------ 1077 !> This scheme keeps the prescribed net radiation constant during the run 1078 !------------------------------------------------------------------------------! 1079 SUBROUTINE radiation_constant 1080 1081 1082 IMPLICIT NONE 1083 1084 INTEGER(iwp) :: i, j, k !< loop indices 1085 REAL(wp) :: exn, & !< Exner functions at surface 1086 pt1 !< potential temperature at first grid level 1087 1088 ! 1089 !-- Calculate value of the Exner function 1090 exn = (surface_pressure / 1000.0_wp )**0.286_wp 1091 ! 1092 !-- Prescribe net radiation and estimate a longwave outgoing radiative 1093 !-- flux (needed in land surface model) 1094 DO i = nxlg, nxrg 1095 DO j = nysg, nyng 1096 k = nzb_s_inner(j,i) 1097 1098 rad_net(j,i) = net_radiation 1099 rad_lw_out(0,j,i) = emissivity * sigma_sb * (pt(k,j,i) * exn)**4 1100 1101 ENDDO 1102 ENDDO 1103 1104 END SUBROUTINE radiation_constant 1059 1105 1060 1106 !------------------------------------------------------------------------------! -
palm/trunk/SOURCE/time_integration.f90
r1852 r1853 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! Adjusted for use with radiation_scheme = constant 22 22 ! 23 23 ! Former revisions: … … 295 295 USE radiation_model_mod, & 296 296 ONLY: dt_radiation, force_radiation_call, radiation, & 297 radiation_clearsky, radiation_ rrtmg, radiation_scheme,&298 skip_time_do_radiation, time_radiation297 radiation_clearsky, radiation_constant, radiation_rrtmg, & 298 radiation_scheme, skip_time_do_radiation, time_radiation 299 299 300 300 USE spectra_mod, & … … 863 863 ELSEIF ( radiation_scheme == 'rrtmg' ) THEN 864 864 CALL radiation_rrtmg 865 ELSE 866 CALL radiation_constant 865 867 ENDIF 866 868
Note: See TracChangeset
for help on using the changeset viewer.