Changeset 12377 for NEMO/trunk/src/TOP/C14/trcwri_c14.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/TOP/C14/trcwri_c14.F90
r10425 r12377 23 23 ! Standard ratio: 1.176E-12 ; Avogadro's nbr = 6.022E+23 at/mol ; bomb C14 traditionally reported as 1.E+26 atoms 24 24 REAL(wp), PARAMETER :: atomc14 = 1.176 * 6.022E-15 ! conversion factor 25 !! * Substitutions 26 # include "do_loop_substitute.h90" 25 27 26 28 27 29 CONTAINS 28 30 29 SUBROUTINE trc_wri_c14 31 SUBROUTINE trc_wri_c14( Kmm ) 30 32 !!--------------------------------------------------------------------- 31 33 !! *** ROUTINE trc_wri_c14 *** … … 33 35 !! ** Purpose : output additional C14 tracers fields 34 36 !!--------------------------------------------------------------------- 37 INTEGER, INTENT(in) :: Kmm ! time level indices 35 38 CHARACTER (len=20) :: cltra ! short title for tracer 36 39 INTEGER :: ji,jj,jk,jn ! dummy loop indexes … … 43 46 ! --------------------------------------- 44 47 cltra = TRIM( ctrcnm(jp_c14) ) ! short title for tracer 45 CALL iom_put( cltra, tr n(:,:,:,jp_c14) )48 CALL iom_put( cltra, tr(:,:,:,jp_c14,Kmm) ) 46 49 47 50 ! compute and write the tracer diagnostic in the file … … 57 60 zz3d(:,:,:) = 0._wp 58 61 ! 59 DO jk = 1, jpkm1 60 DO jj = 1, jpj 61 DO ji = 1, jpi 62 IF( tmask(ji,jj,jk) > 0._wp) THEN 63 z3d (ji,jj,jk) = trn(ji,jj,jk,jp_c14) 64 zz3d(ji,jj,jk) = LOG( z3d(ji,jj,jk) ) 65 ENDIF 66 ENDDO 67 ENDDO 68 ENDDO 62 DO_3D_11_11( 1, jpkm1 ) 63 IF( tmask(ji,jj,jk) > 0._wp) THEN 64 z3d (ji,jj,jk) = tr(ji,jj,jk,jp_c14,Kmm) 65 zz3d(ji,jj,jk) = LOG( z3d(ji,jj,jk) ) 66 ENDIF 67 END_3D 69 68 zres(:,:) = z3d(:,:,1) 70 69 … … 72 71 z2d(:,:) =0._wp 73 72 jk = 1 74 DO jj = 1, jpj 75 DO ji = 1, jpi 76 ztemp = zres(ji,jj) / c14sbc(ji,jj) 77 IF( ztemp > 0._wp .AND. tmask(ji,jj,jk) > 0._wp ) z2d(ji,jj) = LOG( ztemp ) 78 ENDDO 79 ENDDO 73 DO_2D_11_11 74 ztemp = zres(ji,jj) / c14sbc(ji,jj) 75 IF( ztemp > 0._wp .AND. tmask(ji,jj,jk) > 0._wp ) z2d(ji,jj) = LOG( ztemp ) 76 END_2D 80 77 ! 81 78 z3d(:,:,:) = 1.d03 * ( z3d(:,:,:) - 1._wp ) … … 113 110 ENDIF 114 111 IF( iom_use("C14Inv") ) THEN 115 ztemp = glob_sum( 'trcwri_c14', tr n(:,:,:,jp_c14) * cvol(:,:,:) )112 ztemp = glob_sum( 'trcwri_c14', tr(:,:,:,jp_c14,Kmm) * cvol(:,:,:) ) 116 113 ztemp = atomc14 * xdicsur * ztemp 117 114 CALL iom_put( "C14Inv", ztemp ) ! Radiocarbon ocean inventory [10^26 atoms] … … 130 127 #endif 131 128 129 !! * Substitutions 130 # include "do_loop_substitute.h90" 132 131 !!---------------------------------------------------------------------- 133 132 !! NEMO/TOP 4.0 , NEMO Consortium (2018)
Note: See TracChangeset
for help on using the changeset viewer.