Changeset 5433 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90
- Timestamp:
- 2015-06-18T08:37:20+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90
r5260 r5433 8 8 USE oce ! ocean dynamics and tracers variables 9 9 USE dom_oce ! ocean space and time domain 10 USE diainsitutem, ONLY: insitu_t, theta2t10 USE diainsitutem, ONLY: rinsitu_t, theta2t 11 11 USE in_out_manager ! I/O units 12 12 USE iom ! I/0 library … … 28 28 29 29 !! * variables for calculating 25-hourly means 30 REAL(wp),SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: tn_25h , sn_25h, insitu_t_25h30 REAL(wp),SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: tn_25h , sn_25h, rinsitu_t_25h 31 31 REAL(wp),SAVE, ALLOCATABLE, DIMENSION(:,:) :: sshn_25h 32 32 REAL(wp),SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: un_25h , vn_25h , wn_25h … … 36 36 #endif 37 37 #if defined key_zdfgls 38 REAL(wp),SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: mxln_25h38 REAL(wp),SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: rmxln_25h 39 39 #endif 40 40 INTEGER, SAVE :: cnt_25h ! Counter for 25 hour means … … 94 94 CALL ctl_stop( 'dia_25h: unable to allocate sn_25h' ) ; RETURN 95 95 ENDIF 96 ALLOCATE( insitu_t_25h(jpi,jpj,jpk), STAT=ierror )97 IF( ierror > 0 ) THEN 98 CALL ctl_stop( 'dia_25h: unable to allocate insitu_t_25h' ) ; RETURN96 ALLOCATE( rinsitu_t_25h(jpi,jpj,jpk), STAT=ierror ) 97 IF( ierror > 0 ) THEN 98 CALL ctl_stop( 'dia_25h: unable to allocate rinsitu_t_25h' ) ; RETURN 99 99 ENDIF 100 100 ALLOCATE( un_25h(jpi,jpj,jpk), STAT=ierror ) … … 125 125 #endif 126 126 # if defined key_zdfgls 127 ALLOCATE( mxln_25h(jpi,jpj,jpk), STAT=ierror )128 IF( ierror > 0 ) THEN 129 CALL ctl_stop( 'dia_25h: unable to allocate mxln_25h' ) ; RETURN127 ALLOCATE( rmxln_25h(jpi,jpj,jpk), STAT=ierror ) 128 IF( ierror > 0 ) THEN 129 CALL ctl_stop( 'dia_25h: unable to allocate rmxln_25h' ) ; RETURN 130 130 ENDIF 131 131 #endif … … 141 141 sn_25h(:,:,:) = tsb(:,:,:,jp_sal) 142 142 CALL theta2t 143 insitu_t_25h(:,:,:) =insitu_t(:,:,:)143 rinsitu_t_25h(:,:,:) = rinsitu_t(:,:,:) 144 144 sshn_25h(:,:) = sshb(:,:) 145 145 un_25h(:,:,:) = ub(:,:,:) … … 152 152 #endif 153 153 # if defined key_zdfgls 154 mxln_25h(:,:,:) = mxln(:,:,:)154 rmxln_25h(:,:,:) = mxln(:,:,:) 155 155 #endif 156 156 #if defined key_lim3 || defined key_lim2 … … 200 200 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d ! temporary workspace 201 201 REAL(wp), DIMENSION(jpi,jpj,3) :: zwtmb ! temporary workspace 202 INTEGER :: nyear0, nmonth0,nday0 ! start year,month,day202 INTEGER :: iyear0, nimonth0,iday0 ! start year,imonth,day 203 203 204 204 !!---------------------------------------------------------------------- … … 235 235 sn_25h(:,:,:) = sn_25h(:,:,:) + tsn(:,:,:,jp_sal) 236 236 CALL theta2t 237 insitu_t_25h(:,:,:) = insitu_t_25h(:,:,:) +insitu_t(:,:,:)237 rinsitu_t_25h(:,:,:) = rinsitu_t_25h(:,:,:) + rinsitu_t(:,:,:) 238 238 sshn_25h(:,:) = sshn_25h(:,:) + sshn (:,:) 239 239 un_25h(:,:,:) = un_25h(:,:,:) + un(:,:,:) … … 246 246 #endif 247 247 # if defined key_zdfgls 248 mxln_25h(:,:,:) =mxln_25h(:,:,:) + mxln(:,:,:)248 rmxln_25h(:,:,:) = rmxln_25h(:,:,:) + mxln(:,:,:) 249 249 #endif 250 250 cnt_25h = cnt_25h + 1 … … 266 266 tn_25h(:,:,:) = tn_25h(:,:,:) / 25.0_wp 267 267 sn_25h(:,:,:) = sn_25h(:,:,:) / 25.0_wp 268 insitu_t_25h(:,:,:) =insitu_t_25h(:,:,:) / 25.0_wp268 rinsitu_t_25h(:,:,:) = rinsitu_t_25h(:,:,:) / 25.0_wp 269 269 sshn_25h(:,:) = sshn_25h(:,:) / 25.0_wp 270 270 un_25h(:,:,:) = un_25h(:,:,:) / 25.0_wp … … 277 277 #endif 278 278 # if defined key_zdfgls 279 mxln_25h(:,:,:) =mxln_25h(:,:,:) / 25.0_wp279 rmxln_25h(:,:,:) = rmxln_25h(:,:,:) / 25.0_wp 280 280 #endif 281 281 … … 286 286 CALL iom_put("temper25h", zw3d) ! potential temperature 287 287 CALL theta2t ! calculate insitu temp 288 zw3d(:,:,:) = insitu_t_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))288 zw3d(:,:,:) = rinsitu_t_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 289 289 CALL iom_put("tempis25h", zw3d) ! in-situ temperature 290 290 zw3d(:,:,:) = sn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) … … 311 311 #endif 312 312 #if defined key_zdfgls 313 zw3d(:,:,:) = mxln_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))313 zw3d(:,:,:) = rmxln_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 314 314 CALL iom_put( "mxln25h",zw3d) 315 315 #endif … … 319 319 sn_25h(:,:,:) = tsn(:,:,:,jp_sal) 320 320 CALL theta2t 321 insitu_t_25h(:,:,:) =insitu_t(:,:,:)321 rinsitu_t_25h(:,:,:) = rinsitu_t(:,:,:) 322 322 sshn_25h(:,:) = sshn (:,:) 323 323 un_25h(:,:,:) = un(:,:,:) … … 330 330 #endif 331 331 # if defined key_zdfgls 332 mxln_25h(:,:,:) = mxln(:,:,:)332 rmxln_25h(:,:,:) = mxln(:,:,:) 333 333 #endif 334 334 cnt_25h = 1
Note: See TracChangeset
for help on using the changeset viewer.