Changeset 3957
- Timestamp:
- 2013-07-05T16:30:37+02:00 (11 years ago)
- 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 147 147 DO ji = fs_2, fs_jpim1 ! vector opt. 148 148 spgu(ji,jj) = spgu(ji,jj) + ( zpice(ji+1,jj) - zpice(ji,jj) ) / e1u(ji,jj) 149 spgv(ji,jj) = spg u(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) 150 150 END DO 151 151 END DO -
branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90
r3955 r3957 49 49 !!---------------------------------------------------------------------- 50 50 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 54 52 ! 55 53 IF( kt == nit000 ) THEN 56 54 ALLOCATE( amp_pot(jpi,jpj,nb_harmo), & 57 55 & phi_pot(jpi,jpj,nb_harmo), pot_astro(jpi,jpj) ) 58 !59 amp_pot(:,:,:) = 0._wp60 phi_pot(:,:,:) = 0._wp61 pot_astro(:,:) = 0._wp62 56 ENDIF 57 ! 58 amp_pot(:,:,:) = 0._wp 59 phi_pot(:,:,:) = 0._wp 60 pot_astro(:,:) = 0._wp 63 61 ! 64 62 CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo ) … … 91 89 92 90 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) 94 92 DO ji = 1, jpi 95 93 DO jj = 1, jpj … … 100 98 ztmp = v0tide(jk) + utide(jk) + Wave(ntide(jk))%nutide * zlon 101 99 ! 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 104 104 ztmp1 = ztmp1 + zcs * COS( ztmp ) 105 105 ztmp2 = ztmp2 - zcs * SIN( ztmp ) 106 106 zamp = SQRT( ztmp1*ztmp1 + ztmp2*ztmp2 ) 107 107 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 ) ) 110 110 END DO 111 111 END DO -
branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90
r3953 r3957 104 104 ENDIF 105 105 ! 106 IF( nsec_day == NINT(0.5 * rdttra(1)) ) THEN107 !108 CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo )109 !110 kt_tide = kt111 !112 IF(lwp) THEN113 WRITE(numout,*)114 WRITE(numout,*) 'tide_ini : Update of the tidal components at kt=', kt115 WRITE(numout,*) '~~~~~~~~ '116 DO jk = 1, nb_harmo117 WRITE(numout,*) Wave(ntide(jk))%cname_tide, utide(jk), ftide(jk), v0tide(jk), omega_tide(jk)118 END DO119 ENDIF120 !121 ENDIF122 !123 106 END SUBROUTINE tide_init 124 107
Note: See TracChangeset
for help on using the changeset viewer.