Changeset 6517

Timestamp:

2016-05-11T12:12:20+02:00 (8 years ago)

Author:

rfurner

Message:

small fixes to surge flux module

File:

• branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcsurge.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcsurge.F90

@@ -122 +122 @@
 
       !                                            ! compute the surface ocean fluxes using CORE bulk formulea
-      IF( MOD( kt - 1, nn_fsbc ) == 0 )   CALL surge_oce( kt, sf, sst_m, ssu_m, ssv_m, rn_charn_const )
+      IF( MOD( kt - 1, nn_fsbc ) == 0 )   CALL surge_oce( kt, sf, ssu_m, ssv_m, rn_charn_const )
       !
    END SUBROUTINE sbc_surge
    
    
-   SUBROUTINE surge_oce( kt, sf, pst, pu, pv, rn_charn_const )
+   SUBROUTINE surge_oce( kt, sf, pu, pv, rn_charn_const )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE surge_oce  ***
@@ -146 +146 @@
       INTEGER  , INTENT(in   )                 ::   kt    ! time step index
       TYPE(fld), INTENT(inout), DIMENSION(:)   ::   sf    ! input data
-      REAL(wp) , INTENT(in)   , DIMENSION(:,:) ::   pst   ! surface temperature                      [Celcius]
       REAL(wp) , INTENT(in)   , DIMENSION(:,:) ::   pu    ! surface current at U-point (i-component) [m/s]
       REAL(wp) , INTENT(in)   , DIMENSION(:,:) ::   pv    ! surface current at V-point (j-component) [m/s]
@@ -199 +198 @@
             !Iterate
             DO WHILE((abs(Cd(ji,jj)-z_Cd1))>1E-6)
-            z_Cd1=Cd(ji,jj)
-            z_z0=rn_charn_const*z_Cd1*wndm(ji,jj)**2/grav
-            Cd(ji,jj)=(0.41_wp/log(10._wp/z_z0))**2
-            zi=zi+1
+               z_Cd1=Cd(ji,jj)
+               z_z0=rn_charn_const*z_Cd1*wndm(ji,jj)**2/grav
+               Cd(ji,jj)=(0.41_wp/log(10._wp/z_z0))**2
+               zi=zi+1
             ENDDO
          ENDDO