Changeset 3957


Ignore:
Timestamp:
2013-07-05T16:30:37+02:00 (7 years ago)
Author:
acc
Message:

Branch 2013/dev_r3858_NOC_ZTC, #863. Fixes to tidal potential forcing changes to ensure no numerical changes due to the re-styling

Location:
branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90

    r3953 r3957  
    147147               DO ji = fs_2, fs_jpim1   ! vector opt. 
    148148                  spgu(ji,jj) = spgu(ji,jj) + ( zpice(ji+1,jj) - zpice(ji,jj) ) / e1u(ji,jj) 
    149                   spgv(ji,jj) = spgu(ji,jj) + ( zpice(ji,jj+1) - zpice(ji,jj) ) / e2v(ji,jj) 
     149                  spgv(ji,jj) = spgv(ji,jj) + ( zpice(ji,jj+1) - zpice(ji,jj) ) / e2v(ji,jj) 
    150150               END DO 
    151151            END DO 
  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90

    r3955 r3957  
    4949      !!---------------------------------------------------------------------- 
    5050 
    51       IF( kt /= nit000 )   CALL tide_init( kt ) 
    52  
    53       IF( nsec_day == NINT(0.5 * rdttra(1)) ) THEN      ! start a new day 
     51      IF( nsec_day == NINT(0.5_wp * rdttra(1)) ) THEN      ! start a new day 
    5452         ! 
    5553         IF( kt == nit000 ) THEN 
    5654            ALLOCATE( amp_pot(jpi,jpj,nb_harmo),                      & 
    5755               &      phi_pot(jpi,jpj,nb_harmo), pot_astro(jpi,jpj)   ) 
    58             ! 
    59             amp_pot(:,:,:) = 0._wp 
    60             phi_pot(:,:,:) = 0._wp 
    61             pot_astro(:,:) = 0._wp 
    6256         ENDIF 
     57         ! 
     58         amp_pot(:,:,:) = 0._wp 
     59         phi_pot(:,:,:) = 0._wp 
     60         pot_astro(:,:) = 0._wp 
    6361         ! 
    6462         CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo ) 
     
    9189 
    9290      DO jk = 1, nb_harmo 
    93          zcons = 0.7 * Wave(ntide(jk))%equitide * ftide(jk) 
     91         zcons = 0.7_wp * Wave(ntide(jk))%equitide * ftide(jk) 
    9492         DO ji = 1, jpi 
    9593            DO jj = 1, jpj 
     
    10098               ztmp = v0tide(jk) + utide(jk) + Wave(ntide(jk))%nutide * zlon 
    10199               ! le potentiel est composé des effets des astres: 
    102                IF( Wave(ntide(jk))%nutide == 1 )   zcs = zcons * SIN( 2.*zlat ) 
    103                IF( Wave(ntide(jk))%nutide == 2 )   zcs = zcons * COS( zlat )**2 
     100               IF    ( Wave(ntide(jk))%nutide == 1 )  THEN  ;  zcs = zcons * SIN( 2._wp*zlat ) 
     101               ELSEIF( Wave(ntide(jk))%nutide == 2 )  THEN  ;  zcs = zcons * COS( zlat )**2 
     102               ELSE                                         ;  zcs = 0._wp 
     103               ENDIF 
    104104               ztmp1 = ztmp1 + zcs * COS( ztmp ) 
    105105               ztmp2 = ztmp2 - zcs * SIN( ztmp ) 
    106106               zamp = SQRT( ztmp1*ztmp1 + ztmp2*ztmp2 ) 
    107107               amp_pot(ji,jj,jk) = zamp 
    108                phi_pot(ji,jj,jk) = ATAN2( -ztmp2 / MAX( 1.e-10 , zamp ) ,   & 
    109                   &                        ztmp1 / MAX( 1.e-10,  zamp )   ) 
     108               phi_pot(ji,jj,jk) = ATAN2( -ztmp2 / MAX( 1.e-10_wp , zamp ) ,   & 
     109                  &                        ztmp1 / MAX( 1.e-10_wp,  zamp )   ) 
    110110            END DO 
    111111         END DO 
  • branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90

    r3953 r3957  
    104104      ENDIF 
    105105      ! 
    106       IF( nsec_day == NINT(0.5 * rdttra(1)) ) THEN 
    107          ! 
    108          CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo ) 
    109          ! 
    110          kt_tide = kt 
    111          ! 
    112          IF(lwp) THEN 
    113             WRITE(numout,*) 
    114             WRITE(numout,*) 'tide_ini : Update of the tidal components at kt=', kt 
    115             WRITE(numout,*) '~~~~~~~~ ' 
    116             DO jk = 1, nb_harmo 
    117                WRITE(numout,*) Wave(ntide(jk))%cname_tide, utide(jk), ftide(jk), v0tide(jk), omega_tide(jk) 
    118             END DO 
    119          ENDIF 
    120          ! 
    121       ENDIF 
    122       ! 
    123106   END SUBROUTINE tide_init 
    124107      
Note: See TracChangeset for help on using the changeset viewer.