Changeset 23989

Timestamp:

06/05/19 10:23:04 (6 years ago)

Author:

dmenemen

Message:

attempting to run coupled MITgcm/ISSM (step 11)
crashes on initialization

File:

• issm/trunk-jpl/test/NightlyRun/test4003.m (modified) (2 diffs)

issm/trunk-jpl/test/NightlyRun/test4003.m

@@ -345 +345 @@
                   md.initialization.pressure=md.results.TransientSolution(end).Pressure;
                   md.transient.isoceancoupling=1;
-                  md.transient.isgroundingline=1;
+                  md.transient.isgroundingline=0;
                   md.masstransport.requested_outputs={'default','BasalforcingsFloatingiceMeltingRate'};
 
@@ -359 +359 @@
         !\mkdir run
         !\cp ../MITgcm/build/mitgcmuv run
-        !\cp ../MITgcm/input/* run
-
-        % initial salinity
-        S=iniSalt*ones(Nx,Ny,Nz);
-        writebin('run/Salt.bin',S);
-
-        % initial temperature
-        T=iniTheta*ones(Nx,Ny,Nz);
-        writebin('run/Theta.bin',T);
-
-        % initial velocity
-        Z=zeros(Nx,Ny,Nz);
-        writebin('run/Uvel.bin',Z);
-        writebin('run/Vvel.bin',Z);
-
-        % initial sea surface height
-        Z=zeros(Nx,Ny);
-        writebin('run/Etan.bin',Z);
-
-        % salinity boundary conditions
-        S=obcSalt*ones(Ny,Nz);
-        thk=delZ*ones(Nz,1);
-        bot=cumsum(thk);
-        ik=find(bot<=mlDepth);
-        S(:,ik)=mlSalt;
-        writebin('run/OBs.bin',S);
-
-        % temperature boundary conditions
-        T=obcTheta*ones(Ny,Nz);
-        T(:,ik)=mlTheta;
-        writebin('run/OBt.bin',T);
-
-        % zonal velocity boundary conditions
-        U=obcUvel*ones(Ny,Nz);
-        writebin('run/OBu.bin',U);
-
-        % zero boundary conditions
-        Z=zeros(Ny,Nz);
-        writebin('run/zeros.bin',Z);
-
-        % build parameter file data.obcs
-        fidi=fopen('../MITgcm/input/data.obcs','r');
-        fido=fopen('run/data.obcs','w');
-        tline = fgetl(fidi);
-        fprintf(fido,'%s\n',tline);
-        while 1
-            tline = fgetl(fidi);
-            if ~ischar(tline), break, end
-            %do the change here: 
-            if strcmpi(tline,' OB_Iwest = 40*1,'),
-                fprintf(fido,'%s%i%s\n',' OB_Iwest = ',Ny,'*1,');
-                continue;
-            end
-            if strcmpi(tline,' OB_Ieast = 40*-1,'),
-                fprintf(fido,'%s%i%s\n',' OB_Ieast = ',Ny,'*-1,');
-                continue;
-            end
-            fprintf(fido,'%s\n',tline);
-        end
-        %close  files
-        fclose(fidi);
-        fclose(fido);
-
-        %save bathymetry and bedrock in run directory
-        writebin('run/bathymetry.bin',bathymetry);
-        writebin('run/icetopo.bin',draft);
+        !\cp ../MITgcm/install_dngoldberg/verification/shelfice_remeshing/input/* run
         % }}}
         % {{{ generate MITgcm parameter file data