Changeset 5433 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diainsitutem.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/diainsitutem.F90
r5260 r5433 22 22 PUBLIC insitu_tem_alloc ! routines called by step.F90 23 23 24 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) :: insitu_t24 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) :: rinsitu_t 25 25 26 26 !! * Substitutions … … 36 36 ierr = 0 37 37 ! 38 ALLOCATE( insitu_t(jpi,jpj,jpk), STAT=ierr(1) )38 ALLOCATE( rinsitu_t(jpi,jpj,jpk), STAT=ierr(1) ) 39 39 ! 40 40 insitu_tem_alloc = MAXVAL(ierr) … … 54 54 SUBROUTINE theta2t() 55 55 56 INTEGER, PARAMETER :: num_steps=10 ! number of steps in integration57 INTEGER :: step ! iteration counter56 INTEGER, PARAMETER :: inum_steps=10 ! number of steps in integration 57 INTEGER :: jstep ! iteration counter 58 58 INTEGER :: ji, jj, jk ! loop indices 59 59 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zP ! pressure (decibars) … … 76 76 DO ji = 1, jpi 77 77 ! These loops expanded for case where fsdept may be 1D 78 zDP(ji,jj,jk) = fsdept(ji,jj,jk) / real( num_steps)78 zDP(ji,jj,jk) = fsdept(ji,jj,jk) / real(inum_steps) 79 79 END DO 80 80 END DO … … 91 91 zTB(:,:,:) = zT(:,:,:) - zLAPSE(:,:,:) * zDP(:,:,:) 92 92 93 interation: DO step=1,num_steps93 interation: DO jstep=1, inum_steps 94 94 ! Calculate lapse rate (dT/dP) and hence TA 95 95 CALL ATG(zP, zT, zSS, zLAPSE) … … 102 102 END DO interation 103 103 104 insitu_t(:,:,:) = zT(:,:,:) * tmask(:,:,:)105 CALL lbc_lnk( insitu_t, 'T', 1.0)104 rinsitu_t(:,:,:) = zT(:,:,:) * tmask(:,:,:) 105 CALL lbc_lnk( rinsitu_t, 'T', 1.0) 106 106 107 107 END SUBROUTINE theta2t
Note: See TracChangeset
for help on using the changeset viewer.