Changeset 6628 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z
- Timestamp:
- 2016-05-26T11:54:20+02:00 (8 years ago)
- Location:
- trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90
r6140 r6628 74 74 REAL(wp) :: zchl 75 75 REAL(wp) :: zc0 , zc1 , zc2, zc3, z1_dep 76 REAL(wp), POINTER, DIMENSION(:,: ) :: zdepmoy, zetmp1, zetmp2, zetmp3, zetmp4, zqsr100 76 REAL(wp), POINTER, DIMENSION(:,: ) :: zdepmoy, zetmp1, zetmp2, zetmp3, zetmp4 77 REAL(wp), POINTER, DIMENSION(:,: ) :: zqsr100, zqsr_corr 77 78 REAL(wp), POINTER, DIMENSION(:,:,:) :: zpar, ze0, ze1, ze2, ze3 78 79 !!--------------------------------------------------------------------- … … 81 82 ! 82 83 ! Allocate temporary workspace 83 CALL wrk_alloc( jpi, jpj, zqsr100, zdepmoy, zetmp1, zetmp2, zetmp3, zetmp4 ) 84 CALL wrk_alloc( jpi, jpj, jpk, zpar, ze0, ze1, ze2, ze3 ) 84 CALL wrk_alloc( jpi, jpj, zdepmoy, zetmp1, zetmp2, zetmp3, zetmp4 ) 85 CALL wrk_alloc( jpi, jpj, zqsr100, zqsr_corr ) 86 CALL wrk_alloc( jpi, jpj, jpk, zpar , ze0, ze1, ze2, ze3 ) 85 87 86 88 IF( knt == 1 .AND. ln_varpar ) CALL p4z_opt_sbc( kt ) … … 108 110 ! ! -------------------------------------- 109 111 IF( l_trcdm2dc ) THEN ! diurnal cycle 110 ! 1% of qsr to compute euphotic layer 111 zqsr100(:,:) = 0.01 * qsr_mean(:,:) ! daily mean qsr 112 ! 113 CALL p4z_opt_par( kt, qsr_mean, ze1, ze2, ze3 ) 112 ! ! 1% of qsr to compute euphotic layer 113 zqsr100(:,:) = 0.01 * qsr_mean(:,:) ! daily mean qsr 114 ! 115 zqsr_corr(:,:) = qsr_mean(:,:) / ( 1. - fr_i(:,:) + rtrn ) 116 ! 117 CALL p4z_opt_par( kt, zqsr_corr, ze1, ze2, ze3 ) 114 118 ! 115 119 DO jk = 1, nksrp … … 119 123 END DO 120 124 ! 121 CALL p4z_opt_par( kt, qsr, ze1, ze2, ze3 ) 125 zqsr_corr(:,:) = qsr(:,:) / ( 1. - fr_i(:,:) + rtrn ) 126 ! 127 CALL p4z_opt_par( kt, zqsr_corr, ze1, ze2, ze3 ) 122 128 ! 123 129 DO jk = 1, nksrp … … 126 132 ! 127 133 ELSE 128 ! 1% of qsr to compute euphotic layer 129 zqsr100(:,:) = 0.01 * qsr(:,:) 130 ! 131 CALL p4z_opt_par( kt, qsr, ze1, ze2, ze3 ) 134 ! ! 1% of qsr to compute euphotic layer 135 zqsr100(:,:) = 0.01 * qsr(:,:) ! daily mean qsr 136 ! 137 zqsr_corr(:,:) = qsr(:,:) / ( 1. - fr_i(:,:) + rtrn ) 138 ! 139 CALL p4z_opt_par( kt, zqsr_corr, ze1, ze2, ze3 ) 132 140 ! 133 141 DO jk = 1, nksrp … … 218 226 ENDIF 219 227 ! 220 CALL wrk_dealloc( jpi, jpj, zqsr100, zdepmoy, zetmp1, zetmp2, zetmp3, zetmp4 ) 221 CALL wrk_dealloc( jpi, jpj, jpk, zpar, ze0, ze1, ze2, ze3 ) 228 CALL wrk_dealloc( jpi, jpj, zdepmoy, zetmp1, zetmp2, zetmp3, zetmp4 ) 229 CALL wrk_dealloc( jpi, jpj, zqsr100, zqsr_corr ) 230 CALL wrk_dealloc( jpi, jpj, jpk, zpar , ze0, ze1, ze2, ze3 ) 222 231 ! 223 232 IF( nn_timing == 1 ) CALL timing_stop('p4z_opt') -
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90
r6140 r6628 133 133 zval = MAX( 1., zstrn(ji,jj) ) 134 134 zval = 1.5 * zval / ( 12. + zval ) 135 zprbio(ji,jj,jk) = prmax(ji,jj,jk) * zval 135 zprbio(ji,jj,jk) = prmax(ji,jj,jk) * zval * ( 1. - fr_i(ji,jj) ) 136 136 zprdia(ji,jj,jk) = zprbio(ji,jj,jk) 137 137 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.