Changeset 3562
- Timestamp:
- Dec 17, 2024, 4:19:04 PM (5 days ago)
- Location:
- trunk/LMDZ.GENERIC/libf/phystd
- Files:
-
- 1 added
- 3 edited
-
dyn1d/rcm1d.F (modified) (13 diffs)
-
dyn1d/restart1D_mod.F90 (added)
-
iostart.F90 (modified) (1 diff)
-
physiq_mod.F90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F
r3515 r3562 38 38 use phys_state_var_mod, only: phys_state_var_init 39 39 use physiq_mod, only: physiq 40 use restart1D_mod, only: writerestart1D 41 use iostart, only: length 42 use netcdf, only: NF90_OPEN, NF90_NOERR, NF90_NOWRITE, 43 & nf90_strerror,NF90_INQ_VARID, NF90_GET_VAR, 44 & NF90_CLOSE 40 45 implicit none 41 46 … … 78 83 INTEGER ilayer,ilevel,isoil,idt,iq 79 84 LOGICAl firstcall,lastcall 85 LOGICAL restart 86 INTEGER nid_restart1D, nid_restartfi 87 INTEGER controleid, uid, vid, tempid, tsurfid, did, tid 88 INTEGER,ALLOCATABLE :: qid(:), qsurfid(:) 80 89 c 81 90 INTEGER day0 ! date initial (sol ; =0 a Ls=0) … … 122 131 123 132 logical oldcompare, earthhack,saveprofile 133 real controle1D(length) 124 134 125 135 ! added by RW for zlay computation … … 186 196 call system("echo 'INCLUDEDEF=rcm1d.def' >> run.def") 187 197 endif 198 199 ! Check restart 200 call getin("restart",restart) 201 if (restart) then 202 ierr = NF90_OPEN('start1D.nc', NF90_NOWRITE, nid_restart1D) 203 if (ierr.NE.NF90_NOERR) then 204 write(*,*)'readstart1D: problem opening 1D start file' 205 write(*,*)trim(nf90_strerror(ierr)) 206 call abort_physic("readstart1D","Cannot open file",1) 207 endif 208 ierr = NF90_OPEN('startfi.nc', NF90_NOWRITE, nid_restartfi) 209 if (ierr.NE.NF90_NOERR) then 210 write(*,*)'readstart1D: problem opening startfi file' 211 write(*,*)trim(nf90_strerror(ierr)) 212 call abort_physic("readstart1D","Cannot open file",1) 213 endif 214 endif ! restart 188 215 189 216 ! read nq from traceur.def … … 474 501 ENDIF 475 502 476 psurf = -99999. 477 PRINT *,'SURFACE PRESSURE in Pa ?' 478 call getin("psurf",psurf) 479 IF (psurf.eq.-99999.) THEN 480 PRINT *,"STOP. I NEED psurf IN RCM1D.DEF." 481 STOP 482 ELSE 483 PRINT *,"--> psurf = ",psurf 484 ENDIF 503 if (restart) then 504 ierr=NF90_INQ_VARID(nid_restart1D,'controle',controleid) 505 if (ierr==NF90_NOERR) then 506 ierr=NF90_GET_VAR(nid_restart1D,controleid,controle1D) 507 if (ierr/=NF90_NOERR) then 508 PRINT*, 'read restart 1D: Failed loading psurf' 509 CALL abort_physic("readstart1D","Failed to read variable",1) 510 endif 511 endif 512 psurf = controle1D(1) 513 PRINT *,"In restart1D.nc, day0=",controle1D(2), 514 & " and time=",controle1D(3) 515 PRINT*,"(These values are NOT used, ", 516 & "those given in rcm1d.def are used)" 517 else ! restart 518 psurf = -99999. 519 PRINT *,'SURFACE PRESSURE in Pa ?' 520 call getin("psurf",psurf) 521 IF (psurf.eq.-99999.) THEN 522 PRINT *,"STOP. I NEED psurf IN RCM1D.DEF." 523 STOP 524 ELSE 525 PRINT *,"--> psurf = ",psurf 526 ENDIF 527 endif ! restart 485 528 !! we need reference pressures 486 529 pa=psurf/30. … … 494 537 c ------------------------------------ 495 538 c Date (en sols depuis le solstice de printemps) du debut du run 496 day0 = 0 ! default value for day0 497 write(*,*) 'Initial date (in martian sols ; =0 at Ls=0)?' 498 call getin("day0",day0) 499 day=float(day0) 500 write(*,*) " day0 = ",day0 539 !if (restart) then 540 ! ! day 541 ! ierr=NF90_INQ_VARID(nid_restart1D,'day',did) 542 ! if (ierr==NF90_NOERR) then 543 ! ierr=NF90_GET_VAR(nid_restart1D,did,day) 544 ! if (ierr/=NF90_NOERR) then 545 ! PRINT*, 'read restart 1D: Failed loading day' 546 ! CALL abort_physic("readstart1D","Failed to read variable",1) 547 ! endif 548 ! endif 549 550 ! ! time 551 ! ierr=NF90_INQ_VARID(nid_restart1D,'time',tid) 552 ! if (ierr==NF90_NOERR) then 553 ! ierr=NF90_GET_VAR(nid_restart1D,tid,time) 554 ! if (ierr/=NF90_NOERR) then 555 ! PRINT*, 'read restart 1D: Failed loading time' 556 ! CALL abort_physic("readstart1D","Failed to read variable",1) 557 ! endif 558 ! endif 559 560 !else 561 day0 = 0 ! default value for day0 562 write(*,*) 'Initial date (in martian sols ; =0 at Ls=0)?' 563 call getin("day0",day0) 564 day=float(day0) 565 write(*,*) " day0 = ",day0 501 566 c Heure de demarrage 502 time=0 ! default value for time503 write(*,*)'Initial local time (in hours, between 0 and 24)?'504 call getin("time",time)505 write(*,*)" time = ",time506 time=time/24.E+0 ! convert time (hours) to fraction of sol507 567 time=0 ! default value for time 568 write(*,*)'Initial local time (in hours, between 0 and 24)?' 569 call getin("time",time) 570 write(*,*)" time = ",time 571 time=time/24.E+0 ! convert time (hours) to fraction of sol 572 !endif 508 573 509 574 c Discretisation (Definition de la grille et des pas de temps) … … 629 694 630 695 do iq=1,nq 696 697 if (restart) then 698 if (iq.eq.1) then 699 allocate(qid(nq)) 700 allocate(qsurfid(nq)) 701 endif 702 ! read q 703 ierr=NF90_INQ_VARID(nid_restart1D,noms(iq),qid(iq)) 704 if (ierr==NF90_NOERR) then 705 ierr=NF90_GET_VAR(nid_restart1D,qid(iq),q(:,iq), 706 & count=(/1,nlayer,1,1/)) 707 if (ierr/=NF90_NOERR) then 708 PRINT*, 'read restart 1D: Failed loading ', noms(iq) 709 write(*,*)trim(nf90_strerror(ierr)) 710 CALL abort_physic("readstart1D", 711 & "Failed to read variable",1) 712 endif 713 endif 714 ! read qsurf 715 ierr=NF90_INQ_VARID(nid_restartfi,noms(iq),qsurfid(iq)) 716 if (ierr==NF90_NOERR) then 717 ierr=NF90_GET_VAR(nid_restartfi,qsurfid(iq),qsurf(iq)) 718 if (ierr/=NF90_NOERR) then 719 PRINT*, 'read restartfi: Failed loading ', noms(iq) 720 write(*,*)trim(nf90_strerror(ierr)) 721 CALL abort_physic("readstartfi", 722 & "Failed to read variable",1) 723 endif 724 endif 725 else ! restart 631 726 632 txt=""633 write(txt,"(a)") tname(iq)634 write(*,*)" tracer:",trim(txt)727 txt="" 728 write(txt,"(a)") tname(iq) 729 write(*,*)" tracer:",trim(txt) 635 730 636 ! CO2637 if (txt.eq."co2_ice") then638 q(:,iq)=0. ! kg/kg of atmosphere639 qsurf(iq)=0. ! kg/m2 at the surface640 ! Look for a "profile_co2_ice" input file641 open(91,file='profile_co2_ice',status='old',642 & form='formatted',iostat=ierr)643 if (ierr.eq.0) then644 read(91,*) qsurf(iq)645 do ilayer=1,nlayer646 read(91,*) q(ilayer,iq)647 enddo648 else649 write(*,*) "No profile_co2_ice file!"650 endif651 close(91)652 endif ! of if (txt.eq."co2")731 ! CO2 732 if (txt.eq."co2_ice") then 733 q(:,iq)=0. ! kg/kg of atmosphere 734 qsurf(iq)=0. ! kg/m2 at the surface 735 ! Look for a "profile_co2_ice" input file 736 open(91,file='profile_co2_ice',status='old', 737 & form='formatted',iostat=ierr) 738 if (ierr.eq.0) then 739 read(91,*) qsurf(iq) 740 do ilayer=1,nlayer 741 read(91,*) q(ilayer,iq) 742 enddo 743 else 744 write(*,*) "No profile_co2_ice file!" 745 endif 746 close(91) 747 endif ! of if (txt.eq."co2") 653 748 654 ! WATER VAPOUR655 if (txt.eq."h2o_vap") then656 q(:,iq)=0. ! kg/kg of atmosphere657 qsurf(iq)=0. ! kg/m2 at the surface658 ! Look for a "profile_h2o_vap" input file659 open(91,file='profile_h2o_vap',status='old',660 & form='formatted',iostat=ierr)661 if (ierr.eq.0) then662 read(91,*) qsurf(iq)663 do ilayer=1,nlayer664 read(91,*) q(ilayer,iq)665 enddo666 else667 write(*,*) "No profile_h2o_vap file!"668 endif669 close(91)670 endif ! of if (txt.eq."h2o_vap")749 ! WATER VAPOUR 750 if (txt.eq."h2o_vap") then 751 q(:,iq)=0. ! kg/kg of atmosphere 752 qsurf(iq)=0. ! kg/m2 at the surface 753 ! Look for a "profile_h2o_vap" input file 754 open(91,file='profile_h2o_vap',status='old', 755 & form='formatted',iostat=ierr) 756 if (ierr.eq.0) then 757 read(91,*) qsurf(iq) 758 do ilayer=1,nlayer 759 read(91,*) q(ilayer,iq) 760 enddo 761 else 762 write(*,*) "No profile_h2o_vap file!" 763 endif 764 close(91) 765 endif ! of if (txt.eq."h2o_vap") 671 766 672 ! WATER ICE 673 if (txt.eq."h2o_ice") then 674 q(:,iq)=0. ! kg/kg of atmosphere 675 qsurf(iq)=0. ! kg/m2 at the surface 676 ! Look for a "profile_h2o_ice" input file 677 open(91,file='profile_h2o_ice',status='old', 678 & form='formatted',iostat=ierr) 679 if (ierr.eq.0) then 680 read(91,*) qsurf(iq) 681 do ilayer=1,nlayer 682 read(91,*) q(ilayer,iq) 683 enddo 684 else 685 write(*,*) "No profile_h2o_ice file!" 686 endif 687 close(91) 688 endif ! of if (txt.eq."h2o_ice") 689 690 !_vap 691 if((txt .ne. 'h2o_vap') .and. 692 & (index(txt,'_vap' ) .ne. 0)) then 693 q(:,iq)=0. !kg/kg of atmosphere 694 qsurf(iq) = 0. !kg/kg of atmosphere 695 ! Look for a "profile_...._vap" input file 696 tracer_profile_file_name="" 697 tracer_profile_file_name='profile_'//txt 698 open(91,file=tracer_profile_file_name,status='old', 699 & form="formatted",iostat=ierr) 700 if (ierr .eq. 0) then 701 read(91,*)qsurf(iq) 702 do ilayer=1,nlayer 703 read(91,*)q(ilayer,iq) 704 enddo 705 else 706 write(*,*) "No initial profile " 707 write(*,*) " for this tracer :" 708 write(*,*) txt 709 endif 710 close(91) 711 endif ! (txt .eq. "_vap") 712 !_ice 713 if((txt.ne."h2o_ice") .and. 714 & (index(txt,'_ice' ) /= 0)) then 715 q(:,iq)=0. !kg/kg of atmosphere 716 qsurf(iq) = 0. !kg/kg of atmosphere 717 endif ! we only initialize the solid at 0 718 enddo ! of do iq=1,nq 719 767 ! WATER ICE 768 if (txt.eq."h2o_ice") then 769 q(:,iq)=0. ! kg/kg of atmosphere 770 qsurf(iq)=0. ! kg/m2 at the surface 771 ! Look for a "profile_h2o_ice" input file 772 open(91,file='profile_h2o_ice',status='old', 773 & form='formatted',iostat=ierr) 774 if (ierr.eq.0) then 775 read(91,*) qsurf(iq) 776 do ilayer=1,nlayer 777 read(91,*) q(ilayer,iq) 778 enddo 779 else 780 write(*,*) "No profile_h2o_ice file!" 781 endif 782 close(91) 783 endif ! of if (txt.eq."h2o_ice") 784 785 !_vap 786 if((txt .ne. 'h2o_vap') .and. 787 & (index(txt,'_vap' ) .ne. 0)) then 788 q(:,iq)=0. !kg/kg of atmosphere 789 qsurf(iq) = 0. !kg/kg of atmosphere 790 ! Look for a "profile_...._vap" input file 791 tracer_profile_file_name="" 792 tracer_profile_file_name='profile_'//txt 793 open(91,file=tracer_profile_file_name,status='old', 794 & form="formatted",iostat=ierr) 795 if (ierr .eq. 0) then 796 read(91,*)qsurf(iq) 797 do ilayer=1,nlayer 798 read(91,*)q(ilayer,iq) 799 enddo 800 else 801 write(*,*) "No initial profile " 802 write(*,*) " for this tracer :" 803 write(*,*) txt 804 endif 805 close(91) 806 endif ! (txt .eq. "_vap") 807 !_ice 808 if((txt.ne."h2o_ice") .and. 809 & (index(txt,'_ice' ) /= 0)) then 810 q(:,iq)=0. !kg/kg of atmosphere 811 qsurf(iq) = 0. !kg/kg of atmosphere 812 endif ! we only initialize the solid at 0 813 endif ! restart 814 enddo ! of do iq=1,nq 815 720 816 endif ! of tracer 721 817 … … 739 835 w(1:nlayer)=0 740 836 741 c Initialisation des vents au premier pas de temps 742 DO ilayer=1,nlayer 743 u(ilayer)=gru 744 v(ilayer)=grv 745 ENDDO 837 if (restart) then 838 ierr=NF90_INQ_VARID(nid_restart1D,'u',uid) 839 if (ierr==NF90_NOERR) then 840 ierr=NF90_GET_VAR(nid_restart1D,uid,u, 841 & count=(/1,nlayer,1,1/)) 842 if (ierr/=NF90_NOERR) then 843 PRINT*, 'read restart 1D: Failed loading u' 844 write(*,*)trim(nf90_strerror(ierr)) 845 CALL abort_physic("readstart1D", 846 & "Failed to read variable",1) 847 endif 848 endif 849 ierr=NF90_INQ_VARID(nid_restart1D,'v',vid) 850 if (ierr==NF90_NOERR) then 851 ierr=NF90_GET_VAR(nid_restart1D,vid,v, 852 & count=(/1,nlayer,1,1/)) 853 if (ierr/=NF90_NOERR) then 854 PRINT*, 'read restart 1D: Failed loading v' 855 write(*,*)trim(nf90_strerror(ierr)) 856 CALL abort_physic("readstart1D", 857 & "Failed to read variable",1) 858 endif 859 endif 860 else ! restart 861 862 c Initialisation des vents au premier pas de temps 863 DO ilayer=1,nlayer 864 u(ilayer)=gru 865 v(ilayer)=grv 866 ENDDO 867 868 endif ! restart 746 869 747 870 c energie cinetique turbulente … … 859 982 ! endif 860 983 984 861 985 c profil de temperature au premier appel 862 986 c -------------------------------------- 863 987 pks=psurf**rcp 864 988 989 if (restart) then 990 ! read temp 991 ierr=NF90_IN Q_VARID(nid_restart1D,'temp',tempid) 992 if (ierr==NF90_NOERR) then 993 ierr=NF90_GET_VAR(nid_restart1D,tempid,temp, 994 & count=(/1,nlayer,1,1/)) 995 if (ierr/=NF90_NOERR) then 996 PRINT*, 'read restart 1D: Failed loading temp' 997 CALL abort_physic("readstart1D", 998 & "Failed to read variable",1) 999 endif 1000 endif 1001 ! read tsurf 1002 ierr=NF90_INQ_VARID(nid_restartfi,'tsurf',tsurfid) 1003 if (ierr==NF90_NOERR) then 1004 ierr=NF90_GET_VAR(nid_restartfi,tsurfid,tsurf) 1005 if (ierr/=NF90_NOERR) then 1006 PRINT*, 'read restartfi: Failed loading tsurf' 1007 CALL abort_physic("readstartfi", 1008 & "Failed to read variable",1) 1009 endif 1010 endif 1011 else ! restart 1012 865 1013 c altitude en km dans profile: on divise zlay par 1000 866 tmp1(0)=0.E+0 867 DO ilayer=1,nlayer 868 tmp1(ilayer)=zlay(ilayer)/1000.E+0 869 ENDDO 870 call profile(nlayer+1,tmp1,tmp2) 871 872 tsurf(1)=tmp2(0) 873 DO ilayer=1,nlayer 874 temp(ilayer)=tmp2(ilayer) 875 ENDDO 876 print*,"check" 877 PRINT*,"INPUT SURFACE TEMPERATURE",tsurf(1) 878 PRINT*,"INPUT TEMPERATURE PROFILE",temp 1014 tmp1(0)=0.E+0 1015 DO ilayer=1,nlayer 1016 tmp1(ilayer)=zlay(ilayer)/1000.E+0 1017 ENDDO 1018 call profile(nlayer+1,tmp1,tmp2) 1019 1020 tsurf(1)=tmp2(0) 1021 DO ilayer=1,nlayer 1022 temp(ilayer)=tmp2(ilayer) 1023 ENDDO 1024 print*,"check" 1025 PRINT*,"INPUT SURFACE TEMPERATURE",tsurf(1) 1026 PRINT*,"INPUT TEMPERATURE PROFILE",temp 1027 1028 endif ! restart 879 1029 880 1030 c Initialisation albedo / inertie du sol … … 937 1087 ! ----------------- 938 1088 #ifndef MESOSCALE 939 if(tracer.and.photochem) then 1089 if(tracer.and.photochem.and. .not.restart) then 1090 print*, "Calling inichim_1D" 940 1091 call initracer(1,nq) 941 1092 allocate(nametmp(nq)) … … 947 1098 #endif 948 1099 1100 if (restart) then 1101 ierr = NF90_CLOSE(nid_restart1D) 1102 if (ierr/=NF90_NOERR) then 1103 PRINT*, 'read restart1D: Failed closing restart1D.nc' 1104 CALL abort_physic("readstart1D", 1105 & "Failed closing file",1) 1106 endif 1107 ierr = NF90_CLOSE(nid_restartfi) 1108 if (ierr/=NF90_NOERR) then 1109 PRINT*, 'read restartfi: Failed closing restartfi.nc' 1110 CALL abort_physic("readstartfi", 1111 & "Failed closing file",1) 1112 endif 1113 endif 949 1114 950 1115 c Write a "startfi" file … … 953 1118 c It is needed to transfert physics variables to "physiq"... 954 1119 955 call physdem0("startfi.nc",longitude,latitude,nsoilmx,1,llm,nq, 956 & dtphys,real(day0),time,cell_area, 957 & albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe) 958 call physdem1("startfi.nc",nsoilmx,1,llm,nq, 959 & dtphys,time, 960 & tsurf,tsoil,emis,albedo,q2,qsurf, 961 & cloudfrac,totcloudfrac,hice, 962 & rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice) 963 1120 if (.not. restart) then 1121 call physdem0("startfi.nc",longitude,latitude,nsoilmx,1,llm,nq, 1122 & dtphys,real(day0),time,cell_area, 1123 & albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe) 1124 call physdem1("startfi.nc",nsoilmx,1,llm,nq, 1125 & dtphys,time, 1126 & tsurf,tsoil,emis,albedo,q2,qsurf, 1127 & cloudfrac,totcloudfrac,hice, 1128 & rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice) 1129 endif 964 1130 c======================================================================= 965 1131 c BOUCLE TEMPORELLE DU MODELE 1D … … 1098 1264 endif 1099 1265 1266 c Produce the restart1D.nc file (MM) 1267 if (lastcall) then 1268 call writerestart1D('restart1D.nc',nlayer,nsoil,day,time,psurf, 1269 & temp,tsoil,u,v,nq,q) 1270 endif 1100 1271 1101 1272 ENDDO ! fin de la boucle temporelle -
trunk/LMDZ.GENERIC/libf/phystd/iostart.F90
r3552 r3562 1140 1140 END SUBROUTINE put_var_c1 1141 1141 1142 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1143 1144 SUBROUTINE create_restart1D(filename,nid_restart) 1145 USE netcdf, only: NF90_CREATE, NF90_CLOBBER, NF90_64BIT_OFFSET, & 1146 NF90_NOERR, nf90_strerror, & 1147 NF90_PUT_ATT, NF90_GLOBAL, NF90_DEF_DIM, & 1148 NF90_UNLIMITED, NF90_ENDDEF 1149 USE mod_phys_lmdz_para, only: is_master 1150 USE mod_grid_phy_lmdz, only: klon_glo 1151 USE dimphy, only: klev, klevp1 1152 USE tracer_h, only: nqtot 1153 USE comsoil_h, only: nsoilmx 1154 ! USE slab_ice_h, only: noceanmx 1155 USE ocean_slab_mod, ONLY: nslay 1156 1157 IMPLICIT NONE 1158 CHARACTER(LEN=*),INTENT(IN) :: filename 1159 INTEGER,INTENT(INOUT) :: nid_restart 1160 INTEGER :: ierr 1161 1162 IF (is_master) THEN 1163 1164 ierr=NF90_CREATE(filename,IOR(NF90_CLOBBER,NF90_64BIT_OFFSET), & 1165 nid_restart) 1166 IF (ierr/=NF90_NOERR) THEN 1167 write(*,*)'create_restart1D: problem creating file '//trim(filename) 1168 write(*,*)trim(nf90_strerror(ierr)) 1169 CALL abort_physic("create_restart1D","Failed creating file",1) 1170 ENDIF 1171 1172 ierr=NF90_PUT_ATT(nid_restart,NF90_GLOBAL,"title",& 1173 "Physics start file") 1174 IF (ierr/=NF90_NOERR) THEN 1175 write(*,*)'create_restart1D: problem writing title ' 1176 write(*,*)trim(nf90_strerror(ierr)) 1177 ENDIF 1178 1179 ierr=NF90_DEF_DIM(nid_restart,"physical_points",klon_glo,idim2) 1180 IF (ierr/=NF90_NOERR) THEN 1181 write(*,*)'create_restart1D: problem defining physical_points dimension ' 1182 write(*,*)trim(nf90_strerror(ierr)) 1183 CALL abort_physic("create_restart1D","Failed defining physical_points",1) 1184 ENDIF 1185 1186 ierr=NF90_DEF_DIM(nid_restart,"nlayer",klev,idim6) 1187 IF (ierr/=NF90_NOERR) THEN 1188 write(*,*)'create_restart1D: problem defining nlayer dimension ' 1189 write(*,*)trim(nf90_strerror(ierr)) 1190 CALL abort_physic("create_restart1D","Failed defining nlayer",1) 1191 ENDIF 1192 1193 ierr=NF90_ENDDEF(nid_restart) 1194 IF (ierr/=NF90_NOERR) THEN 1195 write(*,*)'create_restart1D: problem ending definition mode ' 1196 write(*,*)trim(nf90_strerror(ierr)) 1197 CALL abort_physic("create_restart1D","Failed ending definition mode",1) 1198 ENDIF 1199 ENDIF 1200 1201 END SUBROUTINE create_restart1D 1202 1142 1203 END MODULE iostart -
trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90
r3522 r3562 792 792 793 793 #ifndef MESOSCALE 794 if (ngrid.ne.1) then ! Note : no need to create a restart file in 1d.794 !if (ngrid.ne.1) then ! Note : no need to create a restart file in 1d. 795 795 call physdem0("restartfi.nc",longitude,latitude,nsoilmx,ngrid,nlayer,nq, & 796 796 ptimestep,pday+nday,time_phys,cell_area, & 797 797 albedo_bareground,inertiedat,zmea,zstd,zsig,zgam,zthe) 798 endif798 !endif 799 799 800 800 #endif … … 2397 2397 #ifndef MESOSCALE 2398 2398 2399 if (ngrid.ne.1) then2399 !if (ngrid.ne.1) then 2400 2400 write(*,*)'PHYSIQ: for physdem ztime_fin =',ztime_fin 2401 2401 … … 2411 2411 cloudfrac,totcloudfrac,hice, & 2412 2412 rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice) 2413 endif2413 !endif 2414 2414 #endif 2415 2415 if(ok_slab_ocean) then
Note: See TracChangeset
for help on using the changeset viewer.
