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

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

devel : introduced kernels for caldyn_vert (optional)

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