Changeset 921 for trunk/NEMO/LIM_SRC_3/limcons.F90
- Timestamp:
- 2008-05-13T10:28:52+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/limcons.F90
r834 r921 42 42 CONTAINS 43 43 44 !===============================================================================44 !=============================================================================== 45 45 46 46 SUBROUTINE lim_column_sum(nsum,xin,xout) 47 ! !!-------------------------------------------------------------------48 ! !! *** ROUTINE lim_column_sum ***49 ! !!50 ! !! ** Purpose : Compute the sum of xin over nsum categories51 ! !!52 ! !! ** Method : Arithmetics53 ! !!54 ! !! ** Action : Gets xin(ji,jj,jl) and computes xout(ji,jj)55 ! !!56 ! !! History :57 ! !! author: William H. Lipscomb, LANL58 ! !! 2.1 ! 04-06 (M. Vancoppenolle) Energy Conservation59 ! !!---------------------------------------------------------------------60 ! !! * Local variables47 ! !!------------------------------------------------------------------- 48 ! !! *** ROUTINE lim_column_sum *** 49 ! !! 50 ! !! ** Purpose : Compute the sum of xin over nsum categories 51 ! !! 52 ! !! ** Method : Arithmetics 53 ! !! 54 ! !! ** Action : Gets xin(ji,jj,jl) and computes xout(ji,jj) 55 ! !! 56 ! !! History : 57 ! !! author: William H. Lipscomb, LANL 58 ! !! 2.1 ! 04-06 (M. Vancoppenolle) Energy Conservation 59 ! !!--------------------------------------------------------------------- 60 ! !! * Local variables 61 61 INTEGER, INTENT(in) :: & 62 62 nsum ! number of categories/layers 63 63 64 64 REAL (wp), DIMENSION(jpi, jpj, jpl), INTENT(IN) :: & 65 65 xin ! input field 66 66 67 67 REAL (wp), DIMENSION(jpi, jpj), INTENT(OUT) :: & 68 68 xout ! output field 69 69 INTEGER :: & 70 71 72 ! !!---------------------------------------------------------------------73 ! WRITE(numout,*) ' lim_column_sum '74 ! WRITE(numout,*) ' ~~~~~~~~~~~~~~ '70 ji, jj, jl ! horizontal indices 71 72 ! !!--------------------------------------------------------------------- 73 ! WRITE(numout,*) ' lim_column_sum ' 74 ! WRITE(numout,*) ' ~~~~~~~~~~~~~~ ' 75 75 76 76 xout(:,:) = 0.00 … … 86 86 END SUBROUTINE lim_column_sum 87 87 88 !===============================================================================88 !=============================================================================== 89 89 90 90 SUBROUTINE lim_column_sum_energy(nsum,nlay,xin,xout) … … 106 106 !! * Local variables 107 107 INTEGER, INTENT(in) :: & 108 109 108 nsum, & !: number of categories 109 nlay !: number of vertical layers 110 110 111 111 REAL (wp), DIMENSION(jpi, jpj, jkmax, jpl), INTENT(IN) :: & 112 112 xin !: input field 113 113 114 114 REAL (wp), DIMENSION(jpi, jpj), INTENT(OUT) :: & 115 115 xout !: output field 116 116 117 117 INTEGER :: & 118 119 120 !!--------------------------------------------------------------------- 121 122 ! WRITE(numout,*) ' lim_column_sum_energy '123 ! WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~ '118 ji, jj, & !: horizontal indices 119 jk, jl !: layer and category indices 120 !!--------------------------------------------------------------------- 121 122 ! WRITE(numout,*) ' lim_column_sum_energy ' 123 ! WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~ ' 124 124 125 125 xout(:,:) = 0.00 … … 137 137 END SUBROUTINE lim_column_sum_energy 138 138 139 !===============================================================================140 139 !=============================================================================== 140 141 141 SUBROUTINE lim_cons_check(x1, x2, max_err, fieldid) 142 142 !!------------------------------------------------------------------- … … 206 206 WRITE (numout,*) ' Point : ', ji, jj 207 207 WRITE (numout,*) ' lat, lon : ', gphit(ji,jj), & 208 208 glamt(ji,jj) 209 209 WRITE (numout,*) ' Initial value : ', x1(ji,jj) 210 210 WRITE (numout,*) ' Final value : ', x2(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.