Changeset 12276 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zligand.F90
- Timestamp:
- 2019-12-20T12:14:26+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/TOP/PISCES/P4Z/p4zligand.F90
r11536 r12276 43 43 INTEGER :: ji, jj, jk 44 44 REAL(wp) :: zlgwp, zlgwpr, zlgwr, zlablgw 45 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zligrem, zligpr, zrligprod 46 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3d 45 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zligrem, zligpr, zligprod 47 46 CHARACTER (len=25) :: charout 48 47 !!--------------------------------------------------------------------- … … 69 68 zligrem(ji,jj,jk) = zlgwr 70 69 zligpr(ji,jj,jk) = zlgwpr 71 z rligprod(ji,jj,jk)= zlgwp70 zligprod(ji,jj,jk) = zlgwp 72 71 ! 73 72 END DO … … 78 77 ! --------------------------------- 79 78 IF( lk_iomput .AND. knt == nrdttrc ) THEN 80 ALLOCATE( zw3d(jpi,jpj,jpk) )81 79 IF( iom_use( "LIGREM" ) ) THEN 82 zw3d(:,:,:) = zligrem(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) 83 CALL iom_put( "LIGREM", zw3d ) 80 zligrem(:,:,jpk) = 0. ; CALL iom_put( "LIGREM", zligrem(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 84 81 ENDIF 85 82 IF( iom_use( "LIGPR" ) ) THEN 86 zw3d(:,:,:) = zligpr(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) 87 CALL iom_put( "LIGPR", zw3d ) 83 zligpr(:,:,jpk) = 0. ; CALL iom_put( "LIGPR" , zligpr(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 88 84 ENDIF 89 85 IF( iom_use( "LPRODR" ) ) THEN 90 zw3d(:,:,:) = zrligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) 91 CALL iom_put( "LPRODR", zw3d ) 86 zligprod(:,:,jpk) = 0. ; CALL iom_put( "LPRODR", zligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 92 87 ENDIF 93 DEALLOCATE( zw3d )94 88 ENDIF 95 89 ! … … 125 119 WRITE(numout,*) '~~~~~~~~~~~~~~~' 126 120 ENDIF 127 REWIND( numnatp_ref ) ! Namelist nampislig in reference namelist : Pisces remineralization 121 122 REWIND( numnatp_ref ) 128 123 READ ( numnatp_ref, nampislig, IOSTAT = ios, ERR = 901) 129 124 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislig in reference namelist' ) 130 REWIND( numnatp_cfg ) ! Namelist nampislig in configuration namelist : Pisces remineralization 125 126 REWIND( numnatp_cfg ) 131 127 READ ( numnatp_cfg, nampislig, IOSTAT = ios, ERR = 902 ) 132 128 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampislig in configuration namelist' )
Note: See TracChangeset
for help on using the changeset viewer.