Changeset 1520 for LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F
- Timestamp:
- May 23, 2011, 1:37:09 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F
r1505 r1520 164 164 165 165 character*80 dynhist_file, dynhistave_file 166 character *20 modname166 character(len=*),parameter :: modname="leapfrog" 167 167 character*80 abort_message 168 168 … … 208 208 itaufin = nday*day_step 209 209 itaufinp1 = itaufin +1 210 modname="leapfrog_p"211 210 212 211 itau = 0 … … 236 235 dq(:,:,:)=0. 237 236 CALL pression ( ip1jmp1, ap, bp, ps, p ) 238 CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) 237 if (disvert_type==1) then 238 CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) 239 else ! we assume that we are in the disvert_type==2 case 240 CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf ) 241 endif 239 242 c$OMP END MASTER 240 243 c----------------------------------------------------------------------- … … 691 694 692 695 c$OMP BARRIER 693 CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta,pks, pk, pkf ) 696 if (disvert_type==1) then 697 CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta,pks, pk, pkf ) 698 else ! we assume that we are in the disvert_type==2 case 699 CALL exner_milieu_p( ip1jmp1, ps, p, beta, pks, pk, pkf ) 700 endif 694 701 c$OMP BARRIER 695 702 jD_cur = jD_ref + day_ini - day_ref … … 1030 1037 CALL pression_p ( ip1jmp1, ap, bp, ps, p ) 1031 1038 c$OMP BARRIER 1032 CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) 1039 if (disvert_type==1) then 1040 CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) 1041 else ! we assume that we are in the disvert_type==2 case 1042 CALL exner_milieu_p( ip1jmp1, ps, p, beta, pks, pk, pkf ) 1043 endif 1033 1044 c$OMP BARRIER 1034 1045
Note: See TracChangeset
for help on using the changeset viewer.