Changeset 2056 for LMDZ5/branches/testing/libf/dyn3dpar/leapfrog_p.F
- Timestamp:
- Jun 11, 2014, 3:46:46 PM (10 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 1998,2000-2023,2025-2029,2032,2034,2036-2049,2051-2055
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/dyn3dpar/leapfrog_p.F
r1999 r2056 8 8 & time_0) 9 9 10 use exner_hyb_m, only: exner_hyb 11 use exner_milieu_m, only: exner_milieu 12 use exner_hyb_p_m, only: exner_hyb_p 13 use exner_milieu_p_m, only: exner_milieu_p 10 14 USE misc_mod 11 15 USE parallel_lmdz … … 149 153 character*10 string10 150 154 151 REAL,SAVE :: alpha(ip1jmp1,llm),beta(ip1jmp1,llm)152 155 REAL,SAVE :: flxw(ip1jmp1,llm) ! flux de masse verticale 153 156 … … 209 212 lafin=.false. 210 213 211 itaufin = nday*day_step 214 if (nday>=0) then 215 itaufin = nday*day_step 216 else 217 itaufin = -nday 218 endif 219 212 220 itaufinp1 = itaufin +1 213 221 … … 241 249 CALL pression ( ip1jmp1, ap, bp, ps, p ) 242 250 if (pressure_exner) then 243 CALL exner_hyb( ip1jmp1, ps, p, alpha,beta,pks, pk, pkf )251 CALL exner_hyb( ip1jmp1, ps, p, pks, pk, pkf ) 244 252 else 245 CALL exner_milieu( ip1jmp1, ps, p, beta,pks, pk, pkf )253 CALL exner_milieu( ip1jmp1, ps, p, pks, pk, pkf ) 246 254 endif 247 255 c$OMP END MASTER … … 705 713 c$OMP BARRIER 706 714 if (pressure_exner) then 707 CALL exner_hyb_p( ip1jmp1, ps, p, alpha,beta,pks, pk, pkf )715 CALL exner_hyb_p( ip1jmp1, ps, p,pks, pk, pkf ) 708 716 else 709 CALL exner_milieu_p( ip1jmp1, ps, p, beta,pks, pk, pkf )717 CALL exner_milieu_p( ip1jmp1, ps, p, pks, pk, pkf ) 710 718 endif 719 ! Appel a geopot ajoute le 2014/05/08 pour garantir la convergence numerique 720 ! avec dyn3dmem 721 CALL geopot ( ip1jmp1, teta , pk , pks, phis , phi ) 711 722 c$OMP BARRIER 712 723 jD_cur = jD_ref + day_ini - day_ref … … 918 929 c$OMP BARRIER 919 930 if (pressure_exner) then 920 CALL exner_hyb_p(ip1jmp1,ps,p, alpha,beta,pks,pk,pkf)931 CALL exner_hyb_p(ip1jmp1,ps,p,pks,pk,pkf) 921 932 else 922 CALL exner_milieu_p(ip1jmp1,ps,p, beta,pks,pk,pkf)933 CALL exner_milieu_p(ip1jmp1,ps,p,pks,pk,pkf) 923 934 endif 924 935 c$OMP BARRIER … … 1059 1070 c$OMP BARRIER 1060 1071 if (pressure_exner) then 1061 CALL exner_hyb_p( ip1jmp1, ps, p, alpha,beta,pks, pk, pkf )1072 CALL exner_hyb_p( ip1jmp1, ps, p, pks, pk, pkf ) 1062 1073 else 1063 CALL exner_milieu_p( ip1jmp1, ps, p, beta,pks, pk, pkf )1074 CALL exner_milieu_p( ip1jmp1, ps, p, pks, pk, pkf ) 1064 1075 endif 1065 1076 c$OMP BARRIER
Note: See TracChangeset
for help on using the changeset viewer.