Changeset 4188 for branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90

Timestamp:

2013-11-14T10:40:55+01:00 (11 years ago)

Author:

agn

Message:

Divide zpsi_[uv]w by box width to give units in m^{2/s, correct grid for psiv_mle, update namelists}

File:

• branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90 (modified) (1 diff)

branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90

@@ -245 +245 @@
          CALL iom_put( "Lf_NHpf" , zLf_NH  )    ! Lf = N H / f
          !
-         zpsim_u(jpi,:) = 0._wp   ;   zpsim_u(:,jpj) = 0._wp      ! set a value where it is undefined
-         zpsim_v(jpi,:) = 0._wp   ;   zpsim_v(:,jpj) = 0._wp      ! before the output
+         ! divide by cross distance to give streamfunction with dimensions m^2/s
+         DO jk = 1, ikmax+1
+            zpsi_uw(:,:,jk) = zpsi_uw(:,:,jk)/e2u(:,:)
+            zpsi_vw(:,:,jk) = zpsi_vw(:,:,jk)/e1v(:,:)
+         END DO
          CALL iom_put( "psiu_mle", zpsi_uw )    ! i-mle streamfunction
          CALL iom_put( "psiv_mle", zpsi_vw )    ! j-mle streamfunction