Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/leapfrog_loc.F
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/leapfrog_loc.F
r5091 r5099 1 ! 1 2 2 ! $Id$ 3 ! 3 4 4 c 5 5 c … … 256 256 !$OMP MASTER 257 257 if (firstcall) then 258 ! 258 259 259 ! ALLOCATE(p(ijb_u:ije_u,llmp1)) 260 260 ! ALLOCATE(pks(ijb_u:ije_u)) … … 579 579 ! call Register_SwapFieldHallo(finvmaold,finvmaold,ip1jmp1,llm, 580 580 ! & jj_Nb_caldyn,0,0,TestRequest) 581 ! 581 582 582 ! do j=1,nqtot 583 583 ! call Register_SwapFieldHallo(q(:,:,j),q(:,:,j),ip1jmp1,llm, 584 584 ! & jj_nb_caldyn,0,0,TestRequest) 585 585 ! enddo 586 ! 586 587 587 ! call Set_Distrib(distrib_caldyn) 588 588 ! call SendRequest(TestRequest) … … 864 864 865 865 ! c Diagnostique de conservation de l'energie : initialisation 866 ! 866 867 867 ! c-jld 868 868 ! c$OMP BARRIER … … 879 879 ! #endif 880 880 ! call SetTag(Request_physic,800) 881 ! 881 882 882 ! call Register_SwapField_u(ucov,ucov,distrib_physic, 883 883 ! * Request_physic,up=2,down=2) 884 ! 884 885 885 ! call Register_SwapField_v(vcov,vcov,distrib_physic, 886 886 ! * Request_physic,up=2,down=2) … … 888 888 ! call Register_SwapField_u(teta,teta,distrib_physic, 889 889 ! * Request_physic,up=2,down=2) 890 ! 890 891 891 ! call Register_SwapField_u(masse,masse,distrib_physic, 892 892 ! * Request_physic,up=1,down=2) … … 894 894 ! call Register_SwapField_u(p,p,distrib_physic, 895 895 ! * Request_physic,up=2,down=2) 896 ! 896 897 897 ! call Register_SwapField_u(pk,pk,distrib_physic, 898 898 ! * Request_physic,up=2,down=2) 899 ! 899 900 900 ! call Register_SwapField_u(phis,phis,distrib_physic, 901 901 ! * Request_physic,up=2,down=2) 902 ! 902 903 903 ! call Register_SwapField_u(phi,phi,distrib_physic, 904 904 ! * Request_physic,up=2,down=2) 905 ! 905 906 906 ! call Register_SwapField_u(w,w,distrib_physic, 907 907 ! * Request_physic,up=2,down=2) 908 ! 908 909 909 ! call Register_SwapField_u(q,q,distrib_physic, 910 910 ! * Request_physic,up=2,down=2) … … 912 912 ! call Register_SwapField_u(flxw,flxw,distrib_physic, 913 913 ! * Request_physic,up=2,down=2) 914 ! 914 915 915 ! call SendRequest(Request_Physic) 916 916 ! c$OMP BARRIER … … 921 921 ! call Set_Distrib(distrib_Physic) 922 922 ! call VTe(VThallo) 923 ! 923 924 924 ! call VTb(VTphysiq) 925 925 ! c$OMP END MASTER … … 1080 1080 ! call Register_SwapField_u(ucov,ucov, 1081 1081 ! * distrib_caldyn,Request_physic) 1082 ! 1082 1083 1083 ! call Register_SwapField_v(vcov,vcov, 1084 1084 ! * distrib_caldyn,Request_physic) 1085 ! 1085 1086 1086 ! call Register_SwapField_u(teta,teta, 1087 1087 ! * distrib_caldyn,Request_physic) 1088 ! 1088 1089 1089 ! call Register_SwapField_u(masse,masse, 1090 1090 ! * distrib_caldyn,Request_physic) … … 1092 1092 ! call Register_SwapField_u(p,p, 1093 1093 ! * distrib_caldyn,Request_physic) 1094 ! 1094 1095 1095 ! call Register_SwapField_u(pk,pk, 1096 1096 ! * distrib_caldyn,Request_physic) 1097 ! 1097 1098 1098 ! call Register_SwapField_u(phis,phis, 1099 1099 ! * distrib_caldyn,Request_physic) 1100 ! 1100 1101 1101 ! call Register_SwapField_u(phi,phi, 1102 1102 ! * distrib_caldyn,Request_physic) 1103 ! 1103 1104 1104 ! call Register_SwapField_u(w,w, 1105 1105 ! * distrib_caldyn,Request_physic) … … 1107 1107 ! call Register_SwapField_u(q,q, 1108 1108 ! * distrib_caldyn,Request_physic) 1109 ! 1109 1110 1110 ! call SendRequest(Request_Physic) 1111 1111 ! c$OMP BARRIER … … 1226 1226 !c$OMP MASTER 1227 1227 ! call suspend_timer(timer_caldyn) 1228 ! 1228 1229 1229 !c print*,'Entree dans la dissipation : Iteration No ',true_itau 1230 1230 !c calcul de l'energie cinetique avant dissipation … … 1277 1277 ! call WriteField_u('dtetadis',dtetadis) 1278 1278 !#endif 1279 ! 1279 1280 1280 !! CALL FTRACE_REGION_END("dissip") 1281 ! 1281 1282 1282 ! ijb=ij_begin 1283 1283 ! ije=ij_end … … 1318 1318 ! call covcont_loc(llm,ucov,vcov,ucont,vcont) 1319 1319 ! call enercin_loc(vcov,ucov,vcont,ucont,ecin) 1320 ! 1320 1321 1321 ! ijb=ij_begin 1322 1322 ! ije=ij_end … … 1349 1349 ! ijb=ij_begin 1350 1350 ! ije=ij_end 1351 ! 1351 1352 1352 ! if (pole_nord) then 1353 1353 !c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) … … 1369 1369 ! ENDDO 1370 1370 ! tpn = SSUM(iim,tppn,1)/apoln 1371 ! 1371 1372 1372 ! DO ij = 1, iip1 1373 1373 ! ps( ij ) = tpn … … 1375 1375 !c$OMP END MASTER 1376 1376 ! endif 1377 ! 1377 1378 1378 ! if (pole_sud) then 1379 1379 !c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) … … 1395 1395 ! ENDDO 1396 1396 ! tps = SSUM(iim,tpps,1)/apols 1397 ! 1397 1398 1398 ! DO ij = 1, iip1 1399 1399 ! ps(ij+ip1jm) = tps … … 1408 1408 1409 1409 ! call stop_timer(timer_dissip) 1410 ! 1410 1411 1411 ! call VTb(VThallo) 1412 1412 !c$OMP END MASTER … … 1738 1738 ! iday = day_ini+itau/day_step 1739 1739 ! time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0 1740 ! 1740 1741 1741 ! IF(time.GT.1.) THEN 1742 1742 ! time = time-1.
Note: See TracChangeset
for help on using the changeset viewer.