Changeset 22
- Timestamp:
- 11/25/15 20:33:43 (8 years ago)
- Location:
- trunk/SOURCES
- Files:
-
- 39 added
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/Fichiers-parametres/Makefile.tof-lsce3130.inc
r14 r22 20 20 # lit le nom du fichier param pour lancer sur r2d2 par une commande : echo jobparam | LBQ.... 21 21 22 JOB = 'job' 22 JOB = 'job' 23 23 24 24 … … 31 31 #MKL_LIB = -L$(DIRLOCAL)/intel/mkl9/lib/32 -lmkl_lapack -lmkl_ia32 -lguide -lpthread 32 32 #MKL_LIB = -L$(DIRLOCAL)/intel/mkl10/lib/32 -lmkl -lguide -lpthread # -lmklapack 33 MKL_LIB = -L$MKLROOT/lib/em64t -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -liomp5 -lpthread 34 export $MKL_LIB 33 34 ifeq ($(mkl_c), 1) 35 MKL_LIB = -L$MKLROOT/lib/em64t -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -liomp5 -lpthread 36 export $MKL_LIB 37 endif 35 38 36 39 IFORT= ifort … … 47 50 F_NETCDF = $(IFORT) $(ARITHM) -c -I$(NCDF_INC) # -traceback -CB #-g #-pg -ipo # -g 48 51 FT = $(IFORT) $(ARITHM) -c -I$(NCDF_INC) # -traceback -CB #-g #-pg -ipo # -g 52 F_90 = -w90 49 53 50 54 -
trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_mod.f90
r15 r22 87 87 namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file 88 88 89 ! dmr moved up since namelist declaration should be in the header 90 91 namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time 92 89 93 428 format(A) 90 94 rewind(num_param) ! pour revenir au debut du fichier param_list.dat … … 145 149 !lecture namelist 146 150 147 namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time148 151 149 152 rewind(num_param) ! pour revenir au debut du fichier param_list.dat -
trunk/SOURCES/Makefile
r5 r22 5 5 #------------------------------------------- 6 6 7 include Fichiers-parametres/Makefile.tof-lsce3130.inc 8 # include Fichiers-parametres/Makefile.catritz-pc226.inc9 # include Fichiers-parametres/Makefile.pc246.inc7 #------------------------------------------- 8 # dmr: Added options for a more flexible Makefile with the two compilers inbedded 9 #------------------------------------------- 10 10 11 # Choice for compiler (ifort = 1, gfortran = 0) 12 ifort ?= 0 13 14 # Choice for librairies (mkl_c = 1 = MKL, = 0 = BLAS) 15 #Â NOTA: mkl_c = 1 with ifort = 0 ignores the MKL and uses BLAS anyhow 16 mkl_c ?= 1 17 18 ifeq ($(ifort),1) 19 include Fichiers-parametres/Makefile.tof-lsce3130.inc 20 else 21 include Fichiers-parametres/Makefile.dmr-lsce3027.inc 22 endif 11 23 12 24 # Compilation (le corps du Makefile) 13 25 #------------------------------------------- 14 include Makefile.grisli.inc 26 27 ifeq ($(ifort),1) 28 include Makefile.grisli.inc 29 else 30 include Makefile.grisli-gfortran.inc 31 endif -
trunk/SOURCES/Makefile.grisli.inc
r18 r22 11 11 #JOB = 'job' 12 12 13 # pour lancer 13 # pour lancer 14 14 # sur pc: nohup `echo job-test | LBqAntar-15` & ou echo job-test | LBqAntar-15 15 15 # sur r2d2 le mettre dans job*sh … … 24 24 mod_ell = remplimat_declar_tabTu_mod.o eq_ellipt_sgbsv_mod-0.2.o 25 25 routine_elliptiques = remplimat-shelves-tabTu.o graphique_L2.o graph-sgbsv.o 26 diagnoshelf = diagno-L2_mod.o 26 diagnoshelf = diagno-L2_mod.o 27 27 #diagnoshelf = fake_diagno_mod.o 28 28 … … 122 122 123 123 # Toy est la version juin 2012 avec Schooffing 124 124 125 toy_recul = toy_declar_3D_grisli_mod.o \ 125 126 declar_io_ncdf_toy_mod_fake.o\ … … 261 262 module_choix-gsearise.o track_gsearise_mod.o 262 263 263 264 264 # listes de dimensionnement (dependent de la geographie) 265 265 #----------------------------- … … 291 291 Dim_gsearise05 = paradim-gsearise05_mod.o geography-gsearise05_mod.o 292 292 293 # Liste_BLAS = band.o reduc.o \ 294 # BLAS/isamax.o BLAS/lsame.o BLAS/scopy.o BLAS/sgbmv.o BLAS/sgemm.o \ 295 # BLAS/sgemv.o BLAS/sger.o BLAS/sscal.o BLAS/sswap.o BLAS/stbsv.o \ 296 # BLAS/strsm.o BLAS/xerbla.o 297 ifeq ($(mkl_c), 0) 293 298 Liste_BLAS = band.o reduc.o \ 294 BLAS/isamax.o BLAS/lsame.o BLAS/scopy.o BLAS/sgbmv.o BLAS/sgemm.o \ 295 BLAS/sgemv.o BLAS/sger.o BLAS/sscal.o BLAS/sswap.o BLAS/stbsv.o \ 299 BLAS/isamax.o BLAS/lsame.o BLAS/scopy.o \ 300 BLAS/slamch.o BLAS/slabad.o BLAS/slange.o \ 301 BLAS/slassq.o BLAS/sgeqrf.o BLAS/sormqr.o \ 302 BLAS/strmv.o BLAS/slarft.o \ 303 BLAS/ilaslr.o BLAS/ilaslc.o BLAS/strmm.o \ 304 BLAS/slarfb.o BLAS/slapy2.o BLAS/slarfg.o BLAS/slarf.o \ 305 BLAS/sgeqr2.o BLAS/sorm2r.o BLAS/snrm2.o \ 306 BLAS/slaqps.o BLAS/slaqp2.o BLAS/sgeqp3.o \ 307 BLAS/slaset.o BLAS/slaisnan.o BLAS/sisnan.o \ 308 BLAS/slascl.o BLAS/sdot.o BLAS/slaic1.o \ 309 BLAS/slarzb.o BLAS/slarzt.o \ 310 BLAS/saxpy.o BLAS/slarz.o BLAS/sormr3.o \ 311 BLAS/sormrz.o BLAS/slatrz.o BLAS/stzrzf.o BLAS/sgelsy.o \ 312 BLAS/sgbmv.o BLAS/sgemm.o BLAS/sgemv.o BLAS/sger.o \ 313 BLAS/sscal.o BLAS/sswap.o BLAS/stbsv.o \ 296 314 BLAS/strsm.o BLAS/xerbla.o 297 315 316 endif 298 317 299 318 # ancienne methode … … 611 630 # nouveaux modules climat C. Dumas Fev 2015 612 631 climat_forcage_mois_mod.o : climat_forcage_mois_mod.f90 613 $(F T) climat_forcage_mois_mod.f90632 $(F_NETCDF) climat_forcage_mois_mod.f90 614 633 615 634 climat-forcage-insolation_mod_oneway.o : climat-forcage-insolation_mod_oneway.f90 616 $(F T) climat-forcage-insolation_mod_oneway.f90635 $(F_NETCDF) climat-forcage-insolation_mod_oneway.f90 617 636 618 637 climat-forcage-insolation_mod.o : climat-forcage-insolation_mod.f90 619 $(F T) climat-forcage-insolation_mod.f90638 $(F_NETCDF) climat-forcage-insolation_mod.f90 620 639 621 640 ablation_mod.o : ablation_mod.f90 … … 631 650 632 651 climat_GrIce2sea_years_mod.o : GrIce2sea_files/climat_GrIce2sea_years_mod.f90 633 $(F T) GrIce2sea_files/climat_GrIce2sea_years_mod.f90652 $(F_NETCDF) GrIce2sea_files/climat_GrIce2sea_years_mod.f90 634 653 635 654 climat-forcage-stat-mois_mod-0.1.o : climat-forcage-stat-mois_mod-0.1.f90 … … 1360 1379 # Routines venant de bibliotheques exterieures 1361 1380 reduc.o : reduc.f#module3d.mod 1362 $(FT) -w90reduc.f1381 $(FT) $(F_90) reduc.f 1363 1382 1364 1383 # band vient de Lapack 1365 1384 band.o: band.f 1366 $(FT) -w90 band.f 1385 $(FT) $(F_90) band.f 1386 1387 ifeq ($(mkl_c), 0) 1388 #BLAS 1389 isamax.o : BLAS/isamax.f 1390 $(FT) BLAS/isamax.f 1391 lsame.o : BLAS/lsame.f 1392 $(FT) BLAS/lsame.f 1393 scopy.o: BLAS/scopy.f 1394 $(FT) BLAS/scopy.f 1395 slamch.o: BLAS/slamch.f 1396 $(FT) BLAS/slamch.f 1397 slabad.o: BLAS/slabad.f 1398 $(FT) BLAS/slabad.f 1399 slange.o: BLAS/slange.f 1400 $(FT) BLAS/slange.f 1401 slassq.o: BLAS/slassq.f 1402 $(FT) BLAS/slassq.f 1403 sgeqrf.o: BLAS/sgeqrf.f 1404 $(FT) BLAS/sgeqrf.f 1405 sormqr.o: BLAS/sormqr.f 1406 $(FT) BLAS/sormqr.f 1407 strmv.o: BLAS/strmv.f 1408 $(FT) BLAS/strmv.f 1409 slarft.o: BLAS/slarft.f 1410 $(FT) BLAS/slarft.f 1411 ilaslr.o: BLAS/ilaslr.f 1412 $(FT) BLAS/ilaslr.f 1413 ilaslc.o: BLAS/ilaslc.f 1414 $(FT) BLAS/ilaslc.f 1415 strmm.o: BLAS/strmm.f 1416 $(FT) BLAS/strmm.f 1417 slarfb.o: BLAS/slarfb.f 1418 $(FT) BLAS/slarfb.f 1419 slapy2.o: BLAS/slapy2.f 1420 $(FT) BLAS/slapy2.f 1421 slarfg.o: BLAS/slarfg.f 1422 $(FT) BLAS/slarfg.f 1423 slarf.o: BLAS/slarf.f 1424 $(FT) BLAS/slarf.f 1425 sgeqr2.o: BLAS/sgeqr2.f 1426 $(FT) BLAS/sgeqr2.f 1427 saxpy.o: BLAS/saxpy.f 1428 $(FT) BLAS/saxpy.f 1429 slarz.o: BLAS/slarz.f 1430 $(FT) BLAS/slarz.f 1431 sorm2r.o: BLAS/sorm2r.f 1432 $(FT) BLAS/sorm2r.f 1433 snrm2.o: BLAS/snrm2.f 1434 $(FT) BLAS/snrm2.f 1435 slaqps.o: BLAS/slaqps.f 1436 $(FT) BLAS/slaqps.f 1437 slaqp2.o: BLAS/slaqp2.f 1438 $(FT) BLAS/slaqp2.f 1439 sgeqp3.o: BLAS/sgeqp3.f 1440 $(FT) BLAS/sgeqp3.f 1441 slaset.o: BLAS/slaset.f 1442 $(FT) BLAS/slaset.f 1443 slaisnan.o: BLAS/slaisnan.f 1444 $(FT) BLAS/slaisnan.f 1445 sisnan.o: BLAS/sisnan.f 1446 $(FT) BLAS/sisnan.f 1447 slascl.o: BLAS/slascl.f 1448 $(FT) BLAS/slascl.f 1449 sdot.o: BLAS/sdot.f 1450 $(FT) BLAS/sdot.f 1451 slaic1.o: BLAS/slaic1.f 1452 $(FT) BLAS/slaic1.f 1453 slarzb.o: BLAS/slarzb.f 1454 $(FT) BLAS/slarzb.f 1455 slarzt.o: BLAS/slarzt.f 1456 $(FT) BLAs/slarzt.f 1457 sormr3.o: BLAS/sormr3.f 1458 $(FT) BLAS/sormr3.f 1459 sormrz.o: BLAS/sormrz.f 1460 $(FT) BLAS/sormrz.f 1461 slatrz.o: BLAS/slatrz.f 1462 $(FT) BLAS/slatrz.f 1463 stzrzf.o: BLAS/stzrzf.f 1464 $(FT) BLAS/stzrzf.f 1465 sgelsy.o: BLAS/sgelsy.f 1466 $(FT) BLAS/sgelsy.f 1467 sgbmv.o: BLAS/sgbmv.f 1468 $(FT) BLAS/sgbmv.f 1469 sgemm.o: BLAS/sgemm.f 1470 $(FT) BLAS/sgemm.f 1471 sgemv.o: BLAS/sgemv.f 1472 $(FT) BLAS/sgemv.f 1473 sger.o: BLAS/sger.f 1474 $(FT) BLAS/sger.f 1475 sscal.o: BLAS/sscal.f 1476 $(FT) BLAS/sscal.f 1477 sswap.o: BLAS/sswap.f 1478 $(FT) BLAS/sswap.f 1479 stbsv.o: BLAS/stbsv.f 1480 $(FT) BLAS/stbsv.f 1481 strsm.o: BLAS/strsm.f 1482 $(FT) BLAS/strsm.f 1483 xerbla.o: BLAS/xerbla.f 1484 $(FT) BLAS/xerbla.f 1485 1486 endif 1367 1487 1368 1488 # routines netcdf … … 1419 1539 1420 1540 toy_retreat_mod.o : Recul_force_grounding_line/toy_retreat_mod.f90 1421 $(F T) Recul_force_grounding_line/toy_retreat_mod.f901541 $(F_NETCDF) Recul_force_grounding_line/toy_retreat_mod.f90 1422 1542 1423 1543 … … 1546 1666 $(Liste_Netcdf) \ 1547 1667 $(routines_communes) steps_time_loop.o \ 1548 $(routine_elliptiques) 1668 $(routine_elliptiques) \ 1669 $(Liste_BLAS) 1549 1670 1550 1671 $(LK) -o ../bin/Hemin-40 \ … … 1559 1680 $(Liste_Netcdf) \ 1560 1681 $(routines_communes) steps_time_loop.o \ 1561 $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) 1682 $(routine_elliptiques) $(NCDF_LIB) $(MKL_LIB) $(Liste_BLAS) 1562 1683 1563 1684 Hemin-15 : $(Dim_hemin15) $(mod_dim_communs) \ … … 1695 1816 1696 1817 clean : 1697 rm -f *.o *.mod *~ */*~ 1818 rm -f *.o *.mod *~ */*~ BLAS/*.o 1698 1819 1699 1820 -
trunk/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90
r19 r22 1030 1030 ! pour les classe 3, sortir les champs 2D necessaires 1031 1031 1032 if (cvar_nc(l) .eq. 3 .and. fait .eq . .FALSE.) then1032 if (cvar_nc(l) .eq. 3 .and. fait .eqv. .FALSE.) then 1033 1033 1034 1034 boucle_var2: do p=1,ntab -
trunk/SOURCES/Recul_force_grounding_line/toy_retreat_mod.f90
r4 r22 809 809 subroutine calc_eps_max 810 810 811 use declar_toy_retreat 811 use declar_toy_retreat, only: 812 812 implicit none 813 813 real :: gamma !< coefficient de flottaison -
trunk/SOURCES/climat_forcage_mois_mod.f90
r18 r22 153 153 namelist/snap_forcage_mois/filtr_t1,filtr_p1 154 154 155 ! dmr moved up since namelist declaration should be in the header 156 namelist/lapse_month/lect_lapse 157 155 158 ! formats pour les ecritures dans 42 156 159 428 format(A) … … 166 169 167 170 ! lapserate 168 namelist/lapse_month/lect_lapse169 171 rewind(num_param) ! pour revenir au debut du fichier param_list.dat 170 172 read(num_param,lapse_month) -
trunk/SOURCES/initial-phy-2.f90
r10 r22 140 140 read(num_param,timesteps) 141 141 142 namelist/timesteps/dtmin,dtmax,dtt,testdiag,tbegin,tend142 ! dmr&aurel ## namelist/timesteps/dtmin,dtmax,dtt,testdiag,tbegin,tend 143 143 write(num_rep_42,428)'!___________________________________________________________' 144 144 write(num_rep_42,428) '×teps ! bloc time steps' -
trunk/SOURCES/minmax-format.f90
r4 r22 25 25 character(len=8) :: ff 26 26 character(len=3) :: f1 27 ! dmr&aurel for gfortran compilation ... 28 character(len=1) :: dumf1 27 29 28 30 integer :: if1,if2,ipos,ipos2,n1,if3 … … 46 48 47 49 f1=ff(2:ipos-1) ! sous string contenant le premier chiffre 48 if1=inum(f1) ! valeur numerique du premier chiffre 50 51 ! dmr&aurel inum is not a std intrinsic 52 dumf1 = f1(1:2) 53 if1 = ichar(dumf1)-ichar('0') 54 ! dmr&aurel if1=inum(f1) ! valeur numerique du premier chiffre 49 55 50 56 print*, 'if1=',if1 … … 52 58 print*, 'f1= ',trim(f1)//'***' 53 59 f1=adjustr(f1) 54 if2=inum(f1) ! valeur numerique 60 ! dmr&aurel inum is not a std intrinsic 61 dumf1 = f1(1:2) 62 if2 = ichar(dumf1)-ichar('0') 63 ! if2=inum(f1) ! valeur numerique 55 64 print*, 'if2=',if2 56 65 … … 76 85 ipos2=index(ff,' ') 77 86 f1=ff(2:ipos2-1) ! sous string contenant le chiffre 78 if1=inum(f1) ! valeur numerique du premier chiffre 87 ! dmr&aurel inum is not a std intrinsic 88 dumf1 = f1(1:2) 89 if1 = ichar(dumf1)-ichar('0') 90 ! if1=inum(f1) ! valeur numerique du premier chiffre 79 91 n1=if1-1 80 92 valmin=int(-(10**n1))+1 -
trunk/SOURCES/spinup_mod.f90
r4 r22 209 209 coef_defmx(i,j) = 1. 210 210 Uxbar(i,j) = Vcol_x(i,j) 211 flgzmx(i,j) = Vcol_x(i,j)211 flgzmx(i,j) = transfer(Vcol_x(i,j),flgzmx(i,j)) 212 212 uxdef(i,j) = 0. 213 213 Ubx(i,j) = Vcol_x(i,j) … … 250 250 coef_defmy(i,j) = 1. 251 251 Uybar(i,j) = Vcol_y(i,j) 252 flgzmy(i,j) = Vcol_y(i,j)252 flgzmy(i,j) = transfer(Vcol_y(i,j),flgzmy(i,j)) 253 253 uydef(i,j) = 0. 254 254 Uby(i,j) = Vcol_y(i,j)
Note: See TracChangeset
for help on using the changeset viewer.