Changeset 9190 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

Timestamp:

2018-01-06T15:18:23+01:00 (6 years ago)

Author:

gm

Message:

dev_merge_2017: OPA_SRC: style only, results unchanged

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90 (modified) (7 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

@@ -1 +1 @@
 MODULE icbclv
-
    !!======================================================================
    !!                       ***  MODULE  icbclv  ***
@@ -11 +10 @@
    !!            -    !  2011-05  (Alderson)       budgets into separate module
    !!----------------------------------------------------------------------
+
    !!----------------------------------------------------------------------
    !!   icb_clv_flx   : transfer input flux of ice into iceberg classes
@@ -45 +45 @@
       !!
       !!----------------------------------------------------------------------
-      INTEGER, INTENT(in)             :: kt
+      INTEGER, INTENT(in) ::   kt
       !
-      REAL(wp)                        :: zcalving_used, zdist, zfact
-      INTEGER                         :: jn, ji, jj                    ! loop counters
-      INTEGER                         :: imx                           ! temporary integer for max berg class
-      LOGICAL, SAVE                   :: ll_first_call = .TRUE.
+      REAL(wp)      ::   zcalving_used, zdist, zfact
+      INTEGER       ::   jn, ji, jj                    ! loop counters
+      INTEGER       ::   imx                           ! temporary integer for max berg class
+      LOGICAL, SAVE ::   ll_first_call = .TRUE.
       !!----------------------------------------------------------------------
       !
@@ -70 +70 @@
          DO jj = 2, jpjm1
             DO ji = 2, jpim1
-               IF( berg_grid%calving(ji,jj) /= 0._wp )                                  &    ! Need units of J
-                  berg_grid%stored_heat(ji,jj) = SUM( berg_grid%stored_ice(ji,jj,:) ) *         &  ! initial stored ice in kg
-                                         berg_grid%calving_hflx(ji,jj) * e1e2t(ji,jj) /   &  ! J/s/m2 x m^2 = J/s
-                                         berg_grid%calving(ji,jj)                            ! /calving in kg/s
+               IF( berg_grid%calving(ji,jj) /= 0._wp )                                          &    ! Need units of J
+                  berg_grid%stored_heat(ji,jj) = SUM( berg_grid%stored_ice(ji,jj,:) ) *         &    ! initial stored ice in kg
+                     &                   berg_grid%calving_hflx(ji,jj) * e1e2t(ji,jj) / berg_grid%calving(ji,jj)   ! J/s/m2 x m^2 
+                     !                                                                                             ! = J/s/calving in kg/s
             END DO
          END DO
@@ -80 +80 @@
       ! assume that all calving flux must be distributed even if distribution array does not sum
       ! to one - this may not be what is intended, but it's what you've got
-      DO jj = 1,jpj
-         DO ji = 1,jpi
+      DO jj = 1, jpj
+         DO ji = 1, jpi
             imx = berg_grid%maxclass(ji,jj)
             zdist = SUM( rn_distribution(1:nclasses) ) / SUM( rn_distribution(1:imx) )
             DO jn = 1, imx
-               berg_grid%stored_ice(ji,jj,jn) = berg_grid%stored_ice(ji,jj,jn) + &
-                                          berg_dt * berg_grid%calving(ji,jj) * rn_distribution(jn) * zdist
+               berg_grid%stored_ice(ji,jj,jn) = berg_grid%stored_ice(ji,jj,jn)     &
+                  &                           + berg_dt * berg_grid%calving(ji,jj) * rn_distribution(jn) * zdist
             END DO
          END DO
@@ -98 +98 @@
       !
    END SUBROUTINE icb_clv_flx
+
 
    SUBROUTINE icb_clv()
@@ -171 +172 @@
       END DO
       !
-      DO jn = 1,nclasses
+      DO jn = 1, nclasses
          CALL lbc_lnk( berg_grid%stored_ice(:,:,jn), 'T', 1._wp )
       END DO