- Timestamp:
- 2020-05-14T21:46:00+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Property svn:externals
-
old new 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@HEAD sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/TRD/trdpen.F90
r10425 r12928 35 35 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: rab_pe ! partial derivatives of PE anomaly with respect to T and S 36 36 37 !! * Substitutions38 # include "vectopt_loop_substitute.h90"39 37 !!---------------------------------------------------------------------- 40 38 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 55 53 56 54 57 SUBROUTINE trd_pen( ptrdx, ptrdy, ktrd, kt, pdt )55 SUBROUTINE trd_pen( ptrdx, ptrdy, ktrd, kt, pdt, Kmm ) 58 56 !!--------------------------------------------------------------------- 59 57 !! *** ROUTINE trd_tra_mng *** … … 66 64 INTEGER , INTENT(in) :: ktrd ! tracer trend index 67 65 INTEGER , INTENT(in) :: kt ! time step index 66 INTEGER , INTENT(in) :: Kmm ! time level index 68 67 REAL(wp) , INTENT(in) :: pdt ! time step [s] 69 68 ! … … 77 76 IF( kt /= nkstp ) THEN ! full eos: set partial derivatives at the 1st call of kt time step 78 77 nkstp = kt 79 CALL eos_pen( ts n, rab_PE, zpe)78 CALL eos_pen( ts(:,:,:,:,Kmm), rab_PE, zpe, Kmm ) 80 79 CALL iom_put( "alphaPE", rab_pe(:,:,:,jp_tem) ) 81 80 CALL iom_put( "betaPE" , rab_pe(:,:,:,jp_sal) ) … … 95 94 IF( ln_linssh ) THEN ! cst volume : adv flux through z=0 surface 96 95 ALLOCATE( z2d(jpi,jpj) ) 97 z2d(:,:) = w n(:,:,1) * ( &98 & - ( rab_n(:,:,1,jp_tem) + rab_pe(:,:,1,jp_tem) ) * ts n(:,:,1,jp_tem) &99 & + ( rab_n(:,:,1,jp_sal) + rab_pe(:,:,1,jp_sal) ) * ts n(:,:,1,jp_sal) &100 & ) / e3t _n(:,:,1)96 z2d(:,:) = ww(:,:,1) * ( & 97 & - ( rab_n(:,:,1,jp_tem) + rab_pe(:,:,1,jp_tem) ) * ts(:,:,1,jp_tem,Kmm) & 98 & + ( rab_n(:,:,1,jp_sal) + rab_pe(:,:,1,jp_sal) ) * ts(:,:,1,jp_sal,Kmm) & 99 & ) / e3t(:,:,1,Kmm) 101 100 CALL iom_put( "petrd_sad" , z2d ) 102 101 DEALLOCATE( z2d ) … … 112 111 CASE ( jptra_bbc ) ; CALL iom_put( "petrd_bbc" , zpe ) ! bottom bound cond (geoth flux) 113 112 CASE ( jptra_atf ) ; CALL iom_put( "petrd_atf" , zpe ) ! asselin time filter (last trend) 114 !IF( ln_linssh ) THEN ! cst volume : ssh term (otherwise include in e3t variation)115 ! ALLOCATE( z2d(jpi,jpj) )116 ! z2d(:,:) = ( ssha(:,:) - sshb(:,:) ) &117 ! & * ( dPE_dt(:,:,1) * tsn(:,:,1,jp_tem) &118 ! & + dPE_ds(:,:,1) * tsn(:,:,1,jp_sal) ) / ( e3t_n(:,:,1) * pdt )119 ! CALL iom_put( "petrd_sad" , z2d )120 ! DEALLOCATE( z2d )121 !ENDIF122 113 ! 123 114 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.