Wmig

From GHER

Jump to: navigation, search

vertical migration velocity for scalar (positive upwards)

added to velocity field before to be computed in solvhc

im main.f90

CALL SCALPD(PRODI,DESTI,PRODII,DESTII,FLCS,FLCB,WMIGTL,NS)


! DE MIGRATION WMIG:
!   AJOUTTER WMIG*L  ( L DE LA DEFINITION DE COORD VERTICALE 1 ICI)
!   A HU3
!$omp do private(i,j,k)
        do k=3,KSUPI
         do j=2,jmaxm1
          do i=2,imaxm1
          HU3I(I,J,K)=HU3I(I,J,K)+WMIGTL
          enddo
         enddo
        enddo
!$omp do private(i,j,k)
        do k=3,KSUPII
         do j=2,jmaxm1
          do i=2,imaxm1
          HU3II(I,J,K)=HU3II(I,J,K)+WMIGTL
          enddo
         enddo
        enddo
!$omp do private(i,j,k)
         do j=2,jmax-1
          do i=2,imax-1
          HU3I(I,J,2)=HU3I(I,J,2)+WMIGTL*ICII(I,J)
          HU3II(I,J,KMAXII)=HU3II(I,J,KMAXII)+WMIGTL*ICII(I,J)
          enddo
         enddo
      CALL SOLVHC(HCI(1,1,1,NS),HCII(1,1,1,NS),FLCS,FLCB, &
                       PRODI,DESTI, &
                       PRODII,DESTII,NS)

....


!
! REMISE DES BONNES VITESSES VERTICALES
!$omp do private(i,j,k)
        do k=3,KSUPI
         do j=2,jmaxm1
          do i=2,imaxm1
          HU3I(I,J,K)=HU3I(I,J,K)-WMIGTL
          enddo
         enddo
        enddo
!$omp do private(i,j,k)
        do k=3,KSUPII
         do j=2,jmaxm1
          do i=2,imaxm1
          HU3II(I,J,K)=HU3II(I,J,K)-WMIGTL
          enddo
         enddo
        enddo
!$omp do private(i,j)
         do j=2,jmax-1
          do i=2,imax-1
          HU3I(I,J,2)=HU3I(I,J,2)-WMIGTL*ICII(I,J)
          HU3II(I,J,KMAXII)=HU3II(I,J,KMAXII)-WMIGTL*ICII(I,J)
          enddo
         enddo
! END TRACERS
Personal tools