source: codes/icosagcm/devel/src/kernels_hex/caldyn_wflux.k90 @ 624

Last change on this file since 624 was 624, checked in by dubos, 7 years ago

devel : introduced kernels for caldyn_vert (optional)
File size: 853 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- caldyn_wflux ----------------------------------
3   DO l = llm-1,1,-1
4      DO ij=ij_omp_begin,ij_omp_end
5         ! cumulate mass flux convergence from top to bottom
6         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
7      END DO
8   END DO
9   DO ij=ij_omp_begin,ij_omp_end
10      dmass_col(ij) = convm(ij,1)
11   END DO
12   DO l = 2,llm
13      DO ij=ij_omp_begin,ij_omp_end
14         ! Compute vertical mass flux (l=1,llm+1 set to zero at init)
15         wflux(ij,l) = mass_bl(ij,l) * dmass_col(ij) - convm(ij,l)
16      END DO
17   END DO
18   ! make sure wflux is up to date
19   !$OMP BARRIER
20   !---------------------------- caldyn_wflux ----------------------------------
21   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.