Changeset 12124 for NEMO/branches
- Timestamp:
- 2019-12-09T15:36:08+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/DIA/diaar5.F90
r12009 r12124 30 30 PUBLIC dia_ar5_hst ! heat/salt transport 31 31 32 REAL(wp) 33 REAL(wp) 32 REAL(wp) :: vol0 ! ocean volume (interior domain) 33 REAL(wp) :: area_tot ! total ocean surface (interior domain) 34 34 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,: ) :: area ! cell surface (interior domain) 35 35 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,: ) :: thick0 ! ocean thickness (interior domain) … … 327 327 !! Called from all advection and/or diffusion routines 328 328 !!---------------------------------------------------------------------- 329 INTEGER , INTENT(in ) ::ktra ! tracer index330 CHARACTER(len=3) , INTENT(in) ::cptr ! transport type 'adv'/'ldf'331 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) ::pua ! 3D input array of advection/diffusion332 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) ::pva ! 3D input array of advection/diffusion333 ! 334 INTEGER 335 REAL(wp), DIMENSION(jpi,jpj) ::z2d329 INTEGER , INTENT(in) :: ktra ! tracer index 330 CHARACTER(len=3) , INTENT(in) :: cptr ! transport type 'adv'/'ldf' 331 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) :: pua ! 3D input array of advection/diffusion 332 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) :: pva ! 3D input array of advection/diffusion 333 ! 334 INTEGER :: ji, jj, jk 335 REAL(wp), DIMENSION(jpi,jpj) :: z2d 336 336 337 337 -
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/DIA/diahth.F90
r11995 r12124 303 303 SUBROUTINE dia_hth_dep( ptem, pdept ) 304 304 ! 305 REAL(wp), INTENT(in) ::ptem306 REAL(wp), DIMENSION(jpi,jpj), INTENT(out) :: pdept307 ! 308 INTEGER :: ji, jj, jk, iid309 REAL(wp) :: zztmp, zzdep310 INTEGER, DIMENSION(jpi,jpj) :: iktem305 REAL(wp), INTENT(in) :: ptem 306 REAL(wp), DIMENSION(jpi,jpj), INTENT(out) :: pdept 307 ! 308 INTEGER :: ji, jj, jk, iid 309 REAL(wp) :: zztmp, zzdep 310 INTEGER, DIMENSION(jpi,jpj) :: iktem 311 311 312 312 ! --------------------------------------- ! … … 349 349 SUBROUTINE dia_hth_htc( pdep, ptn, phtc ) 350 350 ! 351 REAL(wp), INTENT(in) ::pdep ! depth over the heat content352 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) ::ptn353 REAL(wp), DIMENSION(jpi,jpj) , INTENT(inout) ::phtc354 ! 355 INTEGER :: ji, jj, jk, ik356 REAL(wp), DIMENSION(jpi,jpj) :: zthick357 INTEGER , DIMENSION(jpi,jpj) :: ilevel351 REAL(wp), INTENT(in) :: pdep ! depth over the heat content 352 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) :: ptn 353 REAL(wp), DIMENSION(jpi,jpj) , INTENT(inout) :: phtc 354 ! 355 INTEGER :: ji, jj, jk, ik 356 REAL(wp), DIMENSION(jpi,jpj) :: zthick 357 INTEGER , DIMENSION(jpi,jpj) :: ilevel 358 358 359 359 -
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/IOM/iom.F90
r11995 r12124 1730 1730 1731 1731 SUBROUTINE iom_p4d( cdname, pfield4d ) 1732 CHARACTER(LEN=*) , INTENT(in) :: cdname1732 CHARACTER(LEN=*) , INTENT(in) :: cdname 1733 1733 REAL(wp), DIMENSION(:,:,:,:), INTENT(in) :: pfield4d 1734 1734 #if defined key_iomput -
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/SBC/sbcblk.F90
r11995 r12124 918 918 IF( iom_use('evap_ao_cea') .OR. iom_use('hflx_evap_cea') ) THEN 919 919 ztmp(:,:) = zevap(:,:) * ( 1._wp - at_i_b(:,:) ) 920 IF( iom_use('evap_ao_cea' ) )CALL iom_put( 'evap_ao_cea' , ztmp(:,:) * tmask(:,:,1) ) ! ice-free oce evap (cell average)921 IF( iom_use('hflx_evap_cea') )CALL iom_put( 'hflx_evap_cea', ztmp(:,:) * sst_m(:,:) * rcp * tmask(:,:,1) ) ! heat flux from evap (cell average)920 CALL iom_put( 'evap_ao_cea' , ztmp(:,:) * tmask(:,:,1) ) ! ice-free oce evap (cell average) 921 CALL iom_put( 'hflx_evap_cea', ztmp(:,:) * sst_m(:,:) * rcp * tmask(:,:,1) ) ! heat flux from evap (cell average) 922 922 ENDIF 923 923 IF( iom_use('hflx_rain_cea') ) THEN 924 924 ztmp(:,:) = rcp * ( SUM( (ptsu-rt0) * a_i_b, dim=3 ) + sst_m(:,:) * ( 1._wp - at_i_b(:,:) ) ) 925 IF( iom_use('hflx_rain_cea') )CALL iom_put( 'hflx_rain_cea', ( tprecip(:,:) - sprecip(:,:) ) * ztmp(:,:) ) ! heat flux from rain (cell average)925 CALL iom_put( 'hflx_rain_cea', ( tprecip(:,:) - sprecip(:,:) ) * ztmp(:,:) ) ! heat flux from rain (cell average) 926 926 ENDIF 927 927 IF( iom_use('hflx_snow_cea') .OR. iom_use('hflx_snow_ao_cea') .OR. iom_use('hflx_snow_ai_cea') ) THEN … … 930 930 ENDWHERE 931 931 ztmp2(:,:) = sprecip(:,:) * ( ztmp(:,:) - rLfus ) 932 IF( iom_use('hflx_snow_cea') )CALL iom_put('hflx_snow_cea' , ztmp2(:,:) ) ! heat flux from snow (cell average)933 IF( iom_use('hflx_snow_ao_cea') )CALL iom_put('hflx_snow_ao_cea', ztmp2(:,:) * ( 1._wp - zsnw(:,:) ) ) ! heat flux from snow (over ocean)934 IF( iom_use('hflx_snow_ai_cea') )CALL iom_put('hflx_snow_ai_cea', ztmp2(:,:) * zsnw(:,:) ) ! heat flux from snow (over ice)932 CALL iom_put('hflx_snow_cea' , ztmp2(:,:) ) ! heat flux from snow (cell average) 933 CALL iom_put('hflx_snow_ao_cea', ztmp2(:,:) * ( 1._wp - zsnw(:,:) ) ) ! heat flux from snow (over ocean) 934 CALL iom_put('hflx_snow_ai_cea', ztmp2(:,:) * zsnw(:,:) ) ! heat flux from snow (over ice) 935 935 ENDIF 936 936 !
Note: See TracChangeset
for help on using the changeset viewer.