Changeset 9330
- Timestamp:
- 2018-02-16T18:29:04+01:00 (6 years ago)
- Location:
- branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90
r8173 r9330 392 392 zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 393 393 CALL iom_put( "PhytoC25h", zw3d ) ! total phytoplankton 394 zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:))395 CALL iom_put( "PhyCMol25h", zw3d ) ! total phytoplankton in mmol396 zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))397 CALL iom_put( "ZooC25h", zw3d ) ! total zooplankton398 zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:))399 CALL iom_put( "ZooCMol25h", zw3d ) ! total zooplankton in mmol394 !zw3d(:,:,:) = fabm3d_25h(:,:,:,11)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:)) 395 !CALL iom_put( "PhyCMol25h", zw3d ) ! total phytoplankton in mmol 396 !zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 397 !CALL iom_put( "ZooC25h", zw3d ) ! total zooplankton 398 !zw3d(:,:,:) = fabm3d_25h(:,:,:,12)*tmask(:,:,:)*0.0832590939745 + zmdi*(1.0-tmask(:,:,:)) 399 !CALL iom_put( "ZooCMol25h", zw3d ) ! total zooplankton in mmol 400 400 #endif 401 401 -
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/par_fabm.F90
r8049 r9330 28 28 jp_fabm_netp2,jp_fabm_netp3, & 29 29 jp_fabm_netp4,jp_fabm_o3ph, & 30 jp_fabm_o3pc, jp_fabm_xeps 30 jp_fabm_o3pc, jp_fabm_xeps, & 31 jp_fabm_o3c, jp_fabm_o3a, & 32 jp_fabm_r4n, jp_fabm_r6n, & 33 jp_fabm_r8n 31 34 32 35 INTEGER, PUBLIC :: jp_fabmdia_3dout = 12 -
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90
r8050 r9330 104 104 jp_fabm_n5s = fabm_state_index( 'N5_s' ) 105 105 jp_fabm_o2o = fabm_state_index( 'O2_o' ) 106 jp_fabm_o3c = fabm_state_index( 'O3_c' ) 107 jp_fabm_o3a = fabm_state_index( 'O3_bioalk' ) 108 jp_fabm_r4n = fabm_state_index( 'R4_n' ) 109 jp_fabm_r6n = fabm_state_index( 'R6_n' ) 110 jp_fabm_r8n = fabm_state_index( 'R8_n' ) 106 111 !jp_fabm_netp1= fabm_diag_index( 'P1_netP' ) 107 112 !jp_fabm_netp2= fabm_diag_index( 'P2_netP' ) … … 181 186 END IF 182 187 IF( lk_mpp ) CALL mppsync !Ensure field_def_fabm is ready. 188 189 MLD_MAX(:,:) = 0.0 190 PGROW_AVG(:,:) = 0.0 191 PLOSS_AVG(:,:) = 0.0 192 PHYT_AVG(:,:) = 0.0 183 193 184 194 1000 FORMAT (A) -
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90
r7835 r9330 33 33 USE inputs_fabm 34 34 USE vertical_movement_fabm 35 USE zdfmxl 35 36 36 37 !USE fldread ! time interpolation … … 113 114 114 115 CALL st2d_fabm_nxt( kt ) 116 117 CALL asmdiags_fabm 115 118 116 119 IF( l_trdtrc ) CALL wrk_alloc( jpi, jpj, jpk, ztrfabm ) … … 130 133 131 134 END SUBROUTINE trc_sms_fabm 135 136 SUBROUTINE asmdiags_fabm 137 INTEGER :: ji,jj,jk,jkmax 138 139 PGROW_AVG(:,:) = 1.0 140 PLOSS_AVG(:,:) = 1.0 141 PHYT_AVG(:,:) = 0.0 142 143 DO jj = 1, jpj 144 DO ji = 1, jpi 145 ! 146 jkmax = jpk-1 147 DO jk = jpk-1, 1, -1 148 IF ( ( hmld_tref(ji,jj) > gdepw_n(ji,jj,jk) ) .AND. & 149 & ( hmld_tref(ji,jj) <= gdepw_n(ji,jj,jk+1) ) ) THEN 150 hmld_tref(ji,jj) = gdepw_n(ji,jj,jk+1) 151 jkmax = jk 152 ENDIF 153 END DO 154 ! 155 DO jk = 1, jkmax 156 PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) + & 157 & trn(ji,jj,jk,jp_fabm_p1n) + & 158 & trn(ji,jj,jk,jp_fabm_p2n) + & 159 & trn(ji,jj,jk,jp_fabm_p3n) + & 160 & trn(ji,jj,jk,jp_fabm_p4n) 161 END DO 162 163 PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) / REAL(jkmax) 164 165 IF ( hmld_tref(ji,jj) .GT. MLD_MAX(ji,jj) ) THEN 166 MLD_MAX(ji,jj) = hmld_tref(ji,jj) 167 ENDIF 168 ! 169 END DO 170 END DO 171 172 END SUBROUTINE asmdiags_fabm 132 173 133 174 SUBROUTINE compute_fabm() -
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/trc.F90
r8049 r9330 225 225 #endif 226 226 227 #if defined key_fabm 228 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: PGROW_AVG ! Phytoplankton growth for use in ASM code 229 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: PLOSS_AVG ! Phytoplankton loss for use in ASM code 230 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: PHYT_AVG ! Phytoplankton for use in ASM code 231 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: MLD_MAX ! Maximum MLD for use in ASM code 232 #endif 233 227 234 !!---------------------------------------------------------------------- 228 235 !! NEMO/TOP 3.3.1 , NEMO Consortium (2010) … … 253 260 ! FABM <<<+++ 254 261 & ln_trc_sbc(jptra) , ln_trc_cbc(jptra) , ln_trc_obc(jptra) , & 262 & PGROW_AVG(jpi,jpj) , PLOSS_AVG(jpi,jpj) , PHYT_AVG(jpi,jpj) , & 263 & MLD_MAX(jpi,jpj) , & 255 264 #endif 256 265 #if defined key_bdy
Note: See TracChangeset
for help on using the changeset viewer.