Changeset 4332 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2013-12-11T15:38:42+01:00 (10 years ago)

Author:

clem

Message:

update LIM3 to fix remaining bugs. Now working in global and regional config.

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90 (modified) (3 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

@@ -49 +49 @@
    PUBLIC   lim_sbc_tau    ! called by sbc_ice_lim
 
-   REAL(wp)  ::   epsi16 = 1.e-16_wp   ! constant values
    REAL(wp)  ::   rzero  = 0._wp    
    REAL(wp)  ::   rone   = 1._wp
@@ -145 +144 @@
 
             !   computation the solar flux at ocean surface
-            IF (lk_cpl) THEN ! be carfeful: not being tested yet
+            IF (lk_cpl) THEN ! be carfeful: not been tested yet
                ! original line
                !zfcm1 = qsr_tot(ji,jj) + fstric(ji,jj) * at_i(ji,jj)
@@ -188 +187 @@
             qns(ji,jj) = zfcm2 - fdtcn(ji,jj)                        ! non solar heat flux
             !                           ! fdtcn : turbulent oceanic heat flux
-
-            !!gm   this IF prevents the vertorisation of the whole loop
-          !  IF ( ( ji == jiindx ) .AND. ( jj == jjindx) ) THEN
-          !     WRITE(numout,*) ' lim_sbc : heat fluxes '
-          !     WRITE(numout,*) ' qsr       : ', qsr(jiindx,jjindx)
-          !     WRITE(numout,*) ' pfrld     : ', pfrld(jiindx,jjindx)
-          !     WRITE(numout,*) ' fstric    : ', fstric (jiindx,jjindx)
-          !     WRITE(numout,*)
-          !     WRITE(numout,*) ' qns       : ', qns(jiindx,jjindx)
-          !     WRITE(numout,*) ' fdtcn     : ', fdtcn(jiindx,jjindx)
-          !     WRITE(numout,*) ' ifral     : ', ifral
-          !     WRITE(numout,*) ' ial       : ', ial  
-          !     WRITE(numout,*) ' qcmif     : ', qcmif(jiindx,jjindx)
-          !     WRITE(numout,*) ' qldif     : ', qldif(jiindx,jjindx)
-          !     !WRITE(numout,*) ' qcmif / dt: ', qcmif(jiindx,jjindx) * r1_rdtice
-          !     !WRITE(numout,*) ' qldif / dt: ', qldif(jiindx,jjindx) * r1_rdtice
-          !     WRITE(numout,*) ' ifrdv     : ', ifrdv
-          !     WRITE(numout,*) ' qfvbq     : ', qfvbq(jiindx,jjindx)
-          !     WRITE(numout,*) ' qdtcn     : ', qdtcn(jiindx,jjindx)
-          !     !WRITE(numout,*) ' qfvbq / dt: ', qfvbq(jiindx,jjindx) * r1_rdtice
-          !     !WRITE(numout,*) ' qdtcn / dt: ', qdtcn(jiindx,jjindx) * r1_rdtice
-          !     WRITE(numout,*) ' '
-          !     WRITE(numout,*) ' fdtcn     : ', fdtcn(jiindx,jjindx)
-          !     WRITE(numout,*) ' fhmec     : ', fhmec(jiindx,jjindx)
-          !     WRITE(numout,*) ' fheat_mec : ', fheat_mec(jiindx,jjindx)
-          !     WRITE(numout,*) ' fhbri     : ', fhbri(jiindx,jjindx)
-          !     WRITE(numout,*) ' fheat_res : ', fheat_res(jiindx,jjindx)
-          !  ENDIF
-            !!gm   end
          END DO
       END DO